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MEASUREMENTS WITH A WHOLE BODY COUNTER! 


K. G. McNEILL AND R. M. GREEN 


ABSTRACT 


An account is given of work on the measurement of the gamma radioactivity 
of humans. A steel room has been constructed, in which this gamma activity may 
be measured to a statistical accuracy of a few per cent in 30 minutes using a 
5-in. Nal crystal as detector. The average potassium content of thirty 19-year- 
old males living in Ontario has been found to be 0.212% of body weight, and 
the average Cs/K ratio 51 wuc/g K. 


INTRODUCTION 


The increase in radioactivity of the world due to nuclear weapon tests, the 
subsequent ingestion of this activity by humans, and its possible effects on 
human health have in recent years increased the interest in measurements 
of the current radioactivity of the body and other biological material (Burch 
and Spiers 1953; Reines et a/. 1953; Mayneord et al. 1958). 

The gamma activity of the normal human body is largely due to K*, a 
naturally occurring isotope of potassium, and to Cs"? from nuclear fallout; 
1.46-Mev gamma rays are emitted in 11% of the K*® disintegrations, and the 
potassium of a typical adult emits 500 gamma rays of this energy per second. 
Cs'*7 decays to Ba!’ with a half-life of 27 years; 92% of these disintegrations 
are to the isomeric state of Ba!*? which quickly (2.6 minutes) de-excites with 
emission of a 0.66-Mev gamma ray. At the present time the gamma activity 
from this source is of the same order as that from K* in the body. 

Extensive measurements of gamma activity from the body have been made 
by the Argonne National Laboratories group (e.g. Miller and Marinelli 1956), 
at other places in the United States (e.g. Anderson et a/. 1957) and in Britain 
(Burch and Spiers 1953), but so far no extensive investigation has been 
reported on similar work in Canada. Results for different latitudes are, how- 
ever, of value owing to the geographical variation of fallout (Libby 1958). At 
the University of Toronto apparatus has been built which enables the gamma 
activity of a person to be determined in a relatively short time. The purpose 
of this paper is to report on this measurement equipment and to give some 
of the results so far obtained on the activity of persons in Ontario. 


1Manuscript received February 25, 1959. 
Contribution from the Department of Physics, University of Toronto, Toronto, Ontario. 
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THE DETECTION APPARATUS 

In order to measure the low level of activity from the body it is necessary 
to reduce other sources of radiation to as great an extent as possible, for 
otherwise background from cosmic rays, local gamma emitters, and airborne 
contamination will swamp the body activity. For example, using a 5-in. Nal 
crystal the maximum count in the photopeak of the K*® spectrum of a normal 
subject in our standard geometry is only 1/20th of the 1.46-Mev background 
count with no shielding. In the building of radiation shielding, we have 
followed very closely the work and designs of the Argonne National Labora- 
tories group (Miller e¢ al. 1956). A steel room has been constructed with 
internal dimensions 5 ft X6 ftX7 ft 4 in., and with 8 in. thick walls, floor, and 
ceiling. The total weight of the room is 50 tons. For convenience in handling 
during construction, the steel used is in the form of }-in. thick plates, which 
have dimensions ranging from 10 in. to 24 in. wide and from 4 ft 6in. to 
8 ft 9 in. long. The skeleton of the room is of steel girders and angle iron. The 
door, which weighs 3 tons, has such a width that when open a space of 17 
inches is available for entrance into the room. 

For detecting gamma radiation, sodium iodide crystals are used. During 
the initial stages of background measurements a 13-in. diameter X 2-in. 
crystal was employed, but when the background had been sufficiently attenu- 
ated by the increase in the thickness of the steel walls, a change was made 
to a 5-in. diameter X 4-in. crystal. The latter was made from especially 
potassium free iodide, and was mounted in an electrolytic copper can. The 
normal glass window was replaced by a quartz one. These changes, made by 
the Harshaw Chemical Company, served to reduce the inherent radioactivity 
of the crystal assembly. Du Mont 6292 and 6364 photomultiplier tubes are 
used, and the pulses from these fed by conventional electronics to a Radiation 
Instrument Development Laboratories 100-channel pulse height analyzer, 
model 3300. 

For body activity measurements, subjects are seated in a reclining chair 
made of angle iron and canvas webbing, which was found to have no appre- 
ciable activity. The crystal is hung above the subject’s abdomen in a standard 
position with respect to the chair, following the results of experiments made 
by the Argonne group (Miller 1956). 

The whole steel room is enveloped in a plastic bag to help to exclude 
dust-borne activity. Filtered air is blown into the room from the main labora- 
tory, the air of which is itself conditioned and filtered. 


BACKGROUND MEASUREMENTS 

In order that information could be obtained on the attenuation of the 
background by the steel walls of the counting room, measurements were 
taken with the 2-in. crystal, and, latterly, with the 5-in. crystal, as the wall 
thickness was increased. One such measurement was taken for every 1 inch 
of steel added to the walls. The results are illustrated in Table I. It was 
found that although the attenuation was approximately exponential for the 
first few inches, the addition of the seventh and eighth inches of steel made 
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only a small difference to the background count. This is because at this 
thickness the main contributions to the background are from sources in the 
crystal assembly and shielding itself. 


TABLE I 


Attenuation ratios 


No. of counts/minute in 2-in. Nal University of Toronto 
Inches corresponding to: attenuation ratio 
of —_—— --—— ———__——_——— - = 
steel E, > 0.2 Mev Ey, > 0.51 Mev E, > 0.2 Mev E, > 0.51 Mev 
0 2,980 952 1.0 1.0 
1 454 0.48 
2 234 0.25 
3 286 130 0.096 0.137 
+ 169 76 0.057 0.080 
5 114 58 0.038 0.061 
6 86 47 0.029 0.049 
7 82 45 0.028 0.047 
8 &1 45 0.027 0.046 


Peaks in the background spectrum obtained with no steel shielding may 
all be ascribed to gamma rays from RaC, except for the K*® peak at 1.46 
Mev. With maximum shielding, the most obvious peak is at 1.46 Mev, but 
there are also peaks at 0.72 Mev and at about 0.5 Mev. This latter peak 
was found in the gamma spectra of samples of steel, and is due to the fission 
products Ru! and Ru!®. This element is siderophilic (Beamish 1958) and it 
is presumably fission product ruthenium in the shielding steel that is the 
cause of the 0.5-Mev peak. The peak at 0.72 Mev is ascribed to the fission 
products Zr®-Nb® present in dust. As was expected (Miller and Steingraber 
1957), a subsidiary experiment showed that three-quarters of the residual 
Kx4° activity comes from the potassium in the glass of the photomultiplier. 

With the 8-in. steel shielding, the total background in the energy range 
200 kev — 2.0 Mev is 470 c.p.m., and in the range 1 — 2 Mev is 66¢.p.m. 


CALIBRATION OF THE APPARATUS 


The primary purpose of the apparatus is to measure the gamma radio- 
activity of the human body. In order that a particular count in the detector 
might be related to a definite amount of radioactivity in the body, it was 
necessary to calibrate the apparatus using known amounts of radioactive 
substance. The isotope K* was used for this work, for not only will it be 
distributed in the body with the normal potassium content, but also it emits 
a gamma line of 1.53 Mev, very similar in energy to that from K*°. 

Solutions of K*#Cl containing standard activities were obtained from Abbott 
Pharmaceuticals of Oak Ridge. These were administered orally to five subjects 
who had weights varying from 134 to 200 lb, and the apparent activity of 
these subjects was followed for 2 days. At the end of this time it is to be 
expected (Rapoport 1959) that, to a sufficiently high degree of accuracy, the 
K*® will be distributed throughout the body. This expectation was confirmed 
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by the fact that the K*® count, after correction for radioactive decay and 
biological elimination, remained approximately constant after 24 hours. 

Two methods of calibrating the detectors were used. The first method used 
a known mass of ordinary KCI as the standard of radioactivity. Measure- 
ments were made with equal activities of KCl in the five subjects and in 
water solution in a vessel at a standard distance from the crystal. Then 
measurements were made with ordinary KCI solution in the same jar at the 
standard distance. Intercomparison of these three measurements with the 
activity due to normal K in the body gives the mass of ordinary K in the 
body. Corrections had to be made for biological excretion of K#, measured 
to be 1% every 10 hours, and for the different self-absorptions of the dilute 
K#Cl solution and for the concentrated ordinary KCl solution. This latter 
correction was found experimentally by the use of solutions of different 
strength. Measurements were taken of the count obtained from solutions 
containing 750 g, 500 g, and 250 g of KCI dissolved in water to give in each 
case 4 liters of solution. From the results an extrapolation was made to 
infinite dilution. It was found, for example, that with the 750-g sample, the 
self-absorption correction was 5%. Confidence in the experimental extra- 
polation was reinforced by the good agreement with theoretical calculations 
of the expected correction. 

The second method made use of the Abbott Pharmaceutical standardization 
of the K*2. Knowing the amount of K*® in the body at any time (after nuclear 
decay and biological excretion correction), and the corresponding K® photo- 
peak count, the value of A in the relationship count = A X activity can be 
found. A is mainly a geometric term with corrections which are functions of 
photon energy, the variations of which with energy may be found by sub- 
sidiary experiments with different sources. This method, however, depends 
on a knowledge of the decay schemes of both K* and K*, and in both cases 
the percentage of disintegrations which lead to gamma emission is not abso- 
lutely known. Strominger, Hollander, and Seaborg (1958) quote a figure of 
18% for the fraction of K® disintegrations that give 1.53-Mev gamma rays. 
Using values quoted by these same reviewers, the post-1954 figures for the 
specific beta activity of normal K give an average of 28.6 6/sec g, and the 
average of post-1947 results on the y/@ ratio is 0.12. In our calculations, we 
have accepted an average figure of 3.43 y/sec g K. In addition to decay scheme 
uncertainties, the final result using this method involves the calibration errors 
of the K* (44%). 

In fact the calculated results based on the two methods only differed by 
4%, which is within the accuracy of the experiments, and as a result in all 
subsequent calculations the average of the two methods was used. 


RESULTS 
The spectra of 30 male first-year premedical students (average weight 
158 lb) at the University of Toronto were measured, and a typical graph 
of the spectrum of one of them is shown in Fig. 1. The spectrum shows peaks 
at 1.46 Mev and 0.66 Mev, corresponding to K*° and Cs!" respectively. This 
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spectrum is the result of 30 minutes running time on the student, with the 
subtraction of 30 minutes background. Typically such a run gave a count 
in the K*° photopeak of 1200 counts, after the subtraction of a background 
of 780 counts. The statistical error on any individual count is 4%. 
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Fic. 1. The spectrum of a subject with no radioactive history; 30-minute run minus back- 
ground. This subject is calculated to have a body content of 141 g of K and a Cs/K ratio of 


50 puc/g kK. 


In none of the 30 subjects was there evidence for any appreciable activity 
other than that due to K*° and Cs!*?, 

Using the calibration obtained as outlined above, the K contents of the 
students were obtained. The average figure obtained is that (0.212+0.013)% 
of the body weight is potassium, the error being based on uncertainties in 
calibration, geometric errors, etc. As expected, the individual potassium body 
contents show a biological variation, the graph of (number of subjects with 
a given K percentage) vs. (K percentage) being a Gaussian whose mean is 
0.212 and whose standard deviation is 0.023, both figures being percentages 
by weight of potassium in the body. 

Consideration of the spectrum of K* in the body enables the Compton 
spectrum due to K*° in a normal body to be calculated. Subtraction of this 
K*? Compton spectrum from the normal body spectrum leaves at 0.66 Mev 
a photopeak due to Cs'*’, and from this the Cs photopeak count may be 
obtained. Using the calibration, and a knowledge from theory and subsidiary 
experiments of the variation with photon energy of (i) the self-absorption 
of the body and of material between it and the crystal (ii) the capture cross 
section of Nal, and (iii) the probability of a count appearing in the photopeak 
of the energy spectrum, the activity corresponding to this Cs photopeak 
count may be obtained. 

The results show a much greater spread than do the K content measure- 
ments. The average Cs!*7 activity in the 30 students tested was (December, 
1958) 7.8 muc, or, comparing the potassium and cesium contents, 51 pyc/g 
K. The spread of individual activities is shown in Table II, which also gives 
the average liquid milk drunk by the persons in each group. 
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TABLE II 


Average daily milk 


Cs/K, No. intake of group, 
ppc/g K in group imperial pints 
<40 7 Ls 
40-59 14 2.1 
60-79 8 2.1 


>8s0 1 3.5 


DISCUSSION 

The result for the total potassium content of the body, 0.212% of the total 
body weight, is in very good agreement with other determinations of this 
quantity made on subjects from other countries by different research teams 
(Burch and Spiers 1953; Miller and Marinelli 1956; Shohl 1939). This agree- 
ment gives confidence in the results for the cesium content of the 30 sub- 
jects. As stated in the introduction, it might be expected that there is a 
latitude variation in the Cs/K ratio. Toronto (44° N.) results should not be 
expected to vary greatly from results obtained at Chicago (42°N.). In fact 
the average of 13 persons there in June, 1957, was 35 wuc/g K (Miller and 
Marinelli 1957). It will be seen that the Toronto figures are 40° % higher than 
the 1957 Chicago results, but this is perhaps not unexpected owing to the 
rather large number of nuclear bomb tests during the summer of 1958.* 

From the table, it will be seen that there is here evidence for a correlation 
of the Cs'8? content of the body and milk intake, as has been previously 
suggested (Miller and Marinelli 1957). It is of interest to note that (i) the 
person with the lowest Cs/K ratio drank no milk, and (ii) the Cs/K ratio 
for a sample of powdered milk produced in Ontario was measured to be 
31 wuc/g K (McNeill and Trojan 1959). 

It is proposed to follow up this work by testing the same group of students 
at intervals for the next few years so that any change in the activity of this 
control group may be noted. 
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residents is 52 uuc/g K, in good agreement with our figure. 
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POLARIZATION OF RADAR ECHOES FROM AURORA! 
A. KAvVADAS AND D. G. GLAss 


ABSTRACT 


Auroral radar echoes at very high frequency received at antennas sensitive 
to linearly polarized components in directions symmetrical to the plane of 
polarization of the vertically polarized transmitted wave and the ground plane 
indicate that the received wave, in addition to an unpolarized component, 
contains a linear component of polarization tilted in the general direction of 
the earth’s magnetic field lines. 


INTRODUCTION 


Radar echoes from aurora when observed at frequencies in the lower part 
of the very high frequency band are often found to exhibit interesting polari- 
zation effects. In some cases the received signal is found to be equally strong 
on two orthogonally polarized antennas, one of which has the polarization of 
the transmitting antenna. Such echoes are said to be depolarized and the 
degree of depolarization for various circumstances has been investigated by 
several workers (Harang and Landmark 1953, 1954; McNamara and Currie 
1954). The most complete study seems to be that of McNamara and Currie 
(1954), who found that at 56 Mc/sec an appreciable number of echoes showed 
some or complete depolarization. They also found that when the transmitted 
wave was horizontally polarized the ratio of the horizontal component to the 
vertical component of the received wave was a function of range. Echoes less 
depolarized were more frequently observed at ranges of 500 to 900 km, a 
region of the sky where the direction of the incident wave and the direction 
of the earth’s magnetic field were almost perpendicular. Most of the workers 
found that the degree of depolarization decreased with increasing frequency. 

The force of the earth’s magnetic field may influence the polarization of 
the aurorally reflected signals either by means of Faraday rotation or indirectly 
by providing an axis of symmetry of the reflecting clouds of ionization. The 
clouds can be visualized as elongated along the lines of force of the magnetic 
field, and under certain circumstances this may lead to a backscatter cross 
section for the clouds, either singly or collectively, which may vary with the 
plane of polarization. 

This paper concerns an investigation of the polarization of aurorally reflected 
signals which had as its main objective the discovery of any effects which are 
related to the earth’s magnetic field. 


EXPERIMENTAL 


If the change in the state of polarization of the incident wave is associated 
with the angle between the plane of polarization of the incident wave and 

‘Manuscript received January 19, 1959. 

Contribution from the Institute of Upper Atmospheric Physics, University of Saskatchewan, 
Saskatoon, Saskatchewan. Work was performed under Defence Research Board of Canada, 
Contract Nos. PG 69-613013 and PG 69-813002. 
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the direction of the earth’s magnetic field, the received wave may behave in 
a way which depends on whether the echo originates to the east or to the 
west of the observer’s magnetic meridian. Consider an observer located near 
the receiving antennas; to such an observer the lines of force of the earth’s 
magnetic field will appear to be vertical when viewed along the magnetic 
meridian, but they will diverge almost symmetrically, making an angle with 
the vertical along eastward or westward directions. The reflected waves from 
these directions might then contain linearly polarized components tilted away 
from the vertical. 

In order to detect such components, receiving antennas sensitive to linearly 
polarized waves in directions +45° and —45° from the vertical were selected. 
These were a pair of two identical dipoles crossed at their centers and mounted 
in the form of an X. 

At the low frequencies employed ground reflections are almost inevitable. It 
was, therefore, important that the two dipoles were equally affected by these 
reflections. The X antenna satisfied this condition by having its two elements 
symmetrically located with respect to the ground plane. This property of 
symmetry is demonstrated in Fig. 1(@). The electric vectors associated with 
the wave received at y, x are proportional to the sum and difference of the 
same two quantities, the vertical and horizontal component of the received 
wave, v and h respectively. 


~ 
---~> 
—— 


GROUND 





(a) (b) 


Fic. 1. (a) Symmetry properties of X antenna. (6) Experimental arrangement. 


No transmitting equipment was installed. The presence of an I.G.Y. auroral 
radar (McNamara 1958), located in the vicinity, provided a vertically polarized 
transmitted wave at the frequency of 48.2 Mc/sec. 

The complete experimental arrangement is shown in Fig. 1(0). Signals 
from the receiving dipoles y and x were alternately displayed on the two 
halves of the upper trace of a dual beam oscilloscope. On the lower trace of 
the same oscilloscope the signals from the two antennas, marked E and W 
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in Fig. 1(0), were similarly displayed. These two antennas were vertical 
dipoles with corner reflectors, directed in an eastward and westward direction 
respectively, and indicated whether the signal was received from the western 
or the eastern part of the sky. In order to distinguish the traces corresponding 
to the two corner reflector antennas E and W, an identifying pulse appeared 
on the western corner reflector trace. Both traces were locked together by 
means of a synchronizing circuit contained in the switching unit, so that the 
order of appearance of the signals could not change by mistriggerings or 
other interference. Range markers were displayed at 100-km intervals and 
the complete display was sampled by a 16-mm movie camera at a rate of 
2.75 frames per minute. The exposure time was 1.8 seconds per frame. Samples 
of these films may be seen in Fig. 2. 

Because cross talk between the two dipoles would tend to reduce the 
difference in intensity between the signals received at the two antennas, 
particular attention was paid to the signal rejection between the antennas. 
The cross talk was kept at all times lower than —22 db. 

The signals from the dipoles y and x were fed to the inputs 1 and 2 by 
low-loss coaxial cables (Fig. 1(b)) and the impedances were matched at the 
inputs of the two radio-frequency grounded grid amplifiers 1 and 2, instead 
of at the antennas, in order to minimize instabilities that could result from 
varying outdoor weather conditions. 

The switching amplifiers 1 and 2 were constructed with identical charac- 
teristics, their ‘‘on’’ times being those for which the grids were at ground 
potential. Both amplifiers 1 and 2 have a common output connected to 
“receiver 2”’ (Fig. 1(b)). The measured cross talk between these amplifiers 
was —48 db. 

Variations in the characteristics of amplifiers 1 and 2, due to aging of tubes 
or other factors, were checked continuously by means of a short vertical 
antenna that was placed in front of the X antenna at a distance sufficient 
to prevent interference with the receiving antenna field. A delayed pulse was 
impressed on this antenna; this pulse could be used for quick calibrations of 
the receiving equipment and a continuous check on the performance of the 
system. Due to its location with respect to the E, W antennas, this antenna 
could only be seen by the western corner reflector and thus it also supplied 
the identifying pulse. 


OPERATIONAL DETAILS AND OBSERVATIONS 

The receiving equipment was operated as continuously as possible for a 
period of approximately 2.5 months (March 4 to May 18, 1958) during the 
night hours. The number of daytime echoes was too small to warrant an 
attempt to record them. 

The plane of the receiving antennas could be rotated about a vertical axis 
enabling recordings to be made with the antenna facing east, northeast, north, 
and northwest. The results were tabulated with reference to antenna position. 
Of the 513 usable samples obtained, only 24 were received from the north- 
western direction. Severe local interference was mainly responsible for the 
lack of usable echoes from that direction. 





PLATE | 





Fic. 2. Some typical radar echo samples: (a) echo from the west, y < x; (b) echo from 
the east, y > x. 

(a, i), (b, i): signals from the x dipole appear on the left-half trace signals, from the y 
dipole on the right-half. 

(a, ii), (6, ii): signals from the E antenna appear on the left-half trace signals, from the W 
antenna on the right-half. 
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The received signal intensities at the antennas y and x, denoted by y and x 
for the remainder of this paper, were divided into the following three classes 
for each position of the antenna: 

(i) samples for which y > x indicating a linear component of polarization 
tilted in the general direction of the y dipole; 

(ii) samples for which y = x indicating a vertically or horizontally polar- 
ized wave but also a circularly or randomly polarized wave; 

(iii) samples for which y < x indicating again a linear component of polari- 
zation tilted in the general direction of the x dipole. 

Figure 3 shows a histogram of frequency of occurrence of samples with y > x 
and y < x against antenna position. For y > x, the number of samples is 
largest when the antenna faces in an easterly direction. The number decreases 
as the antenna is rotated towards the west. The converse is true for echoes 
with y < x. Because the receiving antenna could not be rotated until it faced 
west, the dotted class interval at the end of the histograms of Fig. 3 represents 
the number of samples with y > x and y < x for echoes received from directions 
west of the magnetic meridian but with the antenna facing north. 

Figure 4 shows a similar plot but for echoes with y = x. The last interval 
in this histogram was again obtained from the north position of the antenna. 
This curve shows a maximum for echoes from the north. 

Figure 5 shows a plot of all echoes with y > x and y < x, obtained when 
the antenna was directed towards north. The direction of arrival of the echoes 





E, E>w,E-WEw, W E, E>WE*WE<W, W 


Fic. 3. Percentage echo occurrence plotted against position of antenna: (a) echoes with 
y > x; (b) echoes with y < x. ; . i 
Fic. 4. Percentage of echoes of equal intensity on both dipoles (y = x) plotted against 


antenna position. 
Fic. 5. Percentage echo occurrence plotted against azimuthal direction of echo arrival. 
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is indicated by the presence of the signal in either the E corner reflector for 
echoes arriving from the east, or the W corner reflector for echoes arriving 
from the west. Five points were obtained for this curve: 

(i) echo only on the trace corresponding to the eastern corner reflector, 
denoted by E; 

(ii) echo on both the eastern and western corner reflectors, denoted by 
E=W,E>W,E < W; 

(iii) echo only on the western corner reflector, denoted by W. This histo- 
gram appears to be quite similar to that obtained in Fig. 4. 

The histograms in Fig. 3 to Fig. 5 indicate that the receiving equipment 
detects a linear component of polarization tilted away from the vertical and 
in the same direction as the general direction of the lines of the earth’s magnetic 
field. The information obtained by this method of observation cannot provide 
the exact direction of the plane of polarization of the reflected wave. Therefore, 
it is not possible to say whether the direction of polarization of the received 
wave is exactly that of the magnetic lines of force. 
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Fic. 6. Ratios of intensities y/x plotted against frequency of occurrence for three position 
of the antenna: (a) antenna facing toward the east; (b) antenna facing toward the northeast; 
(c) antenna facing toward the north. 
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It is possible to find an additional relation between the direction of the 
magnetic lines and the direction of the plane of polarization of the received 
wave. The magnetic lines diverge more rapidly for directions well removed 
from the magnetic meridian than they do for directions near the magnetic 
meridian. If the ratio y/x is plotted against frequency of occurrence, the 
resulting curves should show greater dispersion for echoes received when the 
antenna faces toward the east rather than when it faces north. 

This effect may be seen in the histograms of Fig. 6 (a), (6), and (c). Since 
the number of samples was different for each case, the histograms have been 
normalized by expressing occurrences as percentage of the total number of 
samples received from the east; it is interesting to note that the histogram 
of Fig. 6(a@) peaks in the interval 1.1—-1.3. The displacement of the maximum 
with increasing azimuthal angle is also suggestive of the greater tilt in the 
average linear component of polarization for echoes arriving from that 
direction. 

As mentioned before the number of echoes obtained from the west does not 
allow their statistical treatment and, therefore, no such curves have been 
plotted. 

It has been reported by Harang and Landmark (1954) that sharp echoes 
tend to retain the plane of polarization of the transmitted wave while diffuse 
echoes become depolarized. Most of the received echoes during the present 
investigation could not be separated into diffuse and discrete echoes with 
certainty and no attempt was made to treat them separately. 


CONCLUSIONS AND DISCUSSION 

The results obtained indicate that for a large number of echoes the plane 
of polarization of the transmitted wave is not maintained upon reflection. The 
work of Harang and Landmark (19538, 1954) and McNamara and Currie 
(1954) has already shown that the reflected wave tends to be randomly 
polarized. In addition, the reflected wave is found to contain a linear com- 
ponent of polarization other than the vertical. This component appears to be 
tilted in the same general direction as the earth’s magnetic lines of force. 

This finding does not exclude the presence of other components. It should 
be noted that even if the reflected wave is linearly polarized, the received wave 
will, in general, be elliptically polarized as a result of ground reflections. There- 
fore, the presence of circular components in addition to the linear components 
is almost certain. This is not true if the reflected wave retains its plane of 
polarization (vertical). 

In all positions of the receiving antenna signals y = x have been obtained. 
These may be due to (i) a completely unpolarized signal, (ii) a vertically or 
horizontally polarized signal, and (iii) a circularly polarized signal. The 
experimental arrangement used does not distinguish between these cases. It 
is interesting to note, however, that in all antenna positions in Fig. 6 (@), 
(6), and (c) the frequency of occurrence of ratios y/x lying between 0.9 and 
1.1 is quite large. This suggests that not all echoes contain a tilted component 
of linear polarization, and that a large number of them fall into one or other 
of the groups (i), (ii), and (iii) mentioned above. 
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The general change in the state of polarization of the incident wave may 
be attributed to any one of several factors and possibly to a combination 
of these factors. 

A plane polarized wave, scattered by an assembly of isotropic or non-iso- 
tropic scatterers, may, under certain conditions, change its plane of polari- 
zation. In addition, unpolarized components should be expected in the scattered 
wave, though the depolarizing mechanisms may be different for the two types 
of scatterers. 

Multiple scattering within the assembly, though a second-order effect, will 
tend to depolarize the incident wave whether the scatterers are isotropic or 
anisotropic and the degree of depolarization will depend on the spacing 
between scatterers. Alternatively, if the scatterers have a sufficient electron 
density, they may be asymmetrical to the frequency used and if the degree 
of anisotropy varies randomly from scatterer to scatterer, an appreciable 
part of the incident wave will become depolarized. The percentage of the 
received signal that becomes unpolarized depends on the randomness of the 
degree of anisotropy within the assembly. In addition to the depolarizing 
effects which can be associated with the target, scatter from ground irregu- 
larities or objects near the antenna that may contribute to the depolarization 
of the wave should not be overlooked, although it cannot account for complete 
depolarization. 

The unpolarized components along with the circular components that result 
from direct ground reflections will tend to reduce the difference in intensity 
between the signals received at the two antennas y and x. Thus the actual 
tilt in the plane of polarization of the received wave will be masked to a 
certain extent. The clustering of the ratios y/x around the value 1 in Fig. 6 
may be attributed to this effect. 

Faraday rotation, a basically different mechanism, may also contribute to 
the change in the state of polarization of the incident wave, whenever the 
ionic concentration of the echoing region attains sufficiently high densities. 
This effect, however, is not a likely cause for the rotation of the plane of 
polarization of the reflected wave. The echoes originate in regions where the 
magnetic lines are almost perpendicular to the propagation vector of the wave 
and so the Faraday rotation should be correspondingly slight. Even if such 
an effect were present, the received wave would exhibit linear components of 
polarization in random orientations depending on the depth and degree of 
ionization penetrated by the incident wave. Therefore, one may conclude that 
any residual Faraday rotation that exists cannot account for the observed 
tilt in the plane of polarization of the incident wave. 

As mentioned previously the incident wave may be cross-polarized by 
assemblies of scatterers, some of which are sufficiently elongated along the 
earth’s magnetic field lines to be non-isotropic to the incident wave. Such 
scatterers would have different reflectivities for polarizations parallel to and 
perpendicular to their axes of symmetry. The non-isotropic reflectivity of 
these scatterers may change during the lifetime of the scatterer to the extent 
that it may even become isotropic. Since the scatterers may be formed at 
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different times and places within the assembly, the values of their anisotropic 
reflectivities would be expected to fluctuate about a mean, making possible 
all the observed values of y/x in Fig. 6. 

An attractive alternative to this consideration (Forsyth 1958) may be that 
the tilt in the plane of polarization of the incident wave could result from 
multiple scattering, by scatterers which are elongated along the lines of force 
of the earth’s magnetic field, but are sufficiently large to exhibit individual 
backscattering cross sections which are independent of the radio wave polari- 
zation. Such scatterers will tend to radiate secondary waves with electric 
vectors on a plane parallel to the magnetic lines. Successive scatterings of such 
waves will result in the enhancement of a component with the electric vector 
parallel to the earth’s magnetic lines. In addition, this mechanism will provide 
an abundance of randomly polarized components consistent with the multitude 
of ratios y/x = 1 for all positions of the antenna. 
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THE CALCULATION OF THE ELECTRIC POTENTIAL AND THE 
CAPACITY OF A TORE BY MEANS OF TOROIDAL FUNCTIONS! 


S.C. Lon? 


ABSTRACT 


Expressions in terms of toroidal functions for the electric potential and the 
capacity of a tore with respect to infinity are derived. Some numerical calcula- 
tions of the potential are presented for tores of different sizes in order to illustrate 
the derived results. 

The electric capacity of a tore is also calculated and the numerical results are 
compared with those of a sphere of radius a. 


INTRODUCTION 


The electric potential and the capacity of a tore with respect to infinity 
are usually calculated approximately, assuming that the diameter of the 
conductor is small compared with that of the tore. The exact solution of the 
problem has been treated by Hicks (1881) and Dyson (1893). The latter ex- 
pressed the potential of an anchor ring of circular cross section in terms of an 
infinite series. 

In the present paper, expressions for the potential and the capacity of a 
tore are derived by means of toroidal functions. Some calculations of the 
potential and the capacity are presented for tores of different size in order to 
illustrate the derived results. 

ELECTRIC POTENTIAL OF A TORE 

The potential ¢@ of an electrostatic field always satisfies the differential 
equation known as Laplace’s equation. The general solution of Laplace's 
equation in toroidal co-ordinates is given by Carter and Loh (1958) as follows: 


(1) & = V(cosh u—cosv) D> {dnba(u)+bngn(u)} cos (nv-+an) 
n=0 
where a,, d,, and b, are arbitrary constants, 
p,(u) and q,(u) are known as “toroidal functions’’ of the first and 
second kinds, respectively. 
The toroidal co-ordinates, u, v, w, are related to the well-known Cartesian 
co-ordinates, x, y, z, by the following expressions: 


( a sinh u cos w 


x= 
cosh u—cos v 


] a sinh u sin w 


(2) on eee me 
y cosh u—cos v 


asinv 


cosh u—cos v 


ne 
| 
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When the central plane is a plane of symmetry, a, = 0; therefore using 
eq. (1), the potential (with respect to infinity) outside a charged tore is given 


by: 
= +/ (cosh u—cos 2) > {dnPn(U)+bngn(u) } cos nv. 
n=0 


The constants a, and b, must be so chosen as to make ¢@ = V on the tore 
(where u« = 1, shown in Fig. 1), and @ = 0 at infinity; thus 


(3) b= 0; 
and 
(4) V = (cosh u;—cOos v) AnPn(Ui) COS nv. 





Fic. 1. The tore. 


In order to determine the arbitrary constant a,, we multiply both sides of 
eq. (4) by cos nv dv and integrate from 0 to 2x. Equation (4) yields to 


7 


(5) 5 AuPa (tr) = —— 


Jo V (cosh u;—cOos v) 


V cos nv 





—- dv 


where 6, = 1 when m > 1, but 69 = 3. 
The integral may be evaluated readily by using the following identity: 


(6) { cos nv dv 





—_—— ———- = oC 1). 
o V (cosh u—cos v) 2 aul 


Substituting eq. (6) into eq. (5), we obtain 


(7) _ 2V2V, Qn (tr) 
Tv 


on a(t) ” 


Hence the general expression for the electric potential at any point (u, v) 
outside the tore (uv = u,) having a constant voltage V is 


ts) Pn(u) cos nv. 
¢) 


(8) ¢ = V V (cosh u—cos v) v2 ae "Pal 
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The potential outside a tore given by eq. (8) is calculated and plotted in 
Fig. 2. 
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Fic. 2. Potential outside the tore on the plane v = x. 

The above expression for the potential of a tore at any external point is 
identical with that given by Hicks (1881). This problem was also investigated 
by Dyson (1893) with a different approach, and his results were expressed as 
an infinite series. 

THE ELECTRIC CAPACITY OF A TORE 

In order to determine the capacity of a tore, it is necessary to take the 

surface integral of (1/v)e(0@/0n) over it. 


Let C be the capacity of the tore with respect to infinity, 


; r _, on 0 9 
(9) c=t{ 2ap 1.52. gy 


v dv Ou On 
1/2 “al > 5 1 ) sinh u a 
= 4\/2easinhu Sy a re re 
“—, "cosh u—cos v 12./(cosh u—cosv)*" 


aus _ dba | gu(tts) ae 
++/ (cosh u—cos v) - teins nv dv. 


Putting 
dpn _ _2n+1 


du 2 sinh u (Pat1—cosh u pn), 


(10) 
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and simplifying the expression, eq. (9) becomes 


2n+1 


, , n { 
2 [C’Pndn — PnQn+i t Pn+19n— C’Pndn] fut 


Pur! 








(1) C=ta> in) 


where C’ denotes cosh uw. 
To evaluate eq. (11), we shall employ the following identity: 





9 
(12) Pn+19n — PnQn+1 = 2n+1 . 
Then substituting eq. (12) into eq. (11) yields 
oo i ) 
(13) C= 16ae > 3, 2), 
n=0 n (M1) 


It may be interesting to compare the capacity of a tore with that of a 
sphere having radius a. The capacity of a sphere carrying a total charge Q 
uniformly distributed is well known and can be expressed as 
(14) C, = 4rea. 


Hence the ratio is 


(15) 4S 5 galur) 


T n=0 ” Pn (t1) : 


The ratio given by eq. (15) is plotted in Fig. 3. Instead of using wu as an in- 
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Fic. 3. Capacity of the tore/capacity of the sphere. 
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dependent variable, the dimension ratio r/R (see Fig. 1) has been employed. 
The relation between u and r/R is 


(16) u = arc sech (r/R). 


The capacity of a ‘‘ring without hole’’ with respect to infinity has recently 
been calculated by Bouwkamp (1955), using the method of images. He ex- 
pressed the capacity in terms of the ratio (volume/surface area of the ring). 


CONCLUSIONS 


Expressions for the potential and the capacity of a tore with respect to 
infinity have been derived. The potential distribution in the plane v = 7 has 
been shown (in Fig. 2) to be mainly determined by the size of the tore. For 
tores 0.9 < u < 1.1; that is, the diameter of the conductor is large compared 
with that of the tore, the potential in the plane v = = is nearly constant. 

The capacity of a tore is directly proportional to that of a sphere, and this 
may be seen in Fig. 3. 
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HIGH RESOLUTION RAMAN SPECTROSCOPY OF GASES 


XII. ROTATIONAL SPECTRA OF C,H, AND C.D, AND THE STRUCTURE OF THE 
ETHYLENE MOLECULE! 


J. M. Dow.inc? anv B. P. STOICHEFF 


ABSTRACT 


The pure rotational Raman spectra of C2H, and C2D, were photographed in 
the second order of a 21-ft grating. The resolution achieved was such that 
several lines due to single transitions were observed (and identified) in the spectra 
of both molecules. An analysis of these lines based on the non-rigid asymmetric 
top yielded the following rotational constants for the ground states: 


C2H4: Ao = 4.828+0.009, Bo = 1.0012+0.0009, Cy = 0.8282+0.0004 cm™; 
C2D,4: Ao = 2.432+0.008, Bo = 0.7369+0.0012, Cy = 0.5630+0.0006 cm—. 


The structural parameters obtained from these constants are ro(C==C) = 1.339 
+0.002 A, ro(C—H) = 1.086-0.003 A, and 2 H—C—H = 117°34’+30'. These 
values are in agreement with the recent preliminary results of the electron 
diffraction data. There is also agreement between the above values and those 
obtained from a recent study of the rotation—vibration spectra of C2H, and 
C.D,, although the agreement between the ground state rotational constants is 
not as satisfactory. 


A. INTRODUCTION 


The unique importance of the molecular structure of ethylene for valence 
theory has made this molecule a favorite for spectroscopic studies at high 
resolution. One of the first was the study of the pure rotational Raman 
spectrum by Lewis and Houston (1933). They observed a number of S branch 
lines (AJ = 2, AK = 0) which were interpreted in terms of a near-symmetric 
top molecule. In 1939, a study of the photographic infrared spectrum was made 
by Thompson (1939). Later, Gallaway and Barker (1942) investigated the 
infrared spectra of both C2H, and C.D,. They noted an irregularity in intensity 
of the rotational lines, but the analysis was again based on the symmetric 
top. It is this study which has led to the much quoted parameters r(C=C) 
= 1.353 A, r(C—H) = 1.071 A, and ZHCH = 120°, for the structure of the 
ethylene molecule. 

The pure rotational Raman spectrum was reinvestigated by Romanko, 
Feldman, Stansbury, and McKellar (1954) at much higher dispersion than 
that used by Lewis and Houston. In addition to the ‘‘symmetric’’ top R and 
S branches they also observed another series of lines (called the S’ branch) 
which was due to “asymmetric” top transitions. This was the first report 
of the departure of the spectrum from that of a symmetric top molecule and 
encouraged the study of the ethylene spectrum at even higher resolution. 

While the present investigation was in progress, Allen and Plyler (1958) 
reported their results of a study of the infrared spectra of C2H,s and C:Ds. 
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They observed an alternation in intensity of the A, = 0 and K, = 1 lines of 
the parallel bands studied, thus confirming the planar symmetrical structure 
of ethylene (point group V,). Their analysis of these two parallel bands (up 
to J < 10) yielded the following structure: r(C==C) = 1.337 A, r(C—-H) 
= 1.086 A, and ZH—C—H = 117°22’. These values differ significantly from 
the earlier values quoted above, but are essentially the same as those reported 
by Bartell and Bonham (1957) from the preliminary results of electron 
diffraction data, namely r(C==C) = 1.334 A, r(C—H) = 1.085 A, and 
ZH—C—H = 116°. 

The present paper is a report of the analysis of the rotational Raman 
spectra of C,H, and C.D, photographed with a 21-ft grating. The resolution 
achieved was high enough to warrant analysis based on the non-rigid asym- 
metric top. Since this is the first time that such an analysis has been carried 
out for a Raman spectrum, a fairly detailed account will be given, including 
a review of the theory of the asymmetric top energy levels. 

In earlier papers of this series many examples of the pure rotational spectra 
of linear and symmetric top molecules were shown. All of the spectra were 
simple in appearance; the lines were evenly spaced and had a regular intensity 
distribution. The rotational numbering was unambiguous and the analysis 
straightforward. In general, this simplicity is also found in many of the 
rotation—vibration bands of linear and symmetric top molecules. 

On the other hand, the rotational spectra of asymmetric top molecules 
exhibit such a complex appearance (and the spectra of C2H, and C.D, are 
no exceptions) that the identification of rotational lines is not evident, and 
in general one must resort to other means of analysis. In principle the usual 
procedure is very simple, but in practice it may be tedious and time-con- 
suming. The starting point is a calculation of the spectrum with assumed 
rotational constants. The calculated spectrum is then compared with the 
observed spectrum and the deviations noted. The rotational constants are 
modified until the agreement between the calculated and observed spectra is 
considered satisfactory. 

Under favorable circumstances, particularly when the spectrum is observed 
at high resolution, the above procedure may lead to the identification of 
individual rotational lines, that is, lines due to single transitions. These can 
then be used to determine the rotational constants explicitly. Such a procedure 
is customary in the analysis of the microwave spectra of asymmetric tops. 
Fortunately, several lines in the Raman spectra of C2H, and C.D, reported 
here were identified as single transitions and, as shown below, their analysis 
has led to the evaluation of the rotational constants of these molecules. 


B. EXPERIMENTAL DATA 
The apparatus, including the mirror-type Raman tube, mercury lamps, 
and 21-ft grating spectrograph, has been described in the first paper of this 
series (Stoicheff 1954). As in the earlier work, the present Raman spectra 
were excited by the Hg 4358 line and photographed with a reciprocal linear 
dispersion of 6.75 cm~'/mm on Kodak 103a0 photographic plates. Slit widths 
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of 0.02 and 0.03 mm were used. Second- and third-order iron lines excited in 
a hollow cathode provided standards of wavelength. 

In preliminary experiments, it was found that the rotational lines of CH, 
were somewhat broadened by gas pressures of 1 atmosphere. Sharper lines 
were obtained at a pressure of } atmosphere resulting in considerably better 
resolution of the spectrum. The final photographs of the C2H, and C.D, 
spectra were therefore taken at this pressure. Five photographs of each 
spectrum were obtained with exposure times ranging from 1 to 24 hours. 
Several photographs of higher density were obtained with a cylindrical lens 
placed in front of the photographic plate. 

Examples of the pure rotational spectra of C2H4 and C.D,, showing the 
Stokes branches, are reproduced in Fig. 1. Also, a microphotometer trace of 
part of the C,H, spectrum is shown in Fig. 2. The most prominent feature in 
the spectra of both molecules is the irregularity in line spacing and in intensity 
distribution. Of course, this irregularity is to be expected for the rotational 
spectra of asymmetric top molecules and is therefore an obvious confirmation 
of the asymmetric top structure of ethylene. 

The Raman displacements of the C.H, rotational lines are listed in Table I 
and those of C.D, in Table II. The values are averages of the Stokes and 
anti-Stokes lines measured on five plates. Sharp lines of medium or higher 
intensity have an accuracy of +0.03 cm~! and the remainder have an accuracy 
of about +0.06cm-'. These errors are based on the internal consistency 
amongst the five plates for each molecule. 

One of the problems presented by the irregular spacing of the rotational 
lines was the difficulty in recognizing grating ghosts particularly in the region 
of the exciting line. For this purpose a thorough study of the ghost positions 
and intensities in the region of \4358 A was carried out. Several photographs 
were obtained using one of the lamps as a source and accurate measurements 
of ghost positions were made as in the case of the Raman spectra. A direct 
comparison of the Raman spectra and ghost plates led to the elimination of 
some lines and to the identification of others as blends with grating ghosts. Of 
the latter group, those which fell within +0.2 cm~' of a ghost and appeared 
to have a higher intensity than that expected for the ghost alone were retained 
up to the final stages of analysis. In Tables I and II only those blended lines 
for which an assignment could be made were included (in parentheses). 


C. REVIEW OF THEORY 
1. Rotational Energy Levels 
The rotation energy levels of a prolate symmetric top (B = C) are given by 
(1) F(J, K) = BJ(J+1)+(A —B)K? 


in the usual notation (Mulliken 1955). J can take the values 0, 1, 2,... and 
for each J there are J+1 sublevels with K values 0, 1,2,...J. Each sub- 
level with K # 0 is doubly degenerate. For only a ‘‘slight’’ deviation from the 
symmetric top, K remains a “‘good’’ quantum number and the energy levels 
can be represented by eq. (1) but with B replaced by 3(B+C). Indeed, this 
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TABLE I 
The pure rotational Raman spectrum of ethylene 

















4v obs. Transition Av calc ” 4v obs. Transition Avcale. I* 4v obs. Transition Av calc. 
(9.606)® —$3(1) 9.576 778 (20.191) Ss(4) 20.082 577 Ss(7) 31.176 689 
R:2(9) 9.567 21 Ss(4) 20.392 1346 S0(7) 31.183 295] 
R:(10) 9.615 16 Se(4) 20.170 338 Si0(7) 31.150 175 
R:(11) 9.540 12 S7(4) 20.179 338 Su(7) 31.150 175 
Ss(4) 20.153 123 RK(16) 
(10.508) R,(5) 10.752 270 S9(4) 20.154 288 SK(7) 
Ri(7) 20.011 38 
10.986 Rs(5) 11.001 467 R3(9) 20.625 177 31.922 Ss(7) 31.950 548 
Ri(5) 10.996 540 || Rs(10) 20.555 114 
Rs(5) 10.996 1259 Rs(10) 20.028 114 32.238 S3(7) 32.239 632 
RxK(5) R10) 20.223 388 
Rs(10) 20.210 166 32.917 S2(8) 32.786 585 
(11.249)  Ri(4) 11.281 120 RK(10) Rx(13) 32.523 61 
Rx(5) 11.343 631 Rx(17) 
Rs(5) 11.015 467 21.000 S3(4) 20.999 763 
33.381 Si(8) 33.396 1422 
(12.173) S2(2) 12.186 819 21.966 Re(11) 21.727 87 Ri(11) 33.533 12 
R:(11) 22.084 129 
(12.325) R,(6) 12.368 445 Real 22.055 301 33.975 i777 = ©=©33.980 
K ( 
(12.811) Si(2) 12.721 2338 34.847 S4(8) 34.501 516 
S4(2) 12.801 387 22.529 $2(5) 22.547 722 Se(8) 34.909 416 
Ss(2) 12.881 904 Rs(11) 22.611 87 S7(8) 35.049 416 
Rs(6) 12.871 347 Ss(8) 34.862 304 
Re(6) 12.836 347 (23.162) Si(5) 23.285 773 So0(8) 34.871 710 
Rx(6) 12.835 1042 Ri(8) 23.216 28 Si0(8) 34.829 198 
Rs(6) 12.834 447 || R;(10) 23 .342 58 Sn(8) 34.829 198 
Rx(6) 1] : ; RxK(17) 
|| 23.801 Sa(5) 23.707 1365 SK(8) 
(13.682) — Ss(2) 13.396 819 || Se(5) 23.851 399 
R:(6) 13.455 191 || S7(5) 23.872 399 35 . 869 S2(9) 36.168 537 
Ri(7) 13.869 138 | Ss(5) 23.827 499 S3(8) 35.919 585 
| S9(5) 23.830 214 Ss(8) 35.859 1204 
(14.034)  R(5) 14.035 78 | Si0(5) 23.830 72 Rx(14) = 35.900 20 
Su(5) 23 .830 72 
14.666 Rs(7) 14.744 260 | R;(12) 23.957 234 36.690 Si(9) 36.694 557 
Re(7) 14.665 260 | Rs(12) 23.898 100 
Ri(7) 14.675 353 Rx(12) (38.070)  S,(9) 38.065 1120 
R3(7) 14.675 824 | 
RxK(7) (24.605) S3(5) 24.774 722 (38.531) S0(9) 38.595 396 
Ss(5) 24.208 585 S1(9) 38.824 396 
15.604 $:2(3) 15.651 798 Rs(12) 24.759 67 Ss(9) 38.555 702 
Ri(7) 15.694 323 || So(9) 38.568 301 
R,(8) 15.229 238 | 25.627 S2(6) 25.974 678 Si0(9) 38.508 207 
Ri(13) 25.850 78 Su(9) 38.508 207 
(16.437) Si(3) 16.294 904 || Rs(13) 25.686 182 Sx(9) 
Si(3) 16.446 1233 || RxK(13) 
Ss(3) 16.618 528 | 39.537 S2(10) 39.539 489 
S0(3) 16.494 214 26.695 Si(6) 26.701 1672 S3(9) 39.555 537 
S$7(3) 16.498 214 | Ri(9) 26.552 21 R3(15) 39.420 36 
Ri (6) 16.948 54 | R3(11) 26.235 103 
R,(9) 16.426 77 | 39.940 Si(10) 39.977 1181 
Rs(8) 16.640 196 |} 27.476 Sa(6) 27.320 573 Ss(9) 39.772 480 
Re(8) 16.480 196 Se(6) 27.535 425 
Rx(8) 16.520 642 || S7(6) 27.578 425 || *41.615 S4(10) 41.613 441 
Rs(8) 16.518 275 | Ss(6) 27.498 268 
RK(8) 1 So(6) 27.504 626 42.207 Se(10) 42.279 371 
\| Si0(6) 27.470 134 S7(10) 42.630 371 
(17.410) = -Sa(3) 17.205 798 || Sn (6) 27.470 = 134 || S3(10) 42.254 +9289 
R,(10) 17.446 §=135 | Sx (6) 1 S9(10) 42.277 675 
| Rs(13) 27 .027 52 Si0(10) = 42.186 +=208 
(18.262) R3(8) 18.079 102 Rs(14) 27.549 61 Su(l10) 42.188 208 
R4(11) 18.274 44 RxK(14) Sx(10) 


Rs(9) 18.571 149 || 
Re(9) 18.271 149 || *28.058 Ss(6) 28.064 1336 43 .056 Si(11) 43.250 456 
R:(9) 18.366 214 S2(11) 42.898 442 


| 
Rs(9) 18.363 499 (28.486) S3(6) 28.523 678 S3(10) 43.149 489 
| 


RR(9) 
(29.284) S2(7) 29.388 632 43 .667 Ss(10) 43.679 1030 


(18.826) R4(12) 18.908 79 R3(12) 29.296 34 | 
RR(15) *45.139 Si(11) 45.141 937 
(19.139) S2(4) 19.105 763 | 











R13) 19.350 =.26' || 30.075 —Si(7) 30.069 663 || 45.885  Ss(1l) 45.948 342 
Ri(10) 29.996 16 Ss(11) 45.961 634 
(19.742) —-$,(4) 19.816 1943 || So(1l) 45.995 272 
R14) 19.607. 47: || 31.107. S47) 30.919 1278 Sio(11) 45.877 201 
R15) 19.690 16 Ss(7) 31.222 428 Sul) 45.876 201 
R16) 19.611 29 || S:(7) 31.302 428 || Sx(11) 


“Relative intensities were calculated from symmetric top matrix element with Bo = }(Bo+Co). 


‘Parentheses enclose the Av values of lines blended with grating ghosts. 
Note: Asterisks indicate the Av values of ‘‘single transitions” used for the explicit evaluation of the rotational constants. 
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TABLE I (Concluded) 
The pure rotational Raman spectrum of ethylene 
Avobs. Transition Av calc. I? || Avobs. Transition Avcalc. I Av obs. Transition Av calc. te 
| 
46.515 $i(12) 46.526 952 | 66.265 Si(18) 66.248 411 83.869  Ss(21) 83.837 94 
52(12) 6 . 24 396 S2(18) 66.198 z 
Sill) 46.682 — S7(16) 66.215 192 || (84.484) = So(21) 84.542 72 
ee: ors || (66.599) S3(17) 66.769 202 |] *85.312 Se(22) 85.335 70 
(47.662) Ss(11) 47.567 402 | S16) 66.431 508 || gg.015 —5,(24) 86.042 131 
#48.635  Si(12) 48.649 362 || 67.895 —Se(17 67.742 166 S2(24) 86.036 55 
: p a sy) : Si(23) 85.998 148 
=, 33) 49.588 351 || 68.154  Sa(17) 68.242 325 — Teen 
S2(13) 5 351 || 68.15 § 38.242: : 
Ss(12) 49.617 311 Sio(17) 68.120 111 || 86.256 $3(23) 86.226 68 
Ss(12) 49.633 251 Su(i7) 68.150 111 SK (22) 
Si(i2) 49.574 190 || 68.710 Sv(17) 68.731 139 || 86.676 $x(22) 86.007 59 
$ ® | 10(22 . ¢ 
Su(Q12) 49.571 190 || 69 966 © Si(18) 69.257 160 ao! ‘ane ; 
SK (12) ee mae *87.136  $3(22) 87.131 181 
| OF 44 ¢ 
50.256 S3(12) 50.157 396 || 69.495 Si(19) 69.545 149 Su(22) 87 .055 49 
S112) 50.364 = 311 | S19) 69.510 145 |] [87.811] 
*51.418 Ss5(12) 51.429 845 70.075 S3(18) 7 284 88.588 $6(23) 88.782 57 
Lea Saat Ss(17 70.045 1 $9(22 56 39 
#52.128 Si(13) 52.132 754 || sy} 70.197 166 —— on = 
aE 89.399 Si(25) 89.341 45 
[52.506] | 71.398 S6(18) 71.306 142 | S2(25) 89.336 44 
53.196 Si(14) 53.083 744 | SK (18) aoa Zo. Oe 
So(14 ‘ 31 Ks ‘ we .32 5 
Si3) ya Oe ee ee S:(22) 89.569 70 
¢ 220 | S10 F ) 96 é 
whe aes Su(18) 71.896 = 96 |} 90.296 ~=—-S4(28) «90.380 Ss 64 
Ss(13) Peet asnieucies : et Ss(23) 90.301 114 
S9(13) 229 || 72.673 $:(20) 72.845 290 | $0093) 90/529 40 
S10(13) 176 S20) 72.820 121 |! Se(03 i 
Su(13) 176 || Si(19) 72.626 316 || SK (23) 
SK (13) \| So(18) 72.627 279 |) (90.905)  Su(23) 90.921 40 
*54 287 S7(13) 54.290 280 || 73.562 S3(19) 73.239 145 || 92.628 Si(26) 92.638 83 
[54.775] 5Ss(18) 73.592 374 | S2(26) 92.635 35 
2 49 | i] ~ (6 6 aT) . 
Renae ’ 55.593 286 || *74.155 Sx(18) 74.158 142 | [o. uas os 
55.386 Si(14) «55.593 286 || *74.155 7 7 | Si(25) 92.630 96 
Ss(13) 55.252 323 || [74.988] —Se(19) 74.851 120 || Suh gale 
| | Olde 2.01% te 
Re 249 Si(15 56.371 278 | | 
ae Ss) 56.248 271 || 75.215 Sx (19) (93.266)  S:(23) 93.280 7 
|| (75.52 Ss(19) 75.654 = 238 |! 93. 695 S5(é 3.5 ‘ 
019 © S(14) 80.929 310 || 7-87) Se) TR gg || 98-632 Sul8H), 98.000 120 
—) Se oe Su(19) 75.658 82 . 
Sa(14) 57.134 206 | E é “|| 93.954 Ss(24) 93.915 40 
So(14) 57.248 480 || (76.105) Si(21) 76.142 103 || Si0(24) = 94..265 33 
Si0(14) 56.968 160 |} S2(21) 76.125 101 | ie 
Su(14) 56.977 160 Si(20) 75.990 113 || 94.865 Su(24) 94.803 33 
Sx(lt) || 76.561  $3(20) 76.476 122 || 95.920 Si(27) 95.936 28 
*58.234  S7(14) 58.246 249 | So(19) 76.563 102 S2(27) 95.95 28 
el Si(26) 96.088 35 
58.989 Si(15) 59.038 585 || *77.051 Ss(19) 77.071 136 | S4(26) 95.941 33 
~ aC ¢ 6 | if * (Or ia Q5 
Ss(14) 59.030 667 | 78.088 S3(19) 73.084 120 S6(25) 95.617 37 
59.641 Si(16) 59.660 562 || _. ... “— ; | 96.736 S3(25) 96.798 41 
S2(16) 59.570 235 i (78.663)  Se(20) 78.370 =:101 Pa ‘S9(24) 96.920 46 
60.555 Sx(15) 60.239 271 || (79.026) SK(20) H S9(24) 96.653 92 
So(15) 60.546 220 |) 79.363 © Si(22) 79.443 198 || 97.556 = Ss(25) 97.508 74 
Ss(15) 60.812 426 || S(22) 79.430 83 | Si0(25) 97.986 26 
Sio(15) 60.680 144 Sx(21) 79.718 = 101 || SK (25) 
Su(15) 60.689 144 S4(21) 79 .339 220 || x 
SK(15) | S3(20) 79 : 342 86 || 99.139 Si(28) 
aa S10(20) 79.312 70 || S2(28) 
61.070  So(15) 61.043 183 Sio(2 79.312 | $.(27) 99.318 
a Su(20) 79.438 = 70 || Stay) (90247 
cli Ate i Sil2 99.2 
its sie 320 || *80.474 $620) 80.486 265 || Ss(26) 99.003 
= es S9(20) 80.537 200 | Su(25) 98.724 
se « c nD One 9 | | 
oe Ses gon || 81.870 ‘|[100.042]) Ss(26) 99.987 
SA) SET |) 8876S) 81.005 «8g || 100-8ER SRS) 108. 
*63 .521 S3(16) 63.516 235 || S7(20) 81.968 101 S0(25) -66 
64.201 Ss(16)_ 64.154 192 | go. 817 S,(23) 82.743 69 [100 .962] 0) 101.076 
=), oe = S2(23) 82.734 68 |{ . y 
Su(l6) 64.518 127 Sx(2l) || 102.630 Si(29) 17 
Su(16 5 27 S3(22) 82.968 = 83 | $2(29) 17 
SK (16) Si(22) 82.668 = 78 || S3(28) 22 
series pac aku. cane Ss(21)- 83.012_—‘167 | Si(28) 20 
et nee See ee ‘Swo(21) 83.048 58 | Se(27) 102.371 23 
65.870 Si(17) 65.867 438 Su(21) 83.236 = 58 Su(26) 102.672 21 











°Square brackets enclose the Av values of weak lines whose measured positions may be in error by more than +0.06 cm=! 
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TABLE II 
The pure rotational Raman spectrum of ethylene-d,* 





























S4(9) 26.578 1681 Ri(16) 37 .338 9 


Av obs. Transition Av calc. I || Avobs. Transition Av calc. I j| 4» obs. Transition Av calc. I 
13.107 S2(4) 13.206 1581 || (20.417) Si(7) 20.407 1423 29 .935 S3(10) 29.720 1737 
Ri(6) 13.106 115 | : S2(7) 20.172 1371 R3(14) 29.828 55 
R:(9) 13.208 565 | * S3(6) 20.269 1445 Sp(10) 
Rs(9) 13.119 847 | Ss(6) 20.545 1887 
Ro(9) 13.110 730 | Ri(9) 20.502 47 30.093 Ss(10) 30.063 935 
Rio(9) 13.108 730 | R3(16) 20.428 121 Sio(10) 30.199 616 
| | Ss(10) 30.284 783 
(13.683)  Si(4) 13.678 2562 || 20.866 Rs(17) 20.820 147 Su(10) 30.210 616 
Rs(9) 13.961 366 | Rs(19) 20.976 78 Ri(13) = 30.323 18 
Re(10) 13.539 287 | Rio(15)  =21.036 209 Rs(16) 30.347 74 
R18) 13.954 48 || Rx(18) 30.161 3-118 
| 21.479 Si(7) 21.720 1838 
14.364 Si(4) 14.213 1233 || Rx(11) =. 21.811 = 164 || *30.475 S9(10) 30.478 1174 
S6(4) 14.391 760 R14) 21.077 —-280 | 
S7(4) 14.435 760 | R15) —-.21.367 =. 209 || 31.432 Si(12) 31.556 1492 
Ss(4) 14.369 298 || S2(12) 31.528 976 
So(4) 14.370 447 || ¥22.220 S6(7) 22.283 1006 || Sa(11) 31.317 1484 
R2(8) 14.503 233 || Ss(7) 22.294 1120 Ss(10) 31.580 1585 
Rs(11) 14.279 227 || S9(7) 2 747 | S:(10) =. 331.332» 935 
Rs(17) 14.460 59 || Sio(7) 2 484 
Rx(10) 14.611 «674 |) Su(7) 22. 484 31.801 S3(11) 31.924 1056 
Rs(10) 14.418 449 || Rs(13) 2 144 || 
Ro(10) 14.444 = 596 Riw(16) 22.281 170 || 32.698 Ss(11) 32.606 887 
Rw(10) 14.487 596 |} 1] Sa(11) 32.950 1136 
|| 22.597 Si(8) 22.624 2003 | Sio(1l) 32.871 613 
(14.977) — Sa(4) 15.073 1581 | S$2(8) 22.461 1295 Su(1l) 32.892 613 
Ss(4) 14.833 1849 || S3(7) 22.756 1371 Ri(14) =. 32.686 14 
Re(12) 14.813 180 || Si(7) 22.651 1006 |) Rx(15) 32.493 65 
Rs(13) 15.128 144 || Ro(16) 22.884 170 
Ré(14) 15.230 115 || *33 .319 So(11) 33.281 757 
Rs(15) 15.135 92 || (23.186) $37) 23.384 1225 Rs(17) 33.196 59 
Rs(16) 14.867 74 |] Ri(10) =. 23.006 37 |] 
R15) =. 23.054 =-150 || 33.836 S$i(13) 33.799 912 
15.632 Si(5) 15.953 1605 S2(13) 33.782 896 
S2(5) 15.548 1515 || 23.469 Rw(17) 23.452 137 S412) 33.649 920 
Ri(7) 15.518 83 || 
Rs(10) 15.786 287 || (23.897) S48) 24.164 1178 34.139 S3(12) 34.106 976 
Rs(ll) 15.684 540 |} | Ss(11) 34.158 989 
Roll) 15.787 484 || (24.708) — Si(9) 24.847 1249 || Si(11) 34.268 887 
Rwo(il) 15.769 484 S2(9) 24.738 1217 || 
| S6(8) 24.898 1001 || 35.317 S6(12) 832 
(16.349) Rr(11) 16.068 360 S:(8) 24.953 785 Ri(15) 11 
Si0(8) 24.882 559 || R3(16) 35 
16.910 16.742 1899 Su(8) 24.885 559 | 
17.024 910 Rx(12)—- 24.453 86 || 35.569 Ss(12) 722 
17.126 910 Rs(14) =-:24.897 115 S10(12) 597 
17.002 786 || Ro(17) 24.498 137 Su(12) 597 
17.006 524 Rio(18) 24.503 111 | 
16.974 194 || 36.045 Si(14) 36.044 1249 
16.975 194 || 25.240 S3(8) 25.156 1295 || S2(14) 36.034 818 
16.810 268 || S7(8) 25.500 1001 || S3(13) 36. 26 896 
16.868 288 | S9(8) 25.005 1177 || S4(13) 35.958 1272 
17.180 393 | Ri(ll) 25.483 29 | So(12) += 36.141 1082 
17.102 393 | Ri(16) 25.226 ~—-182 || Rs(18) 36.077 48 
H 
17.651 S2(6) 17.869 *26 169 Ss(8) 26.180 1767 || 36.593 Ss(12) 36.639 1379 
S3(5) 17.704 | Ro(18) =. 26.241 = 111 | 
Ss(5) 17.685 | | 37.027 Si(12) 37.173 832 
Rs(11) 17.796 || 27.034 S:(10) 27.077 1746 || 
Rx(12) 17.632 S2(10) 27.007 1137 | 37.348 So(13) 37.593 774 
| 


(18.241) Si(6) 18.188 R3(13) 27.135 =103 


Ri(8) 17.995 (37 .950) R3(17) 37 .663 41 








i 
Rs(13) 17.992 | 27.433 S3(9) 27.472 +1217 |} 
S0(9) 27.491 975 || (38.233) (15) 38.292 755 
(18.811) Rs(l4) 18.984 187 | Ss(9) 27.619 1191 || S215) 38.286 © 743 
R913) «18.515 319 S0(9) 27.722 794 || Si(14) 38.481 818 
Rw(13) 18.431 319 | Sio(9) 27.535 ~—- 600 | Si(l4) 38.250 778 
Su(9) 27.541 600 | Ss(13) 38.260 1018 
(19.100) — S4(6) 19.245 Rs(15) 27.569 92 || Sw(13) 38.243 571 
R10) =: 19.255 Ri(17) 27.601 98 |) Su(13) 38.332 «571 
R:(13) 19.283 || (38.742) 
(28.067) Ri(12) 27.925 23 || 
(19.806)  Se(6) 19.657 982 | || 39.029 S13) 39.023 848 
S7(6) 19.862 982 || (28.515)  S7(9) 28.399 975 | S9(13) 39.061 679 
Ss(6) 19.643 666 i} 
So(6) 19.654 999 || 28.878 S10) 28.962 1057 || 39.629 Ri(17) 39.638 7 
Si0(6) 19.602 364 Ss(9) 28.916 1121 | 
Su(6) 19.603 364 | || 40.005 Se(14) 40.034 714 
Rs(12) 19.990 180 || (29.228) Si(11) 29.314 1079 | S13) 40.024 774 
Rs(15) 19.808 225 | S11) 29.269 1056 | R18) 40.155 23 


R9(14) 19.917 258 
Rio(14) 19.746 9258 


“Symbols and nomenclature are given in Table I. 
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TABLE II (Concluded) 
The pure rotational Raman spectrum of ethylene-d,* 
4v obs. Transition Ay calc. I 4v obs. Transition Av calc. I 4v obs. Transition Ay calc. I 
40.582 Si(16) 40.540 1021 50.106 Ss(18) 772 60.249 
S216) 40.536 670 
S3(15) 40.686 743 *50 .566 $7(17) 50.546 536 60.697 Ss(22) 
Si(15) = 40.529 1062 || 
*50.961 Ss(18) 435 60.893 $1(25) 
40.701 Ss(14) 40.880 632 So(17) 50.984 484 S$2(25) 
41.160 Ss(14) 41.315 1166 || 51.134 [61.483] 
Sio(14) 40.937 539 
Su(14) 41.095 = 539 51.767 $i(21) 372 62.042 Su(21) 
S2(21) 367 
*42 .038 S9(14) 42.030 948 S3(20) 419 63 .093 S:(26) 
Ri(18) 41.928 5 S4(20) 403 
S6(19) 427 63.974 
*42.412 S15) 42.441 653 Si0(18) 383 
64.536 
42.807 Si(17) 42.788 = 610 52.511 Ss(19) 457 
S17) 42.786 601 $1(18) 383 || [65.104] 
S:(16) 42.901 670 
Si(16) 42.798 641 53.046 S7(18) 480 65.361 $i(27) 
S$7(14) 42.800 714 
53.645 Ss(19) 583 66.119 
42.988 
54.114 $1(22) 485 (66.726) 
43.448 Ss(15) 43.531 708 S2(22) 319 
Ss(15) 43.479 874 S3(21) 367 || [67.150] 
Si0(15) 43.633 502 | S4(21) 530 
S6(20) 378 67 .564 Si(28) 
44.155 Su(15) 43.905 = 502 
54.317 Ss(20) 604 67 .837 
45.055 S$i(18) 816 | S9(18) 653 
S2(18) 537 Si0(19) 345 68.814 
S317) 45.126 601 
Si17) 45.059 864 54.764 Sn(19) 345 69.238 
Se(16) 44.817 594 | 
So(15) 45.029 583 || 55.122 S7(19) 427 || 69.750 Si(29) 
*45.701 Ss(16) 45.697 961 55.444 || 70.309 
S7(15) 45.485 653 | 
55.733 71.154 
45.940 Ss(16) 46.046 = 533 | 
56.050 Ss(20) 345 72.118 Si(30) 
46.272 Sio(16) 46.325 463 | 
|| 56.405 Si(23) 280 || (72.750) 
46.421 S2(23) 276 
\| S3(22) 319 73.306 
46.969 Su(16) 46.764 463 Ss(21) 354 || 
| S6(21) 332 || [73.574] 
47 .278 Si(19) 482 | S9(19) 389 | 
S2(19) 476 | S4(22) 308 || 74.373 $1(31) 
S3(18) 537 | 
S4(18) 515 || [57.207]  S2(20) 378 74.946 
Ss(17) 47.163 536 | Si0(20) 307 | 
|| 75.247 
(47 .699) H SF797 
| || (76.098) 
(47.897) Ss(17) 47.843 576 58.245 Ss(21) 456 | 
Si(16) 48.069 = 594 || || 76.577 Si(32) 
Ss(16) 48.026 800 58.581 Si(24) 361 || 
S2(24) 238 || (78.837) $i (33) 
48.466 Ss(17) 48.580 725 | S3(23) 276 || 
| SK(21) || [81.128] $1(34) 
48.631 I] 
|| 58.763 Ss(22) 462 || [83.331] — $1(35) 
*49 .003 Sio(17) 48.994 424 || Su(20) 307 || 
| || [85.651] $1 (36) 
49.564 S1(20) 637 || [59.013] 
S2(20) 419 | 
S3(19) 476 || 59.336 Si0(21) 
Ss(19) 685 | i| 
S6(18) 480 || [59.629] S7(21) 332 || 
Su(M17) 49.682 424 || S9(20) 517 
“Symbols and nomenclature are given in Table I. 
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representation of the energy levels formed the basis for the analysis of the 
ethylene spectra by Thompson (1939), Houston and Lewis (1933), and 
Gallaway and Barker (1942). 

In the asymmetric top, the K degeneracy is removed: for each value of J 
there are 2J+1 sublevels having energy 


(2) Fyo(J,) = (A+ O)J(J+1)4+3(A —C)E,? (x). 


Here A, B, and C are the usual rotational constants. E,’(«) is a complicated 
function, dependent on J, on the asymmetry parameter « = (2B—A—C)/ 
(A —C), and on r(= K,—K,). K, and K, are the |K| values for the limiting 
prolate and oblate symmetric tops. For a given J, the 2/+1 values of E,7 (x) 
are ordered according to the values of the subscript 7 = K,—K, = —J, —J+1, 
—J+2...+J, with r = —J assigned to the lowest level, r = —/J+1 to 
the next lowest level, and so on. 

The values of E,7(«) cannot be represented by an explicit formula; however, 
several methods of calculation are available, as well as tables which are useful 
for a limited range of J and « values. In the present analysis, the tables of 
Townes and Schawlow (1955) for J < 12 with « given in increments of 0.01 
were used. For J > 12, the method developed by Golden (1948) was used. 
Golden’s method gives quite accurate expressions for E,7(«x) for those levels 
corresponding to the lowest r values (for prolate tops) which satisfy the 
condition (K,/J) « 1. According to the ordering of the sublevels for a given 
J, the sublevels with the smallest K, values are the lowest. It is just these 
levels that show the greatest deviation from the symmetric top levels. More- 
over, since the lowest levels will be the most highly populated, transitions 
from these levels will be expected to contribute the greatest intensity to the 
spectrum. 

Golden’s formulation may be summarized as follows: 


k—1 3-k 
( 2 )ra+n+(3)e 


b— [20+ C0’ +C2(6')?] 





E,” (x) 


The quantity @ is given by 


a 


where 6 is the characteristic value of Mathieu's equation corresponding to 46, 


gett} bere ieee + ole 


a(J) 2I(I+1) 





(3) «+1 


i 
3— UT begets 


and C,; and C; are first- and second-order perturbation corrections. Values 
of C, and C, corresponding to 8 = 0 to 10 are tabulated in Golden’s paper. 
This range of values permits the calculation of E,7(x) up to J = 27 for CoH, 
and up to J = 19 for C.Ds,. 


(J) = 


| 
bol 
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In Golden’s notation, eq. (2) can be written 
(4) Fy(J,) = CJ(J+1)+3(24 —B—C)[b+ (2—C,)0’ —C,(6’)?]. 
For high J values the term 1/[2/(J+1)] in eq. (3) becomes negligible and 


consequently 6’ approaches the value 3}|(x+1)/(3—x)!. Equation (4) then 
takes the following form: 


1 = 2 
(5) Fo(J.) = CI(J+1) +424 —-B-C)b+4(B-0) 2-0) — B= Oc, 


The expressions for Fo(J,) given in eq. (2) and (5) apply to the rigid rotor. 
The energy levels of a non-rigid asymmetric rotor correct to first order (Kivel- 
son and Wilson 1952) may be written 


(6) FU.) = Fot Ai Fo’ t+A2FoI J+) +A? (J+)? +A (I+) (P-*) 
+A (P.4)+A6Fo(P.’). 


The first term, /o(=Fo(J,)), is the energy of the rigid rotor and the remaining 
six terms arise from the effects of centrifugal distortion. The 4,’s are constants 
depending on the usual rotational constants A, B, C and on the distortion 
constants Dy, Dyx, Dx, Rs, Re, and 6,. According to Kivelson and Wilson 
(1953), these distortion constants (and therefore the six coefficients A,) can 
be calculated theoretically, if the rotational constants and harmonic potential 
constants are known. P, is the operator for the component of angular momen- 
tum along the axis of quantization for the limiting symmetric top (the brackets 
indicating average values). The quantities (P,”) can be evaluated according 
to the method of Bragg and Golden (1949). In addition the approximation 
that (P,‘) = (P,*)* is usually assumed. 


2. Symmetry of Rotational Levels 

The symmetry of a rotational level in an asymmetric top molecule is desig- 
nated ++, +-—, —+, or ——, depending on whether the rotational eigen- 
function remains unchanged or changes sign for a 180° rotation of the molecule 
about the twofold axes c and a of the momental ellipsoid. That is, the first 
sign refers to the operation C2° and the second to C,*. Dennison (1931) has 
given general rules for determining the symmetry of any level and Wilson 
(1935) has shown that the statistical weight factors for ++, +—, -+, -—— 
of C,'*H, are 7:3:3:3 and of C.!Dy are 27:18:18:18. 


3. Selection Rules for Rotational Transitions 

Since the ethylene molecule has the plane symmetric form (point group 
V,), it has no permanent dipole moment and, consequently, it cannot exhibit 
a pure rotational spectrum in the microwave or infrared regions. However, 
it does have a pure rotational Raman spectrum. 

According to Placzek and Teller (1933) the Raman selection rule for rota- 
tional transitions of an asymmetric top is 


(7) AJ=0, +1, +2. 


However, there may be further restrictions on account of the symmetry of 
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the polarizability ellipsoid. For example, in C2.H, and C2D, the polarizability 
ellipsoid has the same symmetry as the momental ellipsoid and therefore 
Raman transitions can take place only between levels of the same symmetry, 
that is 

(8) > ee, Fo 8 ey ee eer: 


Since the ethylene molecule deviates only slightly from a symmetric top, 
one might expect that the quantum number K(= K,) will be approximately 
defined and that selection rules for K, might be applicable. For a symmetric 
top (where of course K is a good quantum number) the selection rule for 
pure rotational Raman transitions is AK = 0. Classically this is due to the 
fact that the polarizability ellipsoid is cylindrically symmetric about the 
figure axis. For an asymmetric top, in general, the three principal axes of the 
polarizability ellipsoid would be unequal and Raman transitions with AK, # 0 
would also be allowed. With AJ = 0 the selection rule would be AK, = +2, 
the transition for which AK, = +1 being forbidden by the selection rule 
given in (8). Transitions with AK, = 0 would be the most probable and the 
intensity of transitions for which AK, = +2 would depend on the magnitude 
of the difference of the two polarizability components perpendicular to the 
limiting symmetric top axis. Also according to the symmetric top approxi- 
mation and the value J4/J, ~0.2 for CoH, and J4/I,z ~0.3 for C2D,, the S 
branches would be expected to have considerably higher intensity than the 
R branches. 


4. Equations for Raman Lines 

With the selection rules, eq. (7) and (8), it can be shown that the S branches 
(AJ = 2, AK, = 0) and R branches (AJ = 1, AK, = 0) fall into a number 
of series. The designation for these series is given in Table III along with 
the usual Jx,.x, notation. In the symmetric top limit all the S branches 











TABLE III 
Nomenclature for R and S branches 
Branch Transition Selection rules 
Sensi (J) own Jere “mb 2. 3.60 AJ = 2;AK, = 0 
Son(J) In gens — In, g—n41 Sm TS cied AJ = 2; AK, = 0 
Ron-i(J) Ik tin — Iu toni py ee 4J =1;AK, =0 
Ron(J) Ix fast Je. dan n = 1, 2,3, on AJ = 1;AK, =0 





would coincide to give a single series of lines with a uniform spacing of 4Bo. 
Likewise all of the R branches would also coincide to give another series 
with a uniform spacing of 2Bo. 

The expressions for the wave number displacements of the R and S branch 
lines are obtained by applying the selection rules, AJ = 1,2; AK, = 0, to 
the equations for the energy levels, eq. (2): 


Ave = R(J) = (A+C)(J4+1)4+3(A—C) [Ev 74x) — Ex? (x), 
Avs = S(J) = (A+C)(2I4+3)4+3(A —C) [Ey 7") —E,7(«)]. 


) 
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For high + (or K,) values and slight asymmetry, the above equation can 
be approximated very closely by the symmetric top expansion: 
Ave = (B+C)(J+1), 
Avs = 2(B+C)(J+3/2). 


For low and intermediate 7 (or K,) values and high J eq. (10) is not suit- 
able and one must use the asymmetric top expressions. From eq. (5) the 
equation for the S branches takes on the following form: 

(11) Avs = C(4J+6)+3(24 —B-—C)(b’—b) —}(B—C)(Ci’—Ci) 
ie. 
~ 4(22A—B-—C) 


(10) 


(C2’— C2) 


where the primed quantities refer to the upper state and the unprimed to 
the lower state of the transition. For low K, and sufficiently high J the 
quantity (6’—b) = Ab approaches a constant value. This happens to each 
of the series in pairs, first in S; and S, and then (as J becomes higher) in all 
the other series S; and S,, then S; and Sg..., and so on. The values of Ab 
for large J can be easily obtained from Tables Relating to Mathieu Functions 
(1951). From Golden’s work a similar situation seems to occur for the differ- 
ences (C;’—Ci) = AC; and (C.’—C:) = AC,. Approximate values are found 
by extrapolation. The term containing AC: turns out to be quite small 
(~10-* cm!) for C2H4 and C.D, and can be neglected. With this approxi- 
mation eq. (11) can be written 


(12) Avs = C(4J-+6)+4(24 —B—C)Ab—}(B—C)ACi. 


The general behavior of the S branches can now be inferred from equations 
(10), (11), and (12), for a slightly asymmetric molecule. First of all there 
should be a series of rather broad lines with a spacing of 2(Bo+Cy) corre- 
sponding to the S branch lines for large K, eq. (10). For low J values the S 
branches belonging to the lower K, values will at first follow the 2(By)+C,) 
spacing but as J increases they will ‘‘break off’’ (lowest K, values first) 
gradually increasing and/or decreasing their spacing until they have a spacing 
of 4Cy as given in eq. (12). This behavior is shown (exaggerated) for the first 
six S branches in Fig. 3. 





Fic. 3. The behavior of the S branches for a slightly asymmetric top. (The scale for Avs 
has been exaggerated.) 
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D. ANALYSIS 


1. Determination of Rotational Constants 

The analysis of the observed spectra was carried out as follows. First an 
approximate spectrum was calculated using preliminary values of the rota- 
tional constants. (The constants reported by Romanko ef al. were used as the 
constants obtained by Allen and Plyler were not yet available.) The small 
effects due to centrifugal distortion were also included. These were calculated 
from the force constants given by Crawford, Lancaster, and Inskeep (1953). 
By means of these approximate spectra it was possible to assign almost all 
the lines in the C,H, spectrum and most of those in the C2D, spectrum. (Allen 
and Plyler’s constants were used to calculate the first C.D, spectrum as this 
analysis was carried out after their work was published.) Most of the observed 
lines were found to be blends, that is due to the superposition of lines arising 
from several transitions. However, 14 lines in the spectrum of C.H, and 9 
in that of C.D, were found to be due to single transitions. All were S branch 
lines. These lines were then used to determine the rotational constants making 
use of eq. (9). 

Equation (9) gives the displacement of S branch lines in terms of three 
unknowns, A, C, and x. Direct solution of this equation is not possible since 
x does not appear explicitly. Therefore the following procedure was adopted. 
For each of the single transitions the quantities [E,7+?(x) —E,7(«)] were cal- 
culated for various assumed values of «x. The equations for Av were then 
taken in pairs (each consisting of two different single transitions) to eliminate 
A. The resulting equations then contained C alone for the various values 
of x. This permitted the determination of a value for C as a function of x, 
consistent with the observed Av for each pair of lines. Graphs of these relations 
are given in Fig. 4 for CsH, and in Fig. 5 for C2D4. Ideally all curves should 
intersect at one point whose co-ordinates will be C and «. However, due to 
experimental errors the curves of Figs. 4 and 5 intersect at several points, 
although the intersections are in close proximity. Values of x were determined 
by averaging all intersections except those between curves that were almost 
parallel. The values are given in Table IV. Errors of twice the average devia- 
tion were assigned. These errors correspond to errors of 0.02 cm! in the 
measurement of less sensitive lines and of 0.04cm~! for the more sensitive 


TABLE IV 


Ground state rotational constants (cm~!) and asymmetry parameters of C2H, and C.D, 








Cal Ys 











C.H, 
Infrared, Infrared, 
Raman, Raman, Allen and Raman, Allen and 
present results Romanko et al. Plyler* present results Plyler* 
Ao 4.828+0.009 4.66+0.20t 2.432+0.008 
Bo 1.0012+0.0009 1.008; +0 .006 0.9998 0.7369 +0 .0012 0.7334 
Co 0.8282+0 .0004 0.828,+0.002; 0.8294 0.5630 +0 .0006 0.5636 
K —0.9135+0 .0003 —0.8139+0.0005 


*Allen and Plyler report errors of about +0.0002 cm~! for C2H4 constants and somewhat larger errors for C2Ds4. 
tDeduced from 1/Ao = 1/Co—1/Bo. 
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lines. With x evaluated it was then a simple procedure to recalculate the 
values of [E,-7+2(x) —E,7(x)] for the single transitions and to solve for Ao 
and Cy. The constant By was then calculated from x = (2B—A —C)/(A—C). 





0.9120 0.9130 0.9140 0.9150 
-« 








0.5650 

0.5630 1 $,(17), Sig 7)- 
T S9(i5), S5(8) 
E 
o 0.5610 
o 
© 05590 } 

| 
| $5 (16), Sg(I0) 
0.5570 |-—_+—_—_—+— ———} oe 


0.8125 0.8135 0.8145 0.8155 
—K 


Fic. 4. Determination of « for CoH. 
Fic. 5. Determination of « for C2D,. 


The rotational constants of C2:H4 and C.D, are given in Table IV. Large: 
relative errors are quoted for the constants of C:D, since fewer lines were 
available for the analysis of the C.D, spectrum. It will also be noted that 
the values of Cy quoted in Table IV differ slightly from those obtained from 
the graphs of Cy vs. « (Figs. 4 and 5) but are considered to be more accurate 
than the latter. For comparison, the rotational constants reported by Allen 
and Plyler (1958) and by Romanko, Feldman, Stansbury, and McKellar 
(1954) are included in Table IV. The final values of the calculated centrifugal 
distortion constants are given in Table V. 


TABLE V 


Calculated values* of the centrifugal distortion 
constants for CoH, and C2D, 





C2H,4, cm C2D,4, cm™ 
Dy 1.143 10-6 6.651107 
Dox 2.148X1075 5.822 1076 
Dr 7.960 X 1075 1.916X107-5 
R; —5.028 x 10-6 —1.634X 1078 
Rs —1.788X 1078 —2.249x 1078 
67 2.112107" 1.961 1077 


*The number of significant figures shown is not justified, but 
was used to secure internal consistency in the calculations. 
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2. Discussion of the Analysis and Its Limitations 

Since the values of the constants given in Table IV were obtained from 
an analysis of a small fraction of each of the observed spectra, it would be 
most desirable to confirm these values from the remainder of the spectra. 

An independent check of the value of Cy can be made according to eq. (12) 
from lines of low K, and high J. A graph of Avg vs. J should give a straight 
line (at high J) with a slope of 4Cy) and with an (extrapolated) intercept of 
6C)+3(2A —B—C)Ab—}(B—C)AC,. It is found that for J > 13 the points 
lie on a straight line. The straight line portions of the curves for C.H, and C.D, 
are shown in Figs. 6 and 7. These have been shown on a larger scale by removing 
some of the slope with a plot of Av—3n for C2H4 and Avy—2n for C.D. The 


15.5}-—-+ 
14.5 
13.5 


12.5 


Qv-3n cm-! 


1.5 





14 Ss i <2 22. 24 2 26 


Av-2n cm-'! 





14 16 18 20 22 24 26 28 30 32 34 
n 


Fic. 6. Av. vs. n for CoHy, where » = J for the S; and S2 series; nm = J—1 for the S; and 


S; series;... (In order to provide a larger scale 3n cm™ has been subtracted from each Av 
value.) 

Fic. 7. Av, vs. 2 for C2D4, where n = J for the S; and Se series; nm = J—1 for the S; and 
S, series;... (In order to provide a larger scale 2m cm™ has been subtracted from each Av 
value.) 


running number m on the abscissa is J for the S; and S2 series, /—1 for the 
S3 and S, series, and so on (see Tables I and II). From these graphs values 
of Co(C2H4) = 0.8279+0.0009 cm! and Co(C2D,4) = 0.5634+0.0007 cm=! 
were obtained, in excellent agreement with those in Table IV. Also good 
agreement was found for the intercepts: the observed values are 6.55+0.07 
cm-' and 4.49+0.03 cm~! for CoH, and C.D,, respectively, which may be 
compared with the calculated values of 6.60 cm! and 4.47 cm—'. 

In the earlier work of Romanko et al. (1954) an expression similar to eq. 
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(12) was used for the analysis of the S’ lines (that is the asymmetric top 
lines). It was derived from equations given by Mecke (1933) for the limiting 
oblate top. Nevertheless, it was shown to be a fairly good approximation for 
C.H, and led to a value of Cy = 0.828,+0.002; cm! in agreement with the 
value reported here within the quoted error. 

A final test of the validity of the constants given in Table IV was made 
by calculating the spectra and by comparing them with the observed spectra. 
The results are given in column 3 of Tables I and II. In column 4 of the 
tables are listed the approximate relative intensities of the calculated tran- 
sitions. These are based on the matrix elements of the polarizability for the 
symmetric top (Placzek and Teller 1933) with B = $(Bo+C») and on the 
statistical weights given in Section C(2). For C.Hy, the over-all agreement 
between observed and calculated Av is good. In particular for those lines which 
were subsequently found to be due to single transitions, and were not used 
in the analysis, the agreement is within the experimental error of +0.03 cm —!. 
Furthermore, it was possible to assign all but four of the lines in the CoH, 
spectrum. An indication of the good agreement can also be obtained from 
Fig. 2 where the calculated spectrum is drawn beneath the microphotometer 
trace. For C2D4, however, the over-all agreement is not as good and more- 
over many of the lines remain unassigned. The reasons for this are discussed 
below. 

As just mentioned, all but four of the lines in the spectrum of C,H, were 
assigned with the constants of Table [V and with the selection rules AJ = 1, 2; 
AK, = 0. One of the unassigned lines, however, might possibly arise from a 
transition with AK, = 2. This line is weak and sharp (see Fig. 2) and was 
measured on three plates with an average Av = 33.97;cm™. It lies at the 
calculated position (Av = 33.980 cm!) of the transition 735-717, and in 
this region all other lines of appreciable intensity have been accounted for. All 
of the other transitions with K, = 3<1 (as well as Kg = 2<0) led to 
calculated Av very close to R, S, or ghost lines. Similarly, transitions with 
AK, = 2 amongst the higher K, levels gave lines in crowded regions of the 
spectrum and unfortunately none were detected. Such lines would be expected 
to be less intense because of the Boltzman factor. If the assignment of the 
above line could be established as 73; <— 717, it would serve as a critical check 
on the value of Ao since the Av for transitions with AK, = 2 are mainly 
dependent on Ao. 

The most serious limitation in the preceding analysis is probably the fact 
that only the 11 lowest r levels* (corresponding to Kg < 5) could be calcu- 
lated. Moreover, these were restricted to J < 27 for C.Hy and J < 19 for 
C.D,. The next 7 levels and therefore the next S branches would correspond 
to K, >6. These branches are undoubtedly present in the spectrum of 
C.D, and to a lesser extent in that of C2H;, but with less intensity than 

*In Golden’s paper the second-order corrections (C2) are given for only the eight lowest 
r levels. The next three 7 levels were calculated with the assumption C, = 0. This approxima- 
tion did not produce serious difficulties since it appears that the correction involving C2 is 


almost the same for all energy levels giving rise to a specific S branch and therefore contributes 
very little to any given Ap. 
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the first 11 S branches. As pointed out before, the energy levels of intermediate 
and high 7 levels are expected to show little deviation from the limiting sym- 
metric top levels (for small asymmetry), and transitions between these levels 
would give rise to a series of broad lines with a spacing of 2(Bo+ Cp), according 
to eq. (10). Such a series of lines was indeed found in the C2H, spectrum 
and it is labelled as Sx in Table I and Fig. 2. Some lines of the corresponding 
Rx series were also evident. While some of the lines in the C.D, spectrum 
can be assigned as Sx lines, their identification is not unambiguous and the 
series cannot be followed to the end of the spectrum. This was not unexpected 
in view of the larger asymmetry of the C.D, molecule. There is no doubt 
that the contribution of transitions amongst the 7 levels higher than 11 is 
mainly responsible for the large number of unassigned, and apparently ‘“‘mis- 
fitted’’, lines in the C.D, spectrum. It is also quite possible that some of the 
lines chosen as single transitions are overlapped by weak lines arising from 
these 7 levels and again this occurrence would be more likely for C.D,. 
Finally it should be mentioned that a few of the lines selected as single 
transitions on the basis of the approximate spectrum were found to be over- 
lapped by weaker R or S branch lines as indicated in Tables | and II. These 
weaker lines were so sensitive to the changes in rotational constants that 
their final calculated positions were shifted, by appreciable amounts, to 
within +0.2 cm~! of the lines used in the analysis. However, the neglect 
of these few lines would not alter the analysis or constants of Table IV to 


any appreciable extent. 


E. STRUCTURE OF THE ETHYLENE MOLECULE 
The principal moments of inertia of CoH, and C.D, and the inertial defects 
A = I¢—I4—T are given in Table VI. (The small positive values of the 
inertial defects indicate the planar structure of the molecules.) The three 


TABLE VI 


Ground state moments of inertia and inertial defects for ethylene and 
ethylene-d,, and the structural parameters of the ethylene molecule 











Constant CH, C.D; 
i 3.493-40.007 a.m.u. A? 6.934+0.023 a.m.u. A? 
Ip 16.843+0.015 22 .884+0.037 
It 20.361+0.010 29 .952+0.032 
A 0.025+0.020 0.134+0.054 
ro(C=—60) ro(C—H) H—C—H 
Raman 1.339+0.002 A 1.086+0.003 A 117°34’ +30’ 
Infrared 1.337+0.003 1.086+0.003 117°22’+1° 
Electron diffraction 1.334+0.003 1.085+0.005 116°+1° 


structural parameters of ethylene can be determined from the present data 
by making the assumption that the structures of CoH, and C2D, are identical 
in the v = 0 vibrational levels. This assumption of course neglects the small 
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effects of zero-point vibrations. Only three independent moments of inertia 
are necessary for this determination and the values 12(C2H,), I¢(C2D,), and 
the average of the two values 74 were used. (The values of J; were not 
used in the present work in view of their much larger errors.) The resulting 
parameters are given in Table VI. For comparison, the recent values obtained 
by Bartell and Bonham (1957) by electron diffraction and those obtained 
by Allen and Plyler (1958) by infrared spectroscopy are included in this 
table. 

The three sets of parameters for ethylene as determined by these three 
different experiments are essentially the same. The three sets appear to have 
the same accuracy and are considerably more accurate than all the previous 
determinations of the ethylene structure. All three parameters differ signi- 
ficantly from earlier values with perhaps the change in HCH angle from 
120° to 117° being the most striking. 

It would appear from Table V that the Raman and infrared values are 
in closer agreement; however, the apparent small difference between the 
spectroscopic and electron diffraction values may not be significant because 
of the effects of zero-point vibrations. Of course these effects enter into both 
types of measurements, but they do so in different ways and could result 
in small differences in the structures. An attempt to include these effects 
was made in the case of the electron diffraction results, but this is not yet 
possible for the spectroscopic results. 

In the spectroscopic determination of molecular structures from ground 
state constants of isotopic molecules, it is well known that zero-point vibration 
effects may yield several slightly different sets of values for the parameters, 
depending on the isotopic molecules used in the calculation. Recently, Costain 
(1958) has demonstrated that a consistent set of parameters can be obtained 
from ground state constants by using isotopic substitution on all non-equi- 
valent atoms. This leads to the so-called ‘‘r, structure’ which appears to 
approach more nearly the equilibrium structure than does the “rp structure”’ 
obtained from the usual calculations. In principle ry structures are determined 
directly from the experimentally determined moments of inertia, whereas r, 
structures are determined from differences between moments of inertia. The 
determination of the “‘r, structure’ of ethylene requires the substitution of 
at least one C'? atom by a C™ atom, as well as the H—D substitution.* The 
present data therefore only allow the calculation of the H—H distances in 
ethylene: however, these data do permit the calculation of r,(C—-H) and 
ZHCH for various assumed values of 7,(C==C). The results are given in 
Table VII. Clearly the C—H parameters are insensitive to the assumed range 
of C=C bond lengths. This perhaps indicates little or no real discrepancy 
between the spectroscopic and electron diffraction results. 


*The next best procedure would be to study the spectrum of Dx.C = CHe (C2, symmetry) 
and to make use of the first moment equation to determine the shift in the center of mass 
along the C2 axis and to evaluate the C=C bond length. 





5 
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TABLE VII 
The “‘r,”’ structure of ethylene 


r.(C=C) Soe: r, (C—H), 








Z£H—C—H 
1.335 1.085 116°50’ 
1.337 1.084 116°54’ 


1.339 1.084 117°0’ 
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PHOTO ALPHA REACTIONS IN SILVER! 


J. P. Roatsvic, R. N. H. Hastam, L. D. SKARSGARD,? AND 
E. E. WuscHKeE? 


ABSTRACT 
The reaction Ag'(y, a) Rh! has been investigated by using the secondary 
activity of Rh, Statistical theory was used to obtain information about the 
reactions Ag!7(y, a) Rh" and Ag(y, a)Rh. The present values for the absolute 
yield and cross section are in agreement with previous work using either the 
secondary activity or photographic plates. 
I. INTRODUCTION 
Photo alpha reactions in silver have been investigated both by measurement 
of the secondary activity of the daughter nucleus (Greenberg ef al. 1954; 
Laboulaye and Beudon 1954; Schmouker e¢ al. 1955; Erdés et al. 1955; Erdés 
1957) and by photographic plate work (Greenberg 1957; Toms et a/. 1958). 
The fact that the yield values, and consequently also the cross-section values, 
measured previously at this laboratory by the secondary activity method 
(Greenberg et al. 1954), are much lower than values obtained at other labora- 
tories by the same method and also all values obtained from plate work 
suggests that our previous values are in error. It is probable that the chemical 
separation used to isolate the radioactive Rh'® from the irradiated silver was 
unsatisfactory. It was thus thought worth while to re-examine the reaction 
Ag!(y, a)Rh'” using improved techniques. Another factor making this ex- 
periment desirable is the recent redetermination of the activation curve for 
the reaction Ni*(y, 2) Ni” (Roalsvig, Haslam, and McKenzie 1959), which 
was used as a standard for normalization of the rhodium yield. 


Il. EXPERIMENTAL TECHNIQUE 

1. Irradiation 

The samples, consisting of 80 g of AgNOs, were irradiated in a cylindrical 
lucite container, 5 cm long and 3 cm in diameter. The cylinder was placed in 
a special holder close to the target of the betatron donut, so that the X-ray 
beam passed along its axis. Thin nickel disks of the same thickness and 
diameter were used as monitors and were fastened to the front and back of the 
container. The samples were irradiated at different maximum photon energies 
of the betatron, ranging from 15.5 Mev to 24.0 Mev. The irradiation time 
ranged from 3 hours at 15.5 Mev to 45 minutes at 24.0 Mev. 


2. Separation Method 

At the end of the irradiation period, the rhodium formed in the reaction 
Ag(y, a)Rh was separated from the silver nitrate by the following chemical 
procedure: 
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The silver nitrate was dissolved in 100 ml of distilled water and 10 ml of 
Fe*+ carrier solution (containing 9.64 mg of Fe*+ per ml) were added. This 
was then precipitated with 10 ml of NH,OH, filtered, and washed with dis- 
tilled water. After dissolution of the precipitate in HNO, a similar separation 
was carried out again, adding another 10 ml of Fe** carrier and precipitating 
and washing with NH,OH. After filtration, the precipitate was again dis- 
solved in HNO; and a similar procedure carried out for the third time, this 
time adding 20 ml of Fe*+ and precipitating with 50 ml of HN,OH. After 
filtration, the final precipitate was dried and ignited to Fe.O3. 


3. Counting Procedure 

The final precipitates were counted between two Geiger—Miiller counters 
placed end to end to secure higher counting efficiency. The nickel monitor 
disks were also counted, the front one before and the back one after each 
rhodium count. The counting rates were corrected to the mid-point of the rho- 
dium counting interval by application of the appropriate exponential correc- 
tion for decay. The samples were followed through several half-lives. 


4. Correction Factors 

(a) Geometry and Backscattering 

The samples were all counted in the same geometry, so no correction for 
different geometry factors was necessary. Also, the samples were mounted on 
a very thin Q-film, making corrections for backscattering negligible. 

(b) Self-absorption 

The observed counting rates were corrected for self-absorption in the 
samples, using the formula for self-absorption of Baker and Katz (1953). 

(c) Half-lives 

If the half-lives of Rh! and Ni*” are the same, the ratio of rhodium to 
nickel counts should remain constant for samples irradiated together. The 
ratio was measured through several half-lives, and was found to be constant 
as expected. 

(d) Efficiency of Powder versus Disk 

The counting efficiency of powder is lower than that of the disks, that is, 
if equal masses of nickel powder and nickel disk are given the same dose in 
identical geometry, the observed activity will be higher for the disk than for 
the powder, counted in identical geometry. A test of the relative counting 
efficiency of nickel powder and nickel disks was made, and the ratio obtained 
was used in the analysis of the results to correct for this effect. 

(e) Relative Doses 

Method I: Samples of AgNO, into each of which had been uniformly 
mixed 1 g of nickel powder, were given 15-minute irradiations in the same 
position as was used for the other samples. Nickel disks were again mounted 
in front of and behind the samples. After irradiation the nickel powder was 
removed by dissolving the sample in water and mechanically filtering out the 
nickel. A mass of powder equal to the mass of each of the front and rear disks 
was weighed out, and alternate counts were taken on the activities of the 
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powder and the disks. After correcting for the effect (d), the dose given to the 
powder, and consequently to the silver nitrate, was related to the dose given 
to the two monitor disks. 

Method II: Measurements of the ratio of the activities in the front and rear 
nickel detectors showed that this ratio remained constant over the range of 
energies used in this work. That indicates that the effect of beam spreading 
with energy is largely balanced by changes in the beam absorption with 
energy. The beam was therefore assumed to be absorbed in an exponential 
manner in passing through the sample. In this case N, = Nye~*’” where 


N», = dose given to back detector, 
N; = dose given to front detector, 


, 


be 
i 


Experimental values gave yu’ = 0.168 cm~'. The average dose given to the 
AgNO; sample was then calculated from 

aL 
J Nie “dx 


0 
N av eS ono ae 


L 
The relative doses, as measured by methods I and II, were in good agreement. 


(f) Chemical Separation 

Correction should also be made for the percentage recovery of iron in the 
chemical separation; 40 ml of Fe(NO;)3 solution were introduced, containing 
9.64 mg of Fe** per milliliter. Thus if all iron were recovered, the mass of iron 
should be 386 mg. If the precipitate is in the form of Fe.O; after oxidation, 
its mass should be 552 mg. On the other hand, if the iron is present as FeO, 
the mass of the precipitate should be 497 mg. The actual masses were found 
to be between 510 and 550 mg. Thus the calculated percentage yield range 
from a minimum of 93% to a maximum of 110% depending on the assumption 
made about the form in which the iron is present. The precipitates did differ 
in color, indicating incomplete oxidation in some cases. Thus it is impossible 
to determine the actual chemical yield, and the best assumption seems to be 
to take 100% as the yield. 

(g) Isotopic Abundance 

The counting rates were corrected for the abundance of the isotopes Ag'’® 
and Ni*®. 


effective linear absorption coefficient, 


length of sample. 





III. RESULTS 

1. Activation Curve 

After application of the above-mentioned correction factors, the ratio of the 
yield in a-particles per roentgen per mole Ag!” to the yield in neutrons per 
roentgen per mole Ni®® was found as a function of the maximum betatron 
energy. In order to find the absolute yield values for the Ag!(y, a)Rh'™ 
reaction, the absolute yield values for Ni *(y, 2)Ni57 must be known. The 
activation curve for the latter reaction has recently been redetermined at this 
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laboratory (Roalsvig, Haslam, and McKenzie 1959). Using the values taken 
from this work, the absolute activation curve for Ag!(y, a)Rh!® was calcu- 
lated and is given in Fig. 1. The absolute yield value at 23 Mev was found to 


be 970 a/r mole. 


r permole 


Saturated activity, « particles per 





Fic. 1. Absolute 
values of Laboulaye 


32 Mev. 


2. Cross-section Curve 
The cross-section curve was calculated from the activation curve by using 
the method of Penfold and Leiss (1958), and is given in Fig. 2 (solid curve). 


Fic. 2. 
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The characteristic values are: 
oy, =SS KD 
En = 21.0 Mev 
Tr; = 6.1 Mev 


224 Mev 
{ odE = 0.31 Mev.mb. 
Jo 
3. The Reaction Ag!” (vy, a) Rh! 

This reaction cannot be measured by secondary activity since the daughter 
nucleus Rh'® is stable. However, by using the statistical theory, the cross- 
section curve for this reaction might be determined when the cross-section 
curve for the reaction Ag!(y, a)Rh!® is known. 

For two neighboring isotopes such as Ag!” and Ag!”, the cross sections for 
y-capture are almost equal. Also the sums > F; of the half-widths covering all 
modes of decay of the excited nuclei must be equal. We thus have (Blatt and 
Weisskopf 1952, Chapter VIII): 


aya) _ Fay _ Jan op LK f- "KlSeailt. 
a(y, &) (2) Fa¢2) Ta(2) E,=0 
where 
K =aconstant 
w(e) = the level density of the residual nucleus at an energy e, where ¢ = 
hv—Q-E, 
Q- = the binding energy of the alpha particle 


Sac = the capture cross section for the alpha particle of energy /,. by the 


residual nucleus. 
The ratio {o(y, a)Ag!"}/{o(y, a)Ag'®} has been calculated by Greenberg 


(1957). It decreases fairly linearly from 2.5 at 16 Mev to 1.3 at 24 Mev. The 
calculated cross-section curve for the reaction Ag!""(y, a) Rh! is given in Fig. 
2 (dotted curve). The characteristic values are: 

om = 94 ub 

En = 20.5 Mev 

Ty = 6.3 Mev 


224 Mev 
| odE = 0.54 Mev.mb. 


0 


4. The Total Reaction Ag(y, «)Rh 
The cross section for the total reaction Ag(y, a)Rh was calculated from 
the cross-section curves of the two isotopes. It is also shown in Fig. 2 (dashed 
curve). The characteristic values of this curve are: 
Tn = 15 pb 
En, = 20.8 Mev 
Ty = 6.3 Mev 


24 Mev 
f odE = 0.43 Mev.mb. 
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IV. DISCUSSION 


1. Ag!®(y, a) Rh!% 

The previously reported values for this reaction are given in Table I. 
Laboulaye and Beudon (1954) normalized their relative values according to 
the reported values of Katz et al. (1951) for the reaction Ni®(y, 2)Ni*”. If 
the newer values of Roalsvig, Haslam, and McKenzie (1959) are used for 
this reaction, the yield value of Laboulaye and Beudon at 23 Mev becomes 
1200 a/r mole, which is in reasonable agreement with our value 970 a/r 
mole. The yield value of Greenberg et al. (1954) is as mentioned previously 
far lower than any other value reported. However, the same author (1955) 
obtained a new yield value which is somewhat higher than the present one. 
Schmouker e¢ al. (1955), Erdés et al. (1955), and Erdés (1957) obtained a yield 
of 1360 a/r mole at 32 Mev, which agrees with our value in the order of magni- 
tude, as can be seen by the dashed extension of our activation curve on Fig. 1. 
The present value for the maximum cross section is much lower than that 
reported by Laboulaye and Beudon (1954), by a factor of about two if the 
renormalized curve is used. 

However, our value for the integrated cross section lies in the region between 
the value of Laboulaye and Beudon and that of the Swiss group. 


2. Ag(y, a)Rh 

The previously reported values for the gamma-—alpha reaction in silver are 
given in Table II. By using the photon spectrum as given by Katz and 
Cameron (1951), the ratio of the yields of the reactions Ag(y, a)Rh and Ag! 
(y, a)Rh'® at 23 Mev was calculated from the curves in Fig. 2 to be 1.46. 
Taking the yield value of Ag!®(y, a)Rh!® to be 970 a/r mole, the yield value 
of Ag(y, a) Rh at 23 Mev becomes 1420 a/r mole. 

The present yield value at 23 Mev is somewhat lower than those reported 
by Greenberg (1957) and Toms et a/. (1958) using photographic plates, how- 
ever, of the same order of magnitude. Since events in which neutron or proton 
emission follow the alpha-particle emission would be detected in plate work, 
but not in our work, the present lower value would be expected. 

This might also be seen by a comparison of the cross-section curves of 
Greenberg (1957) and of the present work. Greenberg’s cross-section curve 
is still rising at 24 Mev, while our curve has passed its maximum and is falling 
off rapidly at this energy. 
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STATISTICAL THEORY OF AVALANCHE BREAKDOWN IN 
SILICON! 


R. E. BuRGEss 


ABSTRACT 


The current flowing in silicon junctions under conditions of high electric field 
sufficient to cause avalanche breakdown consists of a sequence of pulses of 
constant magnitude and with randomly distributed spacing and duration. A 
statistical model is presented for this intermittent conduction mechanism in 
which the transition probabilities between the conducting and non-conducting 
states are Markoffian and depend on the excess of the applied voltage over the 
breakdown voltage; power-law and exponential dependences are considered and 
the precise form is shown not to be critical in deriving from these simple con- 
ditions a number of predictions about the macroscopic current-voltage character- 
istic, all of which have been confirmed by experiment. 

It is demonstrated generally that a randomly switched discharge of positive 
resistance has a macroscopic current-voltage characteristic which depends upon 
the d-c. source resistance and can have a region of negative slope. Only with a 
negligibly small source resistance is the true characteristic obtained. 

The observable fluctuations of current and their spectral distribution are 
derived and indicate the noisiness of silicon junctions under avalanche con- 
ditions. 

The influence of illumination on the breakdown statistics and macroscopic 
characteristic are briefly discussed. 


1. INTRODUCTION 


Conventional circuit elements whether linear or non-linear display d-c. 
current-voltage characteristics which are independent of the resistance of the 
e.m.f. source used for the measurement. If, however, the element is time- 
variable (either in a systematic or random fashion), the macroscopic d-c. 
characteristic, i.e. the relation between the average current and average 
voltage, will in general be a function of the source resistance. This implies 
that a detailed specification of external source conditions is required in order 
to give an unequivocal meaning to measured characteristics. 

The present paper is primarily concerned with the properties of avalanche 
diodes which effectively contain a random switch in series with a resistive 
component which may be linear or non-linear but always has a positive 
differential resistance. 

In general if the source e.m.f. is E and its resistance S the mean current 
7 and the mean potential difference (p.d.) 6 across the diode are related by 
E = iS+i where the bar denotes averages over times long enough to provide 
representative values of the current and voltage, i.e. long compared with 
the relaxation times of the process. Then 7 and @ correspond to values 
indicated by meters having inertia. 

The ratio 3/7 is the apparent d-c. resistance Ry of the diode, and the deriva- 
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tive di/di is its apparent differential resistance p. Both of these will, in general, 
be a function of the source resistance and e.m.f. Furthermore, even though 
the resistive components in the circuit are positive, it is possible for the 
apparent differential resistance to be negative under certain conditions. 

These considerations are particularly relevant to the phenomenon of 
avalanche breakdown in strong electric fields (over 210° v/cm) arising in 
silicon avalanche diodes in which the current flow is composed of a random 
sequence of rectangular and unidirectional current pulses of fixed amplitude. 
The random timing and duration of the impulses correspond to the stochastic 
character of the avalanche multiplication process which causes a breakdown 
discharge in localized regions of high electric field (McKay 1954; Newman 
1955; Chynoweth and McKay 1956). The use of these diodes as voltage regu- 
lators and switching elements calls for an examination of the manner in 
which the random character of the conduction process determines the macro- 
scopic characteristics in a given circuit and also the ‘‘noisiness’’ of the device 
as expressed in terms of microscopic fluctuations of the instantaneous operating 
point about its mean position. 

Some of the experimentally observable effects in avalanche diodes are 
shown to stem from the voltage dependence of the switching transition 
probabilities from one state to the other. Physical models are related to these 
statistical parameters and lead to a phenomenological explanation of the 
complicated characteristics observed in avalanche diodes under breakdown 
conditions. The basic statistical model analyzed involves a switch in two 
states—open (state 0) and closed (state 1) and undergoing transitions between 
these states with transition probabilities po; and pio per unit time. In the 
steady state (assuming statistical stationarity) the fraction f of the time that 
the switch is closed is given by 


_f _ bu 
(1.1) race dad 
or 
, — os 
44) f= PotPpio 


The mean frequency of the switch closings (or openings) is 


(1.3) y = Se LL 


where 7; = 1/Pio is the mean duration of the 1-state and 79 = 1/0: is the 
mean duration of the 0-state. 

For the current pulses to be rectangular (i.e. very rapid rise and fall and 
flat-topped) the capacitance across the diode must be small enough for the 
associated charge and discharge time constants to be small compared with 
1/pio and 1/fo1. 

The d-c. macroscopic characteristic relating 7 and 3 as measured by meters 
with time constants larger than 1/v also represents the a-c. characteristic at 
frequencies small compared with v. The irregular variations of the current 
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and voltage will appear as noise superimposed on the average characteristic. 
The differential resistance discussed later is defined as the slope of the 3, 7 
curve at constant S: 


_(@) _(@) /(a 
(1.4) =) 7. 


This definition is chosen because of its theoretical significance and its corre- 
spondence to the manner in which the characteristic is commonly determined, 
i.e. by variation of E with S fixed. 

If 79 and vp are the current and voltage in state 0 and 7; and v,; the values in 
state 1 the macroscopic averages are 

t= fix t(1—f)to b= E-—S[fiit+(1—f)io} 
and the apparent d-c. resistance is 
E 
(1.5) Ry = pS, 
fir+U—f)to 


while the differential resistance is 








1 
1.6) 2S -—§ 
( ° = 3in,+ —\yel dE 
Even in the simple case where the random switch is in series with a fixed 
linear resistance R so that 7) = 0 and 7; = E/(S+R) and even if f is constant 


it is seen that both the d-c. resistance and the differential resistance depend 
on the resistance of the source. Furthermore, if df/dE is sufficiently large, p can 
be negative even though the system contains a linear positive resistance. 
These features are unorthodox since normally the d-c. characteristics of even 
a non-linear resistive device do not depend upon the source resistance used 
for measurement. 


2. CHARACTERISTICS OF A RANDOM AVALANCHE BREAKDOWN 

Let it now be supposed that f and R are both functions of voltage as is the 
case in the intermittent discharge which appears in silicon p—n junctions under 
avalanche conditions. In general the statistics of the switch depend upon the 
voltage across the device in both states. It will be assumed that the diode 
passes a prebreakdown current 7) when the applied p.d. is less than threshold 
value v; for v above v, the current increases proportionally to the excess 
voltage due to the discharge path: 
(2.1) t=% forv<y l= io tg(v—v%) forv > 


where g is the discharge conductance. Then when E > 7+Sio 


5 a Et St— Sin . _ . , g(E—m— Sin) 
7 1+ Sg ae 1+S¢g 


are the p.d. and current in the discharge condition (state 1). 


7 g(E-—% — Sto) 


Hence 
fg(E—v, — Sto) 


9« a N= oe - 
(2.3) l E-Si v tot i+ Sg 
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On the assumption that the parameters g and 7 are independent of E the 
apparent differential resistance is 
(2.4) => le nian 

d[f(E—— Sio)]/dE 
which is a function of both E and S. The value of p when S = 0 will be denoted 
po (the intrinsic p of the diode) and this is seen to be always positive for the 
systems under consideration in which df/dE is non-negative. 

When S > 0, the 7, 6 characteristic can have a negative slope if f is a suffi- 
ciently rapidly increasing function of E and this tendency is more pronounced 
as S is made larger. Clearly a negative differential resistance can only appear 
if S > po/(1—pog). Hence the expected behavior is that, if S is sufficiently 
large, as E (and 7) increase d will show first a maximum and then a minimum 
with the region of negative p lying between these extrema. 

In terms of the transition probabilities we have 


$ <9 
dE Po. dv» 1+gS Pio dv, 
where po; is a function of vp (=E—Sio) while pio is a function of m. It is seen 
that df/dE is zero at f = 0 and | and will have a maximum at some value of 
f which depends upon the precise dependence of foi1/f10 on (E—v) and will 
usually lie near f = 0.5. 

The intrinsic macroscopic characteristic has a slope, 


1 = d{f(E—»)] 


(2.5) oe ee 





which can be appreciably greater than the true conductance g of the avalanche 
discharge path. This intrinsic parameter enables the differential resistance 
for any finite S to be expressed as 

(2.6) p = po(l+g5S)—S. 


In view of the statistical nature of the avalanche process it is clear that 
Po. must be proportional to the rate of influx of carriers entering the high- 
field region where breakdown occurs and also be an increasing function of the 
excess of vp over the breakdown voltage v,. Conversely the switch-off proba- 
bility pio will be a rapidly decreasing function of (v1 —%), since the discharge 
can be maintained with greater certainty when there is an excess of energy 
imparted to the avalanching carriers over that needed just to sustain the 
feedback mechanism provided by the ionization due to the carriers of both 
signs. 

When the discharge region is broken down and taking up a potential differ- 
ence v the resistance (1/g) of the system is assumed to be constant. This 
resistance is composed of the sum of the differential resistance of the dis- 
charge and the series resistance of the spreading regions on either side of the 
discharge; the current constriction due to the extremely localized nature of 
the discharge will probably result in this bulk resistance being the significant 
part of the total. 
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(a) Power-Law Model 

In order to consider the dependence of f on E in more detail we must postu- 
late the dependence of the transition probabilities po; and p19 on the voltages 
vo and v, which exist before each transition. One possible form suggested by 
physical considerations is: 


JS Por = A(vo—%)” 





(2.7) 
Pro = Bivi—%)” 
giving 
1 eee 
£ ~ A (E—vp— Sio)”™ 


which assert that the probability of breakdown (0 — 1) is an increasing func- 
tion of the excess of vp over v% while the probability of the discharge ceasing 
(1 — 0) is a decreasing function of the excess of v; over v%. Thus f is a mono- 
tonically increasing function of £. 

For this model it is found that 


df _ f[U-f)(m+n) 
dE E—n%— Sig 


and hence 


as » = ESL L— fon) 141 

— gf{1+ —f)(n+m)] 

If gS is sufficiently large the macroscopic characteristic will show ‘turnover’ 
(p = 0) at the two values of f given by 

(2.9) p= Mm ltl (mtm—1)'—4(n-+m) /gS)" 

oat ; 2(n+m) : 

which correspond to the maximum and minimum of 3 which appear in the 
macroscopic characteristic. 

The minimum value of po and hence of p for any arbitrary S occurs when 

f = (n+m-+1)/2(n+m) which is midway between the two values for turn- 
over (if they exist). At this point 


_ 4(n+m)—gS(n+-m— 1)*. 





9 = = > 
(2.10) Pp Pmin g(ntm+1) 
Thus if 

: 4(n+m) 
9 
a 3? c(n-+m—1) 


a negative slope will appear in the characteristic which would not be mani- 
fested if S were zero. The slope of the intrinsic characteristic is 


(2.12) . = gf[1t+(m-+n)(1—f)]. 


Experimental data on ionization rates as a function of electric field and 
on the slope of intrinsic characteristic indicates that m-++-n>>1. Thus the value 
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of f for which the slope is a maximum and hence p is a minimum (for any S) 
is just slightly above 3. Experimental observation confirms that the slope 
of the macroscopic characteristic for S = 0 can, when f ~ 3, considerably 
exceed the conductance g corresponding to a sustained discharge (f = 1). 
For instance, if typically m+n = 14 and f = 3 the slope is 4g. Thus the 
apparent sharpness of the breakdown process is enhanced by the relatively 
rapid increase of the duty-cycle ratio f with applied voltage. For m+-n = 14 
the value of f for minimum 9 is 0.536 and a value of S of only 1/3g would be 
adequate to cause the appearance of a negative slope in the 7, 6 characteristic; 
a smaller value of S would preserve a positive slope which, however, would 
be very much greater than g. The importance of the source resistance in 
determining the abruptness of the breakdown in avalanche diodes is over- 
looked in commercial specifications of these devices. 

From equations (1.2), (1.3), and (2.6) it is seen that as the source resistance 
S is increased the mean frequency v of the pulses increases while the mean 
duration 7; and their fractional duration (f = v7,) decrease. As the source 
voltage is increased r, and f increase while v passes through a maximum. 


(b) Exponential Model 

An alternative model for the transition probabilities has the general form 
(2.13) foun = a expla(o—%)’] 

lpr = bexp[—B(u1—%)’] 
where r can be positive or negative. It is seen that the following conditions 
would hold for the parameters: 


ifr >0O a,B>O a<b 
ifr< 0 a, 2 <9 a>b 
These are required to ensure that fo: is an increasing function of vo and that 
pw is a decreasing function of 2. Also if r is positive, po1/p10 must become 
very small compared with unity at voltages just above v,; conversely, if r 
is negative fo1/Pi0 must become very large compared with unity when the 
voltage appreciably exceeds vp. 
The slope of the intrinsic characteristic for these transition probabilities 


is found to be 


: eet oy 
(2.14) - gf +c f)r in (tL 


Po 
and the important term in 7 In(b/a) is seen always to be positive. If, as is 
anticipated, this term is large compared with unity the slope will have its 
maximum for f close to but slightly above 3. At f = 3 
7 1 e.0 
(2.15) -—= & i+¢ la- |. 
Po o “ a 


Typically if r = 1, b/a = 10° the value of 1/po is 4g. It is seen by compari- 
son with the power-law model that the intrinsic slope is closely similar in form 
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if (m+n) is replaced by r In(b/a). It is concluded that (except when f is very 
near to 0 or 1) the details of the model for the transition probabilities are not 
significant in determining the macroscopic characteristics of the device. 

The choice of r would depend upon physical considerations such as the 
exp(—c/F) dependence of the ionization rate on electric field F (cf. Chyno- 
weth 1958) which determines the probabilities of the breakdown initiating 
and ceasing. The use of a negative r may be particularly appropriate where the 
probabilities need to be described accurately in the region of small excess 
voltage (po: very small, p10 very large, f very small). 


3. FLUCTUATIONS 
The variances of the current and potential difference provide measures of 
the over-all ‘‘noisiness’’ of the device: 


(3.1) var i = f(1—f)(i;—%0)?, 
(3.2) var v = S* var 1. 


Both variances have f(1—f) as a factor which arises from the similarity of the 
random switching process to a Bernoulli trial. 
The spectral density of the current fluctuations is given by (see Appendix) 


4rvari — 4f(1—f)(i:—io)*r 


3.3 W,(o) = ——33 = Lege 
( ) i( ) I+w’r’ 1+w'r’ 
where +t = 1/(foitpPio). This simple type of single-relaxation-time spectrum 
had already been inferred for silicon avalanche noise (Burgess 1955). 

Substituting from (2.2) 

; 4f(1—f)r| g(E—v— Sio) |? 
(3.4) Wile ea 
l+w'r 

Although the spectral density has the simple form (1+ w?r?)—! for its fre- 
quency distribution, its dependence on the source e.m.f. E and resistance S 
is complicated since these terms appear directly and also implicitly in f and r. 

Since the noise arises from a gross modulation of the current flow it is 
normally very much greater than the shot noise 2e7 which would be associated 
with the current, the ratio being of the order of (7;—%)7/e which might 
typically lie in the range 10* to 10". Alternatively the noise due to the ava- 
lanche process may be likened to a form of shot noise in which charge gq of 
mean square value 2(7;—7o)?7? is randomly transported at mean rate v so giving 
a low frequency spectral density 2vg? = 4f(1 —f)7r(7;— 7%)? in agreement with 
(3.3). 

4. EFFECT OF ILLUMINATION ON AVALANCHE BREAKDOWN 
If a silicon junction is illuminated with photons of energy sufficient to 


produce carrier pairs, a fraction of these will enter the localized region in 
which avalanche breakdown is occurring. The dominant effect of this will be 


————— 
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to increase fo; in direct proportion to the increase of prebreakdown current 1% 
since each carrier entering the region has a certain probability of initiating 
the discharge, that is, the parameter A in equation (2.7) or a in (2.13) is 
proportional to the carrier influx. The other parameters are not expected to 
be affected significantly; some increase of g due to bulk photoconductivity 
may become important at large values of S. 

The observable effects of illumination are thus expected to be the increased 
frequency of the current pulses by virtue of a decrease of the mean off-time 
to. The resulting increase of f increases 7 whatever the value of the source 
resistance. However, the slope of the intrinsic characteristic may be increased 
or decreased by illumination depending upon whether f is moved closer to or 
further from the value which minimizes po. 


5. CONCLUSIONS 


The statistical treatment of the previous sections leads to the following 
conclusions regarding silicon junction diodes in the avalanche breakdown 
region: 

(i) The macroscopic characteristic is a function of the source resistance and 
as (E—yp) increases the characteristics for various S converge to a single func- 
tion since f — 1 and the uninterrupted discharge is measured. 

(ii) If the source resistance is large enough a negative slope will appear 
between two voltage turnover points in the macroscopic characteristic. The 
minimum value of the slope resistance occurs for a value of f slightly above 
0.5. 

(iii) The mean frequency of the current pulses is an increasing function of 
the source resistance while it displays a maximum as E— yp increases. 

(iv) The amplitude of the current pulses is a linear function of the applied 
voltage (equation 2.2). 

(v) The mean duration 7; of the impulses and the fractional duty cycle f 
are increasing functions of (E—vz) and decreasing functions of S. 

(vi) The slope of the intrinsic characteristic can be appreciably greater 
than the true conductance of the avalanche discharge path. 

(vii) The noise spectral density of the current fluctuations is of the form 
W(0)/(1+w?r?), and is greatly in excess of shot noise (W(0) > 2e7). 

(viii) Illumination of the junction decreases the mean off-time zo and in- 
creases the mean pulse frequency v and the fractional duty cycle f. 

All of these inferences have been confirmed experimentally (Drews 1957). 
They also incidentally indicate that most published experimental data and 
commercial specifications of silicon avalanche diodes are without significance 
in the absence of a statement of the resistance of the voltage source used for 
the measurements. In principle this should be zero; when S is made very small 
there is experimental evidence to indicate that no negative resistance arises 
in the discharge path itself and this is consistent with our adoption of a 


positive g in the model. 








738 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 
REFERENCES 
BurGeEss, R. E. 1955. Brit. J. Appl. Phys. 6, 185. 
—— 1956. Proc. Phys. Soc. B, 69, 1020. 
CuyNnoweTH, A. G. 1958. Phys. Rev. 109, 1537. 
CuyNowETH, A. G. and McKay, K. G. 1956. Phys. Rev. 102, 369. 
Drews, R. E. 1957. M.Sc. Thesis, University of British Columbia, Vancouver, B.C. 
McKay, K. G. 1954. Phys. Rev. 94, 877. 
Newman, R. 1955. Phys. Rev. 100, 700. 
APPENDIX 
SPECTRAL DISTRIBUTION FOR RANDOM TRANSITIONS BETWEEN Two STATES 
The spectral distribution can be readily derived using the ‘‘g-r’’ theory for 
a univariate Markoff process (Burgess 1956). Consider N independent identical 
systems, each capable of existing in two states 0 and 1. Then in terms of the 
transition probabilities po; and pio the rate law for the number of systems a 
in state | is the simple linear equation: 
dn 3 
dt = (\ —n) Pon Pio 
where the first term corresponds to the generation (or forward) rate (g) while 
the latter corresponds to the reverse rate (7). 
The equilibrium value of » given by g(”) = r(n) is 
f= Nf=N in 
PoitPro 
where (cf. equation 1.2) we identify %/N with f, the probability of a system 
being in state 1. 
The relaxation time for perturbations from equilibrium is: 


& ds)" 1 

7 => ee oe eee 

dn dn n=n PoutPro 

and the variance of 7 is given by the Bernoulli result: 
var n = rg (n) = Nf (1). 


Since the fluctuations in 2 derive from a simple relaxation process the 


spectral density is 


For our case of a single discharge path we put NV = 1 and since the current 
change in the transition from the 0-state to the |-state is (¢;—/o), the spectral 
distribution of the current fluctuations is: 


roy — Sf) i= to)" 
(3.3) W (a) = Iter 


by virtue of the one-to-one correspondence between the state variable and 


the current 7. 
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CRYSTAL GROWTH AND THE THERMODYNAMICS 
OF IRREVERSIBLE PROCESSES! 


J. S. KirKALpy 


ABSTRACT 


The principle of minimum rate of entropy production is applied to steady- 
state transport processes in the neighborhood of an alloy crystal face growing into 
its melt. The procedure gives a satisfactory rationale of observed interface mor- 
phology. It is noted that segregation, which occurs in cellular or dendritic growth 
of alloys, is a direct manifestation of the system's attempt to minimize entropy 
production by conserving free energy. The general problems of growth of pure 
and impure single crystals from the melt and vapor are discussed. 


INTRODUCTION 


There are many dissipative systems in nature for which the external con- 
straints and the phenomenological transport equations do not uniquely 
determine the process path. The most familiar of such systems are those 
involving internal phase interfaces. These interfaces introduce further (internal) 
constraints into the system in the form of supersaturation relationships, con- 
tinuity conditions, and heat sources or sinks, but because of the flexibility 
of the interface morphology and the existence (in some cases) of surface 
transport mechanisms, the process path may still not be uniquely deter- 
mined. A familiar example of this is the steady-state unidirectional growth 
of an alloy single crystal into its melt. Evidently, for such a system we can 
satisfy the phenomenological equations and the constraints for a very wide 
range of growth conditions with the presumption of a flat interface. However, 
we know that for many conditions of growth within this range the interface 
actually attains a cellular or dendritic morphology (Tiller and Rutter 1956). 
The question then arises: What thermodynamic principle determines the 
choice of process path (if there is a choice) for a given complex irreversible 
process? 

In the past, various principles have been invoked for the case of crystal 
growth: choice of configuration for optimum heat transfer, the elimination 
of constitutional supercooling (in alloy crystal growth), or, somewhat equi- 
valently, passage through states of minimum free energy. While these perhaps 
contain elements of the truth in their specific applications, they have no 
fundamental foundation. Accordingly we have looked to the thermodynamics 
of irreversible processes for a suitable generalization. 


THE PRINCIPLE OF MINIMUM RATE OF ENTROPY PRODUCTION 
Prigogine (1955) and de Groot (1952a) have stated an important theorem 
applicable to single-phase dissipative systems in the steady state. We para- 
phrase this as follows: In a system which can be described by n independent 
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continuous thermodynamic forces (free energy or thermal gradients, chemical 
affinities), Xi, Xo,...X,, of which Xi, Xs, ...X, are fixed by the boundary 
conditions and a further number n—k of forces, Xz41, Xx+2,... X, are free to 
adjust themselves in conformity with the first k forces, the system approaches 
the steady state defined by the fluxes J; = 0 (¢ = k+1,k+2,...m) witha 
minimum in the rate of entropy production per unit volume, ¢. The J; are 
linearly related to the thermodynamic forces by 

(1) Siu 7 LS, (oF, 2 sit); 

j=1 
and are conjugated in such a way that 


n 


(2) lo= z Ja his 
i=1 

where 7 is the absolute temperature. The quantity To has the form of a 
sum of Joule heats and is accordingly known as the energy dissipation func- 
tion. 

In another formulation, Brown (1958) has shown that for a single-phase, 
multicomponent system in a steady conduction-diffusion—reaction state, the 
field equations can be derived from the variation principle 


. 
(3) | odV = 0 (minimum). 

“ 
This form is most useful for our present purposes. It is to be noted that both 
formulations of the principle of minimum entropy production assume the 
constancy of the phenomenological coefficients and require that the processes 
be near equilibrium for their validity. 

In what follows we plan to interpret these statements in a much wider 
sense than is implied by the mathematical formalism leading to them. In 
particular our variations will be made in the interface morphology without 
establishing the connection between these variations and the mathematical 
ones involved in the precise statements. Furthermore, we extend the principle 
to apply to heterogeneous as well as homogeneous systems. As such, these 
statements become postulates of qualitative behavior, and even though they 
have some theoretical basis and appreciable intuitive compulsion, they must 
stand in the final analysis on the agreement between the deductions from 
them and the observations. 

For non-stationary systems, we extend the steady-state principle by 
assuming the qualitative principle that any dissipative system starting from 
an initial condition with specified external constraints and proceeding towards 
a steady state (equilibrium is considered to be a steady state of order zero) 
proceeds through states of minimum integrated entropy production rate. We 
base this on the fact that one can think of any complex non-steady-state 
process as proceeding towards a continuous series of quasi-stationary states, 
provided the relaxation time for local change in entropy production rate (e.g., 
the time required for a crystal face to change its morphology a significant 
amount) is much less than the relaxation time for the whole process. 


: 
e 
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GROWTH OF ALLOY SINGLE CRYSTALS 


Consider the steady-state growth of a metal single crystal into a binary 
alloy melt. We will ignore for the moment any crystallographic or nucleation 
effects and consider the thermodynamics of a typical volume V (see Fig. 1) 
which migrates at the same constant rate v as the solid—liquid interface. It is 
not necessary for the following qualitative arguments to specify the third 
dimension of Fig. 1. However, in line with the observations, we can think 
of the cuspoidal liquid indentations as the elevation of an array of cell bound- 
aries. 





Fic. 1. Schematic advancing steady-state solid-liquid interface. 


The Variation Principle and Interface Morphology 
The variation principle (3) applied to a growing alloy single crystal is 


(4) 5 | odV ss s| | Be Ket Pe ta K) av+{ > JAuda) = 0 
vr Va 


where q refers to thermal and & to diffusion variables. A; is the chemical 
affinity and J; is the chemical reaction rate for reaction ¢t. Solidification of a 
binary alloy involves two reactions: components | and 2 from liquid (L) to 
solid (S). Because these are surface reactions we have transposed the corre- 
sponding volume integral into a surface one. 

The thermodynamic forces are (de Groot 19520) 


r 
(5) X, = -* 

(6) X, = -—TV (#), 
and 

(7) Ay = pir“. 


and if we ignore the Soret and Dutour effects (due to cross terms in the 
phenomenological transport matrix), the fluxes are 


(8) J, —AVI 

where A is the conventional thermal conductivity, 

(9) ha Les 
t=1 

and 


(10) J, = up,/M,, 
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where p; and M; are the masses per unit volume in the solid alloy and the 
atomic masses, respectively. If v is taken constant, the element of area in (4) 
must always be taken at right angles to the average growth direction. 

For dilute ideal binary solutions the integrands in (4) reduce approximately 
to (see Appendix J) 


— K po2, DRee 2 
(11) ae (VT) MC, (VC) 
and 
‘ I x v | Pi L s pe % 8 | 
) = — : Sei _ a 9 Lo 
(12) Ga 7 as JA; — (mi mi) 7, (ne Me ) 


where D is the diffusion coefficient, and R is the gas constant per mole. The 
concentration of solute, Cj, is measured in mole fractions. To apply the variation 
principle we set our boundary conditions so that the liquid metal has an 
initial composition of Cy and along the direction of growth fix two temperatures 
7; and 7» which include the melting range of the alloy. We then proceed to 
move these fixed temperature points at an equal and constant velocity v. 
Shortly, a steady state will be attained in the volume V, which for argument’s 
sake has a morphology as indicated in Fig. 1. Now we mentally vary the 
interface from the steady configuration to a virtual one involving a flat inter- 
face (dotted) and note thereby the change in the total rate of entropy production 
in V. This is a convenient variation since a tractable mathematical solution 
of the transport equations for a flat interface with the same growth rate is 
already known. 

Observationally, for a given alloy single crystal, cells appear only if the 
velocity exceeds a certain critical value (Walton et al. 1959). For conditions 
only slightly more severe than the critical ones, the cells are quite featureless, 
and crystallographic effects appear only in the extent and orientation of the 
cell boundaries. This suggests for these conditions that the variation of the 
reaction rate term in (4) is of secondary importance and that explicit contri- 
butions will only appear when the cells become sufficiently elongated that 
large areas of other low index interfaces are bared. We will accordingly neglect 
this term to the first order in a variation between a flat and slightly corrugated 
interface. 

The contributions of the conduction and diffusion terms to o can be evaluated 
in terms of known coefficients and in the neighborhood of the interface, near 
the critical conditions, they turn out to be of the same magnitude. For the 
particular example examined in Appendix I, the contribution per unit volume 
of liquid was three times greater in the diffusion field than in the thermal 
field. Comparing the qualitative liquid solute distribution of Fig. 2 with that 
for the corresponding flat interface we conclude that for the non-planar con- 
figuration the effective volume of integration in the liquid is relatively small, 
tending to a lowered total entropy production rate. On the other hand, the 
lateral gradients introduced tend to cancel this to an undetermined extent. 
The entropy production due to diffusion in the solid is negligible because of 
the low diffusion coefficients. 
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PROFILE 


Fic. 2. Representative steady-state liquid solute distribution near an advancing solid- 
liquid interface. 


Finally, referring to Fig. 2 and paying close attention to the projection 
profile and the requirement of substantial interfacial equilibrium, it seems 
plausible that in varying from a flat to a corrugated interface, the average 
latent heat source is moved closer to the sink at 7) as indicated in Fig. 1, 
with a corresponding decrease in the integrated conduction entropy pro- 
duction. It must, however, be noted that the decrease in entropy production 
will arise mainly through the replacement of some of the solid by liquid 
(Ks ~ 2K,) near the cell cusps. The change in effective volume of integration 
therefore appears to be much smaller than in the diffusion field so the con- 
tribution to the total entropy production change may be negligible. Some 
numerical estimates in favor of this conclusion are given at the end of 
Appendix ITI. 

Unfortunately, the qualitative nature of these statements precludes any 
unequivocal statement as to their validity. Accordingly, it is worth while 
to seek another less direct but more rigorous description of the entropy 
evolution of the test volume. 

Consider the total steady-state entropy balance for our volume V (de 


Groot 1952c), 


. . 3 . 
(13) \ Sav =~ | div (So= Qe wedi tv) dV+ odV =0, 

Jy ot Jy 1 Jy 
where s is the entropy per unit volume and V is the velocity of the center of 
mass relative to V. For an ideal solution in the approximation that the total 
latent heat production shows insignificant variations with concentration (or 
temperature) and that the temperatures 7) and 7% lie in regions where the 
components of J, in the growth direction vanish, the variation of the entropy 
production rate between a corrugated and flat interface becomes (see Appendix 


III) 
(14) A | odV A | divevdV =A | sv-nda = waas(T)), 
wy Le va 


where s(7\) is the average entropy per unit volume at temperature 7) and 
v = |v| is the speed of advance of the interface. Since As is experimentally 
always positive, this demonstrates the low relative entropy production of a 
non-planar configuration. Furthermore, since the variation on the right-hand 
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side of (14) is carried out isothermally and since the enthalpy per unit volume 
in the solid is relatively insensitive to local composition and defect structure 
changes, this becomes 


_vadf(Ts) 
~ 





(15) af odV= 
Vv 


where f is the average Gibb’s free energy per unit volume. Thus, the principle 
of minimum entropy production is manifested in cellular growth by the 
conservation of free energy in the solid solute segregation and (possibly) in 
the strain energy of defect structures. 

It is interesting to note that (15) implies a dominant role of the diffusive 
entropy production in the determination of interface morphology, for to the 
extent that 


‘div J 

(16) af £dV ~0, 
which was established in Appendix III, there is no mechanism available to 
create chemical free energy from thermal energy at the rate required by (15). 

While it is not clear as yet how to interpret the phenomenon in a precise 
mechanistic way, there appears to be an important connection between the 
principle of minimum entropy production and Hamilton’s principle. For a 
single particle system in a non-conservative force field this latter is (Whit- 
taker 1944) 





ell 3 
(17) | (67+ >> Q,6q,)dt = 0 (minimum), 


to r=] 


3 
where 7 is the total kinetic energy and }* Q,6q, is the work done by all the 
r=1 
forces in an arbitrary displacement of the co-ordinates g,. In general the 
Q, are time-dependent. In an alloy system in the steady state the forces 
associated with diffusion gradients and the change in chemical potential at 
a phase boundary can be assumed to be derived from an average potential 
P, which is a function of the co-ordinates only, so we can write approximately 


et . 3 
(18) 6 | | =P) “+ | - ora |at == () (minimum). 
J « vant 

This can be interpreted to mean that the particle follows a trajectory with 
the highest possible value of the potential P throughout consistent with the 
constraints. With suitable averaging it also seems plausible to suppose that 
an assembly of particles all starting out at some base potential at time ¢o 
should attain the highest possible potential or free energy at some later time 
t; in agreement with the inferences of the principle of minimum entropy 
production. 

The qualitative result of the analysis leading to equation (15) should be 
re-emphasized since it is perhaps a novel one. It states that this dissipative 
system passes through states of maximum free energy consistent with the 


ee 
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constraints. This points up the danger of applying the equilibrium principle 
of minimum free energy to any part of a dissipative system unless the process 
is an extremely slow one. 

While the tendency towards maximum free energy states has only been 
demonstrated for the steady-state alloy solidification process, it appears to 
have a qualitative applicability to all dissipative processes in nature. The 
possible extension of these ideas to other phenomena will be outlined in a 
later section. 


The Solute Distribution 

Whatever the steady-state interface morphology, there remains the necessity 
that (substantial) equilibrium be maintained at every point on the inter- 
face. This is only possible along the sides of the cells (with their continuously 
decreasing temperature) if a solute distribution can be provided which con- 
stitutionally matches the interface temperature distribution. 

With a starting alloy of normal constitution and of sufficiently high con- 
centration, this is automatically supplied by the pile-up of solute on the 
liquid side of the interface (see Fig. 2). If we presume, in line with the theoretical 
discussion above, that the cells extend themselves to the maximum consistent 
with the constraints so as to assure the optimum production of chemical and 
strain energy in the final solid, the diffusion field at the tip of the projections 
must satisfy the continuity relation that the actual thermal field (V7), and 
the constitutional thermal field (V7)-, derived from the solute distribution, 


be congruent, i.e., 
1 - =) 
(19) (V7), = (VWI)c = (=) (VO)1, 
dC/ 
where (d7‘/dC), is the liquidus slope obtained from the constitution diagram. 
In the terminology of the Chalmers’ school this means that constitutional 
supercooling must be eliminated at the cell tips. The limiting case in which 
(VC), for a flat interface satisfies (19) defines the condition for onset of 
cellular growth (Walton eft al. 1955). 
An estimate of the solute distribution for long cells and for (V7)s5 > (V7), 
can be obtained by imposing on the solution of the linear diffusion equation 


(20) Ci(x) ~ Aew""/¥+ Co, 


which originates at a cell cusp, the condition that at the latter point 


(21) (V7) & (VT)s = (42) wor. 
The calculation of the value of x in (20) which satisfies (19) gives an estimate 
of the cell length which is of the order of that observed (see Appendix I1). 

It remains to make an estimate of the cell spacing, d. The controlling 
constraint is the requirement that there be appreciable lateral mass transfer 
to assure the solute balance in the system. The condition for this is the 
requirement that an atom be able to move from the tip of the cell to an 
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intermediate point in a time of the same order that it takes the interface 
to move a similar distance. Hence, 






(22) euake 
2 







This leads to values of the spacing within a factor of about two of those 
observed (cf. Tiller and Rutter 1956, Table V) and shows the right qualitative 
variation with the variable v. This estimate does not, however, account for 
the variation with (V7); observed by Rutter and Chalmers (1953) and Tiller 
and Rutter (1956). 








DENDRITE FORMATION 


Dendritic Growth into Alloy Melts 

In the discussion thus far we have ignored any crystallographic effects at 
the interface. This is a reasonable approximation for slight changes in the 
profile of the interface since the supersaturation (or supercooling) required 
to maintain steady-state growth is relatively constant over the whole inter- 
face. The observations bear this out since for conditions slightly more severe 
than the critical ones defined by equation 19, the corrugations or cells are 
quite smooth and featureless, and the crystallography appears only in the 
orientation and length of the corrugation or cell boundaries. We naturally 
expect crystallographic effects to appear as soon as the cusps are deep enough 
to bare an appreciable area of different low index plane to the liquid. This is 
because the supersaturation conditions are different for different crystal faces. 
According to the principle of minimum entropy production, for a given 
velocity preferred growth will occur on those faces which have the lowest 
chemical entropy production as indicated by relation (12). This can rationalize 











































the dependence of cell boundary orientation and the tendency towards side- 
branching on equivalent preferred planes when the cell cusps become deep 
enough. As the growth conditions become more severe, the field components 
at right angles to a primary preferred growth direction of a cubic crystal 
approach a condition wherein (19) can be satisfied in that direction and there 
will be then a tendency towards side-branching from the newly bared faces 
of low dissipation. Branching, however, is incompatible with local steady 
state so an unstable condition arises. 

An overeall steady state of minimum entropy production can still be main- 
tained if an available coupled mechanism of decreasing local entropy pro- 
duction allows the system to enter an oscillatory mode. Figure 3, showing a 


== 5 br | 


Fic. 3. Representative cyclic configuration near an advancing dendrite tip. 





series of qualitative pictures equally spaced in time of a dendrite tip, indicates 
that this is in fact what happens. The decreasing local entropy production 
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which sustains the oscillation presumably occurs in the build-up of solute 
and latent heat which occurs as the branches and dendrite cores make their 
alternating forward spurts. It is further necessary to presume that a region 
of constitutional supercooling oscillates between regions A and B on Fig. 3 to 
anticipate each spurt so that the sudden latent heat build-up can be accom- 
modated. 

Since the local entropy production is now a function of time it is obvious 
that (3) can no longer describe the details of the process. However, the prin- 
ciple of minimum entropy production might be preserved by forming a double 
variation integral, 

(23) &ffo dV dt = 0 (minimum). 


If ¢ could then be expressed as a suitable quadratic function of independent 
field variables and their time rates of change then the equivalent Lagrangian 
field equations would give the required oscillatory solutions. The mathematical 
complexities, of course, preclude any detailed procedure of this type, but in 
principle it would seem to represent the situation adequately. 

The frequency of oscillation is determined by the relaxation time of the 
diffusion field since this is the slowest process in the system. This is the time 
to transmit an interaction across the dendrite spacing d, so the frequency 


D 
© yp ~~ 5 
(24) oes 


and by (22) the branch spacing 


(25) a= : ~ (S)-4 ~d 


as is usually observed. 

The application of the total entropy balance (13) to a variation between 
a dendritic and a flat interface is complicated by the fact that (15) may no 
longer be a very good approximation because of the concentration dependence 
of the latent heat production. However, there is little doubt that (15) still 
qualitatively describes the process and that we can interpret the stored free 
energy in the final solid as an indication of the system’s tendency towards 
minimum entropy production. 

In the previous section we proposed that no constitutional supercooling 
occurs at the tips of cells on a solidifying interface and in this one we have 
found it necessary to propose a region of supercooling oscillating between 
the dendrite tips and the adjacent interdendritic spaces in steady-state growth. 
Plaskett and Winegard (1958) have demonstrated the existence of consti- 
tutional supercooling ahead of a slowly growing, non-steady-state dendritic 
interface by observing the development of equiaxed grains in that region. It 
may be that the nucleating regions are the pockets of supercooling required 
by our analysis of dendritic growth, or it may be that a more extensive, 
continuous region of supercooling develops ahead of a non-steady interface. In 
either case, the extent of these regions cannot be appreciably greater than 
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the diffusion length D/v which was about 10 microns for their conditions. 
Accordingly, nucleation of the new grains must have occurred very close 
to the growing interface. 


Dendritic Growth into Pure Melts 

Weinberg and Chalmers (1951) have demonstrated that dendritic growth 
only occurs in pure melts if there is an appreciable region of liquid super- 
cooling ahead of the interface. A typical observation in commercial pure 
ingots is a top surface skin of dendrites which presumably derives from the 
supercooling caused by radiation and convection from this uninsulated sur- 
face. This process is almost inherently a non-steady-state one because of the 
experimental difficulties in producing a moving zone of supercooling. 

As before, there is a choice of solidification paths compatible with the 
constraints and phenomenological equations (thermal conduction and chemical 
reaction only). In this case the only important constraint is that the latent 
heat production and the coupled thermal field be compatible with the super- 
cooled surroundings and boundaries. The interface morphology is accordingly 
very flexible under these conditions. 

If we examine the dominant thermal conduction term in terms of an equi- 
valent quasi-stationary process described by the variation principle (3) we 
form the qualitative conclusion that dendritic growth reduces the entropy pro- 
duction as compared to a flat interface by distributing the latent heat source 
as intimately as possible within the sink, shortening the transport paths 
and lessening the average thermal gradients. The corresponding increase in 
the available energy state is provided by the increase of surface energy. 


GROWTH OF CRYSTALS FROM THE VAPOR 

We generally note in sublimation processes a rather greater sensitivity to 
crystallography than in solid—liquid reactions. This is due in part to the 
great sensitivity of growth rates to the supersaturation and to the existence 
of different growth mechanisms and saturation conditions on different crystal 
faces. Nonetheless, the same thermodynamic principles should apply in 
solid-gas reactions as above. Application of these principles may often be 
simplified by assuming (with some justification) a quasi-stationary state in 
the transport fields, since in this case (3) applies. 

While we reserve a full discussion of this complicated problem to a tater 
contribution, it is interesting to note that the existence of efficient surface 
transport mechanisms may result in the freeing of certain of the transport 
forces. This means, as a result of Prigogine’s theorem, that at the steady 
state, the conjugate fluxes vanish. If we have two components of surface 
flow on a metal whisker or crystal with low entropy production rate, then 
two components of the external diffusion and (or) convection (depending on 
whether there is an impure or pure vapor environment) fields may be freed. 
These will evidently be the two angular components in spherical co-ordinates 
so the diffusion or convection fields can be (approximately) spherically sym- 
metric, irrespective of the interface morphology. The rationalization of the 
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high observed growth rates of metal whiskers (Sears 1955; Coleman and 
Sears 1957) and snow crystal dendrite branches requires just this sort of 
mass transport configuration. 

It is interesting to consider in this respect the experiment of Bunn (1949). 
In two-dimensional crystallization from a solution he observed a circular 
solute distribution about an approximately square crystal. In the discussions, 
these results were discounted by F. C. Frank because a surface diffusion 
mechanism was unable to rationalize them. In the light of the present work 
it would seem attractive to postulate an alternative surface transport mecha- 
nism—perhaps some sort of monolayer viscous transport—which would explain 


the results. 


FURTHER APPLICATIONS OF THE PRINCIPLE OF MINIMUM ENTROPY 
PRODUCTION 

The principle of minimum entropy production appears to have qualitative 
applicability to all dissipative processes in nature. This is particularly evident 
if we regard as an approximate corollary the statement that ‘a dissipative 
system passes through states of maximum available energy compatible with 
the constraints’’. In the case of a system proceeding towards a steady state 
this is exemplified by the build-up of potential energy and the maintenance 
of high dynamic efficiency in the steady state. Good examples of this are the 
thermal diffusion cell, the steady-state alloy growth considered here, the 
vortex in the bath tub drain, or the history of the adult biological organism. 

In the many common processes which are proceeding towards equilibrium 
we invariably note the tendency to pass through states of high available 
energy. Ostwald’s rule for chemical reactions (Darken and Gurry 1953a), which 
notes that many reactions pass preferentially through metastable states before 
reaching the ground state, appears to be closely related to this principle. The 
many exceptions to the rule appear to express nothing more than the rate 
dependence of the entropy production, i.e., fast reactions will tend to go 
through higher energy states than slow ones and very slow ones will tend to 
go directly to the ground state. Many examples of such a rate dependence 
are familiar in metallurgical transformations. 

Finally, one can qualitatively extend the application of the principle to 
the evolution of both the universe and the species by noting the local tendency 
towards organization as universal dissipation proceeds. 


SUMMARY 


The principle of minimum entropy production (or energy dissipation) has 
been applied to the dissipative processes involved in metal crystal growth 
and this procedure leads to a satisfactory rationale of surface morphology. It 
is proposed that in solidification from an alloy melt, diffusion supplies the 
dominant source of entropy production, and the system tends to minimize 
this by arranging the morphology in such a wavy as to conserve available 
energy through solute segregation. In a supercooled pure melt, the morphology 
is arranged for low entropy production by distributing the latent heat source 
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uniformly within the sink, i.e., by a dendritic formation which exhibits 
corresponding states of high surface energy. 
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APPENDIX I 





To simplify the entropy production expression, 
: | 
(26) fo = git+o2 i a, J,-X,+ > JX, ’ 


we obtain directly from (5) and (8) 


(27) a, = . (V7)°. 

7 
For a binary alloy 
(28) Ji = Ly Xit+LyX2 = Li(Xi— Xe) 
and 
(29) Jo = Ly XitLnX, = —Ly(Xi— X2) 
since (de Groot 1952d) 
(30) Ly = -—Ly = —La Lo. 
Thus 


(31) 





03 = Ly(X1— X2)?/T. 





i 
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Neglecting the 7 dependence of equation (7) and using Fick’s relation 


(32) J; = —DVG, 

we obtain for a dilute, ideal solution, 

(33) a= DE WC)" 
1 


The units of the concentration C, are mole fractions. Since R is the gas constant 
g 

per mole, we convert a2 to unit volume and finally obtain equation 11. To com- 

pare a; and oz at the critical condition defined by equation 19 we obtain 


_ | K, DR (26)" | 1)? 
(34) =| aT /, (VT). 


Using the critical values for the 0.03% tin in lead alloy given in Tiller and 
Rutter’s Table V and the values for lead, pp = 10.5 g/cm’, Mz = 207.2 g/mole, 
and Ky, = 0.04 cal/sec cm °C, we find that on the liquid side of the critical flat 
interface, the diffusion term per unit volume is about three times the con- 
duction one. 


APPENDIX II 
To make an estimate of cell or dendrite lengths we assume the liquid solute 
distribution (20) 
C~ CotdA e7*/™* 


originating at point B on Fig. 2 and impose the condition 


(35) vC(0) = w1)s(26) 
dl L 


to determine 1. We then solve for the value x = x’, which satisfies 


| es 
(36) VC(x’) = (V7 ($6) . 


This procedure should lead to reasonable estimates provided (V7) 5 > (V7) . 
Taking as typical values those given by Tiller and Rutter (1956, pp. 105 

and 116), we have: 

Co = 0.01% tin in lead, v = 0.0167 cm/sec, (VT), = 5.4°C/cm, D=7X10~ 

cm?*/sec, and from handbooks Z = 5.89 cal/g, Ks(Pb) = 0.08 cal/sec cm °C, 

K,(Pb) = 0.04 cal/sec cm °C, p(Pb) = 10.5 g/cm’, p(Sn) = 7.3 g/cm’. 

We obtain (V7) from the continuity relation for a flat interface, 


(37) K3(V7)s = Ky (VT). +upl. 
This gives (V7) 5 = 15.6° C/cm and satisfies our condition that (V7) 35> (V7) . 
Carrying out the calculation we obtain A = 0.022°% and x’ = 45 microns. 


This is in order of magnitude agreement with the cell profiles shown in Tiller 
and Rutter’s Fig. 27. 
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APPENDIX III 
For the simplification to equation (14) of the expression for the difference 

in total entropy balance between a flat and corrugated interface, we note first 
that for steady state and for reference temperature points 7; and T> well 
removed from the interface 

Pa div Re ony 
(38) af “dV = af divi mde gy = 0, 

V ot V t 


Accordingly we can write 


(39) af oav = af Mw Xaveaf Seay, 
Vv V V 


and it remains to indicate under what conditions the last term on the right- 
hand side is negligible compared with the first. 

We can write the heat balance at constant pressure 
(40) divJ, = —C +0 

ot 

where Q represents the heat sources and sinks in the system. This will include 
the latent heat source, chemical free energy converted to heat in the reaction 
and diffusion fields, and a sink corresponding to conversion of thermal energy 
into chemical energy. This latter could arise from thermal diffusion effects, 
but because of the relatively small magnitude of such effects, this term can 
be considered as negligible. We can now invoke the steady-state condition 
that 


(41) 


divide by 7 and integrate to give 


r div Ben tee 7 { Copp al i, qd V } pul, 
af T adV = 4 » T ax d| +4 . T +4 tT da, 


where g represents the non-latent chemical heat sources and L is the latent 


(42 


o 


heat. 
We now consider each of the terms on the right-hand side separately. A 

non-zero value of the first one would arise mainly from the exchange of some 

liquid and solid in the solidification zone giving rise to a change AC,. We 

can, therefore, approximate this term for the system considered in Appendix 

I] by 

_ AGon aT 


; cal deg sCcCCc, 
1 Ox / b 


(43 ) Ay a 


/ 


fav ~ £1107 - 


In this estimate we have used AC, = 0.001 cal/g (Kelley 1935), pz: = 10.5 
g/cm, v 0.0167 cm/sec, and d7°/dx ~ 10° C/em. The effective volume of 
integration is probably zero to the first order, but as a conservative estimate 
we have represented it by one tenth of a solidification zone of 45 microns 


width. 
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The second term represents the conversion of free energy in the diffusion 
field into enthalpy. For ideal solutions this is zero since the enthalpy of 
mixing is zero. More generally it may contribute another (usually fractional) 
term to the right-hand side of (14). For a regular solution (Darken and 
Gurry 19530) the difference in the rate of heat evolution in the diffusion 
field per degree absolute is estimated to be 


(44) Ay ~ a Ez Whe cq ~ +1.5X — cal /deg sec 
where a generally lies between +2 (we have chosen a = +1 here) and AC; is 
the concentration difference between the interface and the alloy concentration 
(0.00014 in this calculation). The bracketed quantity is the rate of heat 
evolution for a flat interface and the difference between this and that for a 
non-planar interface, A: (which is zero to the first order), is taken as one 
tenth of the whole. 

Finally, the latent heat Z depends on the interface concentration (or 
temperature). We therefore make the estimate 





da ~ +1X10° = cal deg sec. 


pe AG ak 
7 


(45) A; T did, 
In this we have taken Z = 5.89 (14(C,) cal/g and AC, = 1/404. (See Appendix 
II for calculation of A.) The mean change in one or the other direction is 
probably about }.1, but to the first order the positive changes will cancel 
the negative ones so we have reduced this even more. 

For comparison we estimate the right-hand side of (15) as 


(46) —vAf =~ 1X10 = cal/deg sec. 

re 1 
In evaluating this we have assumed an ideal solution and a segregation of 
0.01% existing in one twentieth of the solid volume. 

As an alternative we can estimate the maximum change in entropy production 
due to change in the diffusion fields by means of equation 12. This change 
should be of the order of one-half the entropy production existing in a volume 
near a flat interface which is equal to the actual volume of the solidification 
zone. For an ideal solution this gives 


4 


r 


| DRps 
2 MC, S 


” 


(47) 


. 
Cy)" | dV ~04AX 10 cal /deg sec 
« 
in rough agreement with (46). 
lor comparison, the corresponding conduction term is 


cal deg sec 


(48) AK wr) fav ~0.7x10°° 


i l 
Hlere, the main change would appear to be through the possible exchange of 
some liquid for solid with a corresponding change in the thermal conductivity 
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We have taken AK = 0.04 cal/seccm, V7’ = 15.6° C/cm and the effective 
volume as one tenth of the solidification zone. The term appears to be small 


compared with the diffusion term (47). 
The general result of these calculations is to indicate that the quantity 


is made up of a sum of terms of variable sign which to the first order vanishes. 
Our estimates of the second-order variations suggest that in ideal solutions 
and most others as well, these terms can be neglected in comparison to (46) 
and that equations (14) or (15) give an accurate description of the solidifi- 
cation process. Even if there exist cases where the latter terms are not negligible, 
the inference remains valid that solid solute segregation attending non-planar 
growth is a measure of the system’s attainment of minimum entropy pro- 
duction, since the segregation term in (14), like the entropy production change, 


is always positive. 





SPINS OF EXCITED STATES OF Os!**! 


W. J. Kinc ano M. W. Jouns 


ABSTRACT 


Seven excited states of Os!** populated through the beta decay of Re'’* have 
been studied by gamma-gamma angular correlation experiments. The data for 
the various cascades fit the following theoretical functions: (all gamma-ray 
energies are in kev and the description following each cascade refers to the 
first transition) 478-155 (2-2-0, 99.7% E2, 6 +ve), 931-155 (0-2-0), 672-633 
(2-2-0), 97% M1, 6 —ve or 3-2-0, 94% M1, 6 —ve or conceivably 4-2-0), 
828-633 (2-2-0, 98% M1, 6 —ve or conceivably 3-2-0, 80% M1, 6 —ve), 
1132-633 (0-2-0), and 1308-633 (2-2-0, 98% M1, 5 —ve). Some attenuation 
occurred in the correlations involving the 155-kev gamma ray. The K-conversion 
coefficients for the 478- and 633-kev transitions have been measured as 0.023 
+0.003 and 0.010+0.002 respectively, confirming the £2 character for these 
transitions obtained by angular correlation. These results, coupled with data 
concerning log ft values and gamma-ray intensities, lead to the following spin 
and parity assignments for levels in Os!88: 155 (2+), 633 (2+), 1086 (0+), 
1306 (2+ or 3+), 1461 (2+), 1765 (0+), 1941 (2+), and 1958 (1+ or 2+). 

The 631-137 correlation data in Os'8* was found to fit an attenuated 2-2-0 
function with the first transition 99% E2 and 6 —ve. This result supports the 
established decay scheme for this nucleus. 


INTRODUCTION 


Seventeen-hour Re'*’ has been the subject of many investigations but no 
well-established decay scheme was proposed until 1956. In that year, workers 
in this laboratory (Johns et al. 1956) and at the Bartol Research Foundation 
(Potnis 1956) independently arrived at essentially the decay scheme of Fig. 1. 
An investigation of the shapes of the beta-spectra feeding the 0- and 0.155- 
Mev levels of Os'*§ (Nielsen and Nielsen 1958) revealed that both spectra 
show departures from the allowed shape which establish them as non-unique 
first-forbidden transitions. Since the spins of these two levels are known to 
be 0+ and 2+, one is forced to assign spin 1 and negative parity to the 
ground state of Re'’’. From the gamma-ray intensity balance, Johns et ai. 
(1956) found that all of the excited states of Os'’’ shown in Fig. 1 are fed 
by first-forbidden beta transitions. Hence one is led to expect that the spins 
of all levels of Os!** excited by the beta decay of Re'’® must be 0, 1, 2, or 3 
with positive parity. Levels with higher spins can only be excited indirectly 
through cascading gamma rays. 

Both angular correlation (Potnis ef al, 1956) and internal conversion 
measurements (Nielsen and Nielsen 1958; Johns et a/. 1956) led to a 2+ 
assignment for the 0.633-Mev level, but nothing was known concerning the 


spins of the weakly populated higher levels. This paper describes a series of 


angular correlation experiments whose results lead to an assignment of spins 
to five of these states. 

The beta decay of Re!’® populates levels in Os!'** at 0.187 and 0.768 Mev 
above the ground state (Metzger and Hill 1951) and leads to the 631-137 


‘Manuscript received January 6, 1959. ; 
Contribution from Hamilton College, McMaster University, Hamilton, Ontario. 
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1958 (0.012%) 
1802.7 (0.027%) 
1783.0 (0.015%) 
1307.9 (0.02%) 


672.6 (0.09%) 


453.6 (0.04%) 


2+ 633.1 


(9.0% e-) 
(9.0% ¥) 


° 
o 
3 


2+ 155.0 


0.0 


Fic. 1. Decay scheme of Re!*’, 


cascade* in 0.025% of the disintegrations (Porter et a/. 1955; Johns et al. 
1956). The angular correlation of this cascade was studied by Hurley and 
Jastram (1954) but their results are available only in abstract form. In view 
of this fact, it was felt worth while to re-examine this cascade as part of the 
present investigation. It has been clearly shown by many workers that the 
0.137-Mev transition is pure E2 and that hence the cascade of interest is of 
the type J—2-0, with all the available evidence suggesting that J = 2. 
Recent investigations of the decay of Ir'*® and Ir'*’ by Diamond and 
Hollander (1958) have identified a number of the higher spin states belonging 
to the ground-state rotational sequence 0, 2, 4, 6... in Os'8® and Os!8*, A 
similar sequence has been established in Os! (Scharff-Goldhaber et al. 1958) 
and the presence of two low-lying spin 2 states has been confirmed in all 
these nuclei and in Os!® by the Coulomb excitation experiments of Barloutaud 
et al. (1958). The last-named workers also reinterpret earlier data available 
concerning Os!” (Baggerley et al. 1955) to make it fit the pattern established 
from the other osmium nuclei. None of the states with spin greater than two 
have been identified in the rhenium decay processes, but all of this data will 


*In the interests of brevity, gamma-ray cascades will be simply labelled by the energies of 
the two gamma rays in kev without explicit mention of units. 
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be used in the discussion of the results of the present experiments and is 
summarized in Fig. 7. 
APPARATUS AND METHOD 

The gamma-ray angular correlation spectrometer consisted of two scin- 
tillation detectors mounted in the same horizontal plane with the source 
located 7 cm from the front face of the detectors. Each detector consisted 
of a 1} in. X 1 in. Nal(TI) crystal coupled to an R.C.A. 5819 photomultiplier. 
The angle between the fixed and movable detectors could be altered by rotating 
the latter in a horizontal plane about a vertical axis through the source. 
Coincidences between photons detected in the two detectors were observed 
with a Beil-Graham-—Petch (1952) coincidence circuit modified in the following 
manner. The pulse from each photomultiplier was inverted, clipped by means 
of a 200-ohm cable to 4X 10~-§ seconds, and then amplified with two Hewlett- 
Packard wide-band amplifiers (460A and 460B) in cascade. The amplified 
pulses were fed to a 404A pentode limiter whose output led to the fast coinci- 
dence circuit. The resolving time of this unit was adjustable from 2x10-° 
to 2X 10-8 seconds, with a value of 8X 10~* seconds being commonly employed. 
The output of the fast coincidence unit recorded coincidences between all 
pairs of gamma rays detected in the scintillators. Energy selection was 
performed in both side channels by means of single-channel pulse height 
analyzers (Tracerlab RL1-6). The outputs of these channels and the fast 
coincidence unit were fed to a triple coincidence unit of resolving time 5X 10-7 
seconds whose output yielded coincidences between pulses of selected heights. 
Under the conditions of the experiment, the maximum usable source strength 
was determined by the number of pulses reaching the side-channel analyzers 
rather than by the accidental coincidence rate. Wherever possible, unwanted 
low-energy pulses were prevented from reaching these analyzers by removing 
the associated X rays or gamma rays with 1 g/cm? lead or 0.6 g/cm? tin 
absorbers placed over the face of the detector. 

Sources were prepared from spectroscopically pure natural rhenium or 
rhenium enriched to 98.2% in the heavier isotope. The solid sources con- 
sisted of about 1 mg of polycrystalline metallic powder while the liquid 
sources consisted of rhenium dissolved in hydrogen peroxide. The liquid and 
solid sources gave essentially the same anisotropy for all the Re'** cascades 
but there was a noticeable difference between the two for the Re!*® cascade. 
Most of the Re'*’ data was obtained from metallic sources, irradiated in 
sealed quartz vials about 10 mm long and 1 mm in diameter, and used in an 
unopened condition. The primary beta rays were absorbed by inserting the 
source vial in a hole drilled along the axis of a carbon cylinder mounted 
concentrically with the axis of rotation of the spectrometer. The position of 
this cylinder was adjusted to locate the small source accurately at the center 
of the circle defined by the intersection of the horizontal axes of the crystals. 
The correct vertical position was first determined by visual inspection and 
then the lateral position adjusted until the counting rate in the movable 
detector was constant to within 1% as the detector was rotated from the 
90° to the 270° position. 
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For each cascade studied, an experimental resolution curve was obtained 
by adjusting a variable delay line in the fast coincidence unit. The gamma-— 
gamma correlation data was then obtained with the circuit set at the delay 
corresponding to the peak of this curve. The chance rate was measured either 
by inserting a sufficient delay in one channel leading to the fast coincidence 
circuit or by inserting a large delay in one of the side channels. Under the 
conditions of these experiments the chance rate was never more than 5% of 
the true coincidence rate. 

Before studying a cascade, each side-channel spectrum was determined and 
the “‘windows”’ set on the desired portions of the photopeaks of the two 
gamma rays involved using the same geometry as that in which the correlation 
data was to be recorded. With these window settings, the number of coinci- 
dence counts and the number of counts in the movable counter obtained in 
15 minutes were recorded for a sequence of angles ranging from 90° to 180° 
and back in 22.5-degree steps. Following this sequence, the side channels 
were recalibrated to check for drifts, and, in the rare cases where appreciable 
drifts had occurred, the data was rejected. For most of the cascades, the 
sequence was repeated until several thousand coincidences had been recorded 
at each of the five angles. For each angular position, the ratio R(@) = C,(0)/N(@) 
was formed. In this expression, N(@) is the cumulated number of counts 
recorded in the movable detector for a given angular position while C,(@) is 
the cumulated number of observed coincidences at the same position corrected 
for the small number of chance coincidences present. The function R(@) 


embodies the experimental information concerning the angular correlation 
of the cascade under study. It can be seen that this method of treating the 
data not only corrects for decay but also introduces a first-order correction 
for errors in centering of the source on the axis and for small window drifts. 
Using the least-squares method outlined by Rose (1953), the function R(@) 
was used to determine the A‘ coefficients in the expansion 


R(6)/C = 1+ A4P,(cos 6) + AP, (cos 6) 


where C is a constant of proportionality of no further importance. 
These At coefficients were corrected for the finite angular resolution of 
1e detectors in the manner described by Rose (1953), using the functions 
explicitly evaluated by him for a source-to-crystal-face distance of 7 em and 
3 in. X 1 in. detectors.* The A,'s so obtained represent the true coefficients 


"It was pointed out by the reviewer that the correction factors for the finite size of the 
detectors give ay: Kose are not strictly applicable to an experiment in which only the photo 
peak pulses are counted. To answer this criticism, the authors have determined the correction 
factors experimentally, using the method des« ribed by Lawson and Frauenfelder (1953) but 
setting the pulse height analyzer to accept only the photopeaks of the gamma rays from the 
collimator. ‘The solid augle correction factors were determined for radiations of energy 0.129 
Mey (Os'*), 0.411 Mev (Au), and 1.43 Mey (Co®), As might be expected, the experimental 
factors obtained were slightly le ss than those computed by Rose; the differences ranging from 
zero at 0.129 Mev to a maximum of 0.2% for the J,/Jo ratio and 0.8% for the J4/Jo ratio at 
an energy of 1.4 Mev. These measurements indicate that the A, values in Table I should be 
reduced by U4Y% while the A, values should be reduced by amounts ranging from 0.8 to 
1.2%. Since other uncertainties in the measured coefficients are of the order of 10%, the 
ntroduced by the use of Kose’s coefficients are trivial and do not alter any 


of the conclusions reached 


systematic errors 
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for the selected pair of gamma rays combined with whatever contribution 
from other gamma rays was included in the windows chosen. A discussion of 
these effects is included in the section dealing with experimental results. 
EXPERIMENTAL RESULTS 
Figure 2 presents the gamma-ray spectrum of Re'®® as observed in either 
scintillation detector. This fails to reveal the detailed structure of the spectrum 
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Fic, 2. Scintillation counter spectrum of the gamma rays of Os!**. Lower curve-—-gamma 


rays in coincidence with the 0.155-Mev radiation. 

obtained with magnetic spectrometers (Johns ef a/. 1956) but is sufficiently 
well defined to enable one to pick out any cascade of Fig. | without too much 
interference from the others. The windows used in the angular correlation 
studies are indicated above the photopeaks. The weak 0.672-Mev peak is 
masked by the much stronger 0.633-Mev peak in the single-channel spec- 
trum; its intensity, deduced from the coincidence spectrum, agrees with the 
value obtained with the magnetic spectrometer. 

The experimental data for the cascades of Os!’ are presented in Fig. 3 
and in Table I. Each figure shows the experimental points, the least-squares 
fit to these points, and the least-squares fit to the data corrected tor the 
finite solid angle of the detectors. The differences between the two curves are 
not large and small errors in these corrections would not significantly alter 


the conclusions which are presented in the following pages. Table | gives the 
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Fic. 3. Experimental angular correlation functions for the gamma rays of Os!**, ‘The 
broken line represents the least-squares fit to the experimental points. The solid line is the 
correlation function corrected for the finite solid angle of the detectors. 


solid-angle-corrected coefficients for each cascade together with pertinent 
information concerning the experimental conditions. It also presents the 
coefficients for the ‘‘pure’’ cascades obtained by subtracting the effect of 
interfering cascades. These background-corrected cascades are shown in Fig. 
5. The difficulty in estimating the magnitude of such backgrounds of course 
leads to rather large uncertainties in the “‘pure’’ coefficients; in some cases, 
these large uncertainties made it impossible to reach unique spin assignments 
The standard deviations on the experimental coefficients quoted in Table I 
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TABLE I 


Angular correlation functions of Os!88 and Os!86 


Energies of Windows Numberof Average 
cascading used counts rate Intensity of 
gamma rays, Type of (see at each per interfering 
kev source Fig. 2) position minute cascades Ae 


478-155 Os!88 = Liquid Aand B 5500 20 10% due to —0.004+0.013 +0.199+0 
Solid 4500 20 931-155 etc. —0.030+0.029 +0.213+0 
Total 10000 20 —0.016+0.010 +0.205+0 
Pure 478-155 —0.04840.010 +0.162+0 
931-155 Os!88 = Solid A and F 1500 4 20% due to +0.271+40.032 588 +0 

(1000+ )-155 
Pure 931-155 +0.31340.060 +0.640+0 
(1000+ )-155 Solid A and L 1300 +0.10140.060 +0.379+0 

Os188 
28-633 Os!88 Liquid Cand E 4500 g 12.5% dueto +0.193+40.008 +0.133+0 
1308-633 
Solid 4500 12.5% dueto +0.216+40.019 +0.146+0 
1132-633 

Total 9000 4 +0.205+0.010 +0.139+0 
Pure 828-633 +0.222+0.02 +0.026+0 
672-633 Os'!88 Solid D and D 5000 17% 828-633 +0.131+40.023 +0.124+0 

8% 1132-633 

8% 1308-633 
Pure 672-633 +0.103+0.05 +0.06 +0 
1308-633 Os!88 Solid C and H Pure 1308-633 +0.180+0.005 +0.017+0 
1132-633 Os!88 Solid Cand G Q 2 50% 1308-633 +0.153+40.037 +0.479+0 
Pure 1132-633 +0.13 +0.08 0.94 +0 
631-137 Os!8¢ = Solid 137 d None —0.026+0.003 +0.061+0 

Photopeak 
and all 
pulses above 
200 kev 3000 2 None —0.056+0.003 +0.123+0.004 


Liquid 


were computed using equation (30) of Rose’s (1953) article. In using this 
equation, it is assumed that the standard deviations in the individual points 
determine only the relative weight of the points in the least-squares analysis, 
but that there is no uncertainty in the position of the points themselves. This 
assumption seems to be reasonable when the number of coincidences at each 
point is large. In such cases, the standard deviations of the coefficients calcu- 
lated in this manner give a true indication of the sensitivity of the function 
to be fitted to variations in the coefficients. However, when the number of 
coincidences recorded at each position is small, the uncertainties in the 
positions of the individual points may well be the dominant factors in deter- 
mining the errors in the coefficients. The appropriate expression for calculating 
the standard deviations in the coefficients is then equation (20) of Rose's 
paper. In the analysis of the results, both equations were used and found to 
give essentially the same standard deviations except for the 672-633 and 
1308-633 cascades; for these, the limits of error deduced from equation (20) 
were several times larger than those deduced from equation (30). A casual 
inspection of Fig. 3 shows that the least-squares fit to the experimental 
points for these two cases is considerably better than one would expect on 
the basis of the statistical accuracy of the individual points. This suggests 
that larger limits of error are more realistic than those quoted in Table | for 


these two cascades. 
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The 478-155 and the 828-633 cascades were studied with both solid and 
liquid sources. Since both types of source gave essentially the same correlation 
function, the two sets of data were averaged in interpretation as shown in 
Table I. 

Although Potnis et al. (1956) had demonstrated from angular correlation 
studies that the 0.633-Mev level was 2+ in character and Johns et al. (1956) 
had measured a, for the 0.478-Mev radiation, the discrepancies in the mixing 
ratio for the 0.478-Mev radiation between the two sets of results made it 
seem worth while to remeasure the internal conversion coefficients of the 
0.478- and 0.633-Mev transitions. These were obtained by using a source of 
rhenium enriched in the heavier isotope and a 50-cm beta-ray spectrometer. 
The values of 0.023+0.003 and 0.010+0.002 obtained for the K-conversion 
coefficients of the 0.478-Mev and 0.633-Mev transitions respectively are in 
good agreement with the theoretical values for pure E2 transitions and the 
recent experimental values of Nielsen and Nielsen (1958). 

The 631-137 cascade of Os!*® was examined with both solid and liquid 
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sources. In order to prevent overloading of the side-channel analyzers by X 
rays and low-energy gamma rays, absorbers of 0.9 g/cm? lead and 0.6 g/cm? 
tin were used over the high-energy and low-energy detectors respectively. 
About 3000 coincidences were recorded at each angle for the liquid source 
and about 8000 coincidences per position for the solid source at an average 
coincidence rate of 2 per minute. The results of these experiments, presented 
in Fig. 4, (a) and (0), show clearly the existence of attenuation due to the 
reorientation of the nuclei in the relatively long-lived intermediate state 
(8X10-!9 seconds; McGowan 1951). 


DISCUSSION OF EXPERIMENTAL RESULTS 

From Fig. 1, it is evident that the only cascades of Os'*’ present in appre- 
ciable intensity involve the ground state and one of the 2+ states at 0.155 
and 0.633 Mev. Hence the second component of all the cascades studied is 
pure E2 and the experimental functions may be compared with theoretical 
curves for cascades of the form J-2-0 where J is restricted to the integral 
values from 0 to 3 by the fact that the ground state of Re!S$ is 1— in character. 
The choice J = 4 is also possible for levels excited from higher-lying states 
of spin 2 or 3. 

The coefficients 1, can be evaluated as the product of two terms (Bieden- 
harn and Rose 1953), each of which depends on only one transition in the 
cascade. 

A, = AM A™ wm [A (1—0,)] x [42 (1-—0,)] 
where Q is the fractional intensity of the higher-order multipole in the tran- 
sition. Q is equal to 6?/(1+6*) where 6 is the ratio of the reduced matrix 
element for the higher multipole to that for the lower multipole in the same 


I 


transition. In all of the cascades of concern here Q2 0, AS? = —0.5976, and 
A‘ = —1.069 so that the theoretical A, values are given directly in terms 
of AS’ and Q;. Arns and Wiedenbeck (1958) have published a convenient 
table of graphs showing .1$" and A? as a function of Q for all cases of interest. 
A comparison of experiment with theory is made by dividing the experimental 
Ay and Aq coefficients by the values of LS? and AY? given above and then 


examining these graphs to find a pair of A$” and 1S values tor a given spin 


sequence consistent with these experimental values and a common Q value 
The lifetime of the 0.155-Mev level has been found to be 6.5 X LO~’ seconds 
(Sunyar 1954). This lifetime is long enough that one might expect to find 


some attenuation of the correlation through nuclear reorientation in the 
intermediate state. The effect of this attenuation is to replace the .l, co- 
efhcients by A,G, where the values of G, (v = 2,4...) are the same tor all 
cascades involving this intermediate state. Steffen (1955) points out that the 
attenuation is a function only of the lifetime of this state and the environment 
in which the nucleus finds itself. While one might expect some attenuation 
effects for cascades involving the 0.155-Mev level, none should appear tor 
those involving the much shorter-lived 0.633-Mev level. 

The 931-155 coincidence data must be corrected for contributions from 


gamma rays of energy greater than | Mev which tend to smear out the 
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observed correlation pattern. Even without this correction, it is obvious that 
this cascade can only be associated with a 0-2-0 spin sequence since none of 
the other possibilities can give such a large A, coefficient (Fig. 5, (A), and 
Table I). 


478-155 
(8) 


931-155 


(A) 


120 150 180 
8 


5. Comparison of experimental and theoretical angular correlation functions in 
The dotted curves present the experimental data corrected for the finite solid angle 
detectors. The solid curves are the “pure’’ experimental angular correlation functions 
orrections for finite solid angle and for interfering cascades have been made, The broken 
s are theoretical functions 

931-155. The “pure” experimental function and the theoretical function for a 0-2-0 

can be matched exactly with attenuation factors Gy = 0.87 and G4 = 0.56 

fhe “pure” experimental function when corrected for attenuation using the factors 
above yields the heavy broken curve. It agrees closely with the theoretical curve shown 
2-2-4) cascade with VY; = 99.7%, 6 ve 

1142 644. The theoretical function is for a 0-2-0 cascade 

672 63%. The theoretical curve is a 2-2-9 correlation function with Q; = 0.03, 6 —ve 
828 644. The theoretical curve is a 2 2-0 correlation function with Q; = 0.01, 6 ~ ve 
}308 633. The theoretical curve for a 2-2-0 cascade with VQ, = 0.018 and 6 ve fits 


xperimental dati exactly 
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value of A'Y and A‘) as a function of Q for the cascades of Fig. 5 

(A) The spin sequence 1-2-0 

(B) The sequence 3-2-0 

(C) The sequence 2-2-0 applied to the 478-155 cascade 

(1)) The sequence 2-2-0 applied to the 1808-633 cascade 

(f}) The sequence 2-2-0 applied to the 828-633 cascade 


(BF) The sequence 2-2-0 applied to the 672-6335 cascade 





766 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


However, the “‘pure’’ experimental coefficients are still not as large as the 
theoretical 0-2-0 coefficients (42 = 0.358, A, = 1.142) and one is tempted 
to regard the difference as due to attenuation of the type discussed above. To 
bring the experimental coefficients in agreement with the theoretical ones 
requires that G2 = 0.87+0.16 and Gs, = 0.56+0.06. The validity of this 
assumption can be tested by applying these factors to the 478-155 cascade 
since it should be attenuated in the same way. 

It is estimated that 10% of the observed 478-155 coincidence recorded 
may be due to the effect of higher-energy gamma rays in coincidence with 
the 0.155-Mev radiation. The pure correlation function which is presented 
in Fig. 5, (B), and whose coefficients are tabulated in Table I corresponds 
to 1$? = +0.080+0.020 and A‘? = —0.152+0.015. An inspection of Fig. 6 
shows that these values rule out the cascades 1-2-0 and 3-2-0 and suggests 
that the cascade may be an attenuated 2-2-0). The A$” coefficient is con- 
sistent with Q; = 0.15 or 1.0 while the A‘” coefficient yields Q; = 0.46. 
Application of the attenuation factors Gz and G,; deduced from the 931-155 
cascade alters 4§ and A‘ to +0.092+0.03 and —0.27+0.03 respectively. 
While the change in A$” makes no significant change in the predicted Q 
values the attenuation-corrected A{” coefficient is consistent with Q; = 0.88 
+0.12, thus removing the inconsistency present with the uncorrected co- 
efficients. Figure 5, (B), presents both the correlation function for the attenua- 
tion-corrected coefficients and the theoretical 2-2-0 function for Q = 0.997 
and 6 +ve. The agreement is excellent and leads to the conclusion that the 
intensity of the dipole component in the 0.478-Mev transition is less than 
0.5%. Support for this conclusion is obtained from the measured internal 
conversion coefficient for this transition which sets an upper limit of 16° on 
the intensity of the dipole component. 

The 478-155 cascade has also been studied by Potnis et al. (1956). Although 
these workers reached the same conclusions concerning the 0.478-Mev 
radiation as those above, they obtained the unattenuated correlation function 
with solid metallic sources. The detailed shape of the 478-155 experimental 
function is very sensitive to the amount of the 931-155 cascade accepted by 
the apparatus because of the strong anisotropy of the latter cascade. In the 
present experiments, the “‘window”’ used for the 0.478-Mev gamma ray accepted 
only about half of its photopeak. Even with this window, it is estimated that 
10% of the coincidences recorded were due to higher-energy gamma rays. 
With the wider window used by Potnis and his co-workers, it is conceivable 
that 30% of the coincidences were due to these radiations. Such a contribution 
is quite sufficient to explain the discrepancies between their experimental 
curves and ours. It is possible that the full correlation obtained by Potnis 
et al. was fortuitous, and that their correlation function would have been in 
close agreement with ours had they corrected for the coincidences caused by 
higher-energy gamma rays. On the other hand, it is also quite reasonable to 
attribute the differences between their work and the results found here to 
the differences in source environment in the two experiments. 

The 1308-633 cascade requires no background correction so that its inter- 
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pretation is relatively simple. The data leads to values of —0.30+0.01 and 
—0.016+0.005 for A$? and A$” respectively. The 2-2-0 spin sequence is the 
only physically acceptable cascade which can simultaneously give a large 
negative A$” and a small negative A$” coefficient (Fig. 5, (F), and 6, (D)). 
Thus, according to the data, the 1.308-Mev transition must be a dipole— 
quadrupole mixture with Q; ~ 0.02 and 6 —ve and the spin of the 1.941-Mev 
level is 2. By arbitrarily increasing the limits of error on A§” to 0.03, the data 
can be fitted to a 3-2-0 correlation function with Q; = 0.12 and 6 —ve. While 
the spin choice of 3 for the 1.941-Mev level cannot be rigorously excluded 
on the basis of the angular correlation data, it seems to be effectively ruled 
out by the fact that the log fot value (8.2) for beta transitions to this level 
is appropriate to a first-forbidden non-unique beta decay and considerably 
too small for a first-forbidden transition of the unique type. Thus it is con- 
cluded that the 1.941-Mev level has spin 2 and even parity. 

The large A, coefficient associated with the 1132-633 experimental corre- 
lation function suggests that one of the cascades contributing to this data is 
0-2-0 in character. Since the 1308-633 cascade is definitely not of this type, 
the 1132-633 cascade must be. Subtraction of the 1308-633 cascade back- 
ground from the experimental data leads to the ‘‘pure”’ correlation coefficients 
of Table I and the solid curve of Fig. 5, (C). Although these do not represent 
the complete anisotropy expected for a pure 0-2-0 cascade the differences 
can easily be due to uncertainties in the magnitude of the background from 
the 1308-633 cascade. It is concluded that the level at 1.765 Mev has zero 
spin. Independent evidence for the existence of a zero-spin level above 1.3 
Mev is afforded by the strongly anisotropic correlation function obtained 
for the coincidences between all gamma rays above 1 Mev and the 0.155- 
Mev transition. 

The 828-633 experimental correlation function has been corrected for the 
background of coincidences due to two higher-energy cascades as shown in 
Table I and Fig. 5, (E). While the “pure” 4.2 coefficient is quite insensitive 
to the magnitude of this background, the “pure” A, coefficient depends on 
it strongly and can be driven negative if the intensity of this background 
is made large enough. The “pure” coefficients of Table I are consistent only 
with a 2-2-0 cascade in which Q, is less than 0.02 and 6 is negative (Fig. 6, 
(E)). However, if the estimated background intensity is increased from the 
25% upper limit used to the rather implausible value of 33%, the data can 
also be interpreted in terms of a 3-2-0 cascade with Q ~ 0.2 (6 —ve). The 
spin 3 assignment for the 1.461-Mev level is inconsistent with other evidence 
available. The log fot value for the beta transitions to this level has been 
established as 8.4 from the gamma-ray intensity balance. While failure to 
detect weak gamma rays makes this value a lower limit, it is inconceivable 
that the careful search for weak gamma rays in the energy region 0.2 to 
0.7 Mev carried on earlier in this laboratory (Johns et al. 1956) could have 
missed radiation of sufficient intensity to raise this value to more than 8.7. 
The value 8.4 lies between the well-established 8.6 and 8.1 log fot values for 
the first-forbidden non-unique transitions to the 0.155-Mev and ground 
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states of Os'** and is much too small for the first-forbidden unique transition 
required if the 1.461-Mev level is 3+ in character. Although the absence of 
the ground-state transition from this level seems to suggest a spin assignment 
greater than 2, the upper limit of 0.02% set on the intensity of this transition 
(Johns et a/. 1956) is quite inadequate to help one in deciding between the 
2+ and the 3+ choice. In the light of the angular correlation data and the 
fot value for the beta transitions to this level, it is concluded that the 1.461- 
Mev level has spin 2 and positive parity. 

One third of the observed 672-633 coincidences may be due to background 
from other cascades, and consequently the “‘pure’’ coefficients tabulated in 
Table I show rather large experimental errors. In fact when the errors are 
considered the coefficients (Fig. 6, (F)) are consistent with the following three 
possible cascades, in order of decreasing probability: 2-2-0 with Q,; = 0.03 
(6 —ve), 3-2-0 with Q, = 0.06 (6 —ve), and 4-2-0. 

It is impossible to distinguish between these three choices on the basis of 
the angular correlation measurements alone. If the spin 4 assignment for the 
1.306-Mev level is accepted, the beta transition feeding this level must be 
of negligible intensity and undetected gamma rays of a total intensity of 
0.09% must be present to feed this level. Since possible transitions from the 
1.461- or 1.941-Mev levels to this state are masked by other stronger radia- 
tions, it is impossible to test this assumption directly. However, in view of 
the small energy release involved in such radiations and the general intensity 
pattern of the decay process of Fig. 1, it seems very unlikely that radiations 
of sufficient intensity to feed this level in .09% of the disintegrations can 
occur. There seems to be no way of deciding between the assignments 2+ 
and 3+ from either the angular correlation data or the ft value for the 
beta transition to the 1.306-Mev level. Nor can the presence of a 1.806-Mev 
gamma ray be used to reach a decision, since such a radiation fits into the 
decay scheme in three places. The assignments actually made in Fig. 1 are 
based on coincidence results and on the known doublet line profile of this 
gamma ray as seen with the 50-cm beta-ray spectrometer (Johns et a/. 1956). 
No evidence concerning the existence of a 1.306-Mev ground-state transition 
exists at the present time. Therefore, one can only conclude that the spin 
of the 1.306-Mev level is 2 or 3 with positive parity with a reservation that 
the spin 4 assignment is remotely possible. 

In principle, it should be possible to use the angular correlation technique 
to determine the spin of the 1.958-Mev level by studying the 1802-155 
cascade. However, all cascades involving the 0.155-Mev radiation proved 
difficult to handle because very few of the 0.155-Mev photons are coincidence 
makers while their high intensity results in overloading in the side-channel 
analyzers unless weak sources are used. On the other hand, most of the 
0.633-Mev photons are potential coincidence makers and strong sources can 
be used when the 0.155-Mev photons are removed with lead absorbers. It is 
estimated that the maximum achievable coincidence rate for the 1802-155 
cascade would have been an order of magnitude smaller than for any of the 
other cascades studied. One can draw certain conclusions concerning the 
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spin of the 1.958-Mev level from the decay scheme of Fig. 1. The fact that 
the beta transition to this level has a log fot value of 8.0 characteristic of a 
first-forbidden transition leads to the conclusion that this state has positive 
parity. The fact that the gamma-ray transitions from this level to the zero- 
spin ground state and the 2+ state at 0.155 Mev are roughly equal in intensity 
restricts the spin assignment of the 1.958-Mev level to the choices 1 or 2. 
The 631-137 cascade of Os'** must be of the form J—2-0 with J restricted 
to the values 1 or 2 by virtue of the fact that the gamma-ray transitions 
from the 0.768-Mev level to the ground state and the 2+ state at 0.137 Mev 
are roughly equal in intensity. Since the measured A, coefficients for both 
solid and liquid sources are positive while the A 4 coefficient for a 1-2-0 cascade 
is zero or negative no matter what mixing occurs, the data of Table I and 
Figs. 4 and 6 shows clearly that the angular correlation pattern must be inter- 
preted as a 2-2-0 cascade. The fact that the coefficients for the liquid and 
solid sources are different suggests that the attenuation due to nuclear re- 
orientation in the intermediate state occurs with one or both types of sources. 
The experimental solid-source data yields A$” = +0.044 which gives Q = 0.12 


(6 —ve) or Q = 0.99 (6 +ve) and A$” = —0.057 consistent with Q = 0.18. 
The experimental liquid-source data yields A$” = +0.094 associated with 
Q = 0.15 (6 —ve) or Q = 0.995 (6 +ve) and A$” = —0.115 consistent with 


Q = 0.36. Steffen (1955) points out that attenuation factors are not smaller 
than the “hard core’’ value of 0.2 for Gz and 0.11 for Gs. Therefore the 4 §” 
and A‘” coefficients corrected for attenuation (i.e., A4$?/G2 and A‘$”/G4) must 
lie between the experimental values for the solid curve and +0.22 and —0.52 
respectively. The range of values allowed for A$” limits Q to values from 
0.16 to 0.26 (6 —ve) or 0.99 to 1.00 (6 +ve or —ve) while the range for 4‘) 
limits Q to values between 0.36 and 1.00. 

The only value of Q consistent with both the 43” and A$” ranges is 
0.99+0.01 (6 +ve or —ve). This choice imposes essentially no limits on the 
G» factor and requires G, to be 0.19 and 0.38 for the solid and liquid sources 
respectively. On this basis, it is concluded that the spin of the 0.768-Mev 
level is certainly 2+, and that the 0.631-Mev transition is about 99% E2. 

No internal conversion data exists to check these conclusions but independ- 
ent angular correlation data has been provided by Hurley and Jastram (1954). 
These workers, using electrochemically deposited sources (Jastram 1958), 
found an unattenuated 2-2-0 angular correlation pattern for this cascade 
corresponding to Q = 0.99 with 6 —ve. 
INTERPRETATION 

Figure 7 presents the information available concerning low-lying levels in 
four even-even osmium nuclides together with additional details concerning 
the beta-ray and gamma-ray transitions involved in the decay of Re'*® and 
Re!88, The following observations concerning the level systematics suggest that 
these nuclei fall in the transition region between spheroidal nuclei with well- 
developed rotational bands and nearly spherical nuclei with a vibrational 
level structure. 
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Fic. 7. Excited states in Os!8®, Os!88, Os!%, and Os!*. The decay schemes for Re'’* and 
Re!88 show the positions of the levels, the log ft values of the beta groups feeding them, and 
the multipole classifications of the gamma rays wherever they are known. The decay scheme 
of Re!8§ is drawn to illustrate the interpretation in the discussion. The data for the high 
spin states of Os!8* and Os!88 and for the levels of Os!® and Os!® is taken from the refer- 


ences discussed in the introduction. 


(a) The spacings of the 0, 2,4,6... levels of the ground-state rotational 
bands show increasing departures from those expected for pure rotational 
levels as the mass number increases; thus, while in Os!8* the 4+ and 6+ 
states are respectively 4% and 10% too low, the departures amount to 12% 
and 20% in Os!®. 

(b) The excitation energy of the second 2+ state (hereafter referred to 
as 2’+) which is observed in all four nuclei decreases rapidly with increasing 
mass number and has almost reached the position expected for a two-phonon 
vibrational state when the mass reaches 192. 

(c) The ratio of the ft values for the beta transitions from Re'®* to the 
two lowest states of Os'** is 1.6+0.4 in good agreement with the value 2.0 
expected for a pure rotational structure. However, in the Re'®® decay, the 
corresponding ratio, 3.20.6, already shows the effect of departures from pure 
rotational behavior. 

(d) The ratio of the reduced E2 transition probabilities for the 2’ — 2 
and 2’ — 0 transitions is found to be 2.70.6 and 2.6+0.3 in Os!®* and Os!88 
respectively. These ratios are to be compared with a value of 1.43 predicted 
for a pure rotational structure and a value of infinity for a pure vibration. All 
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these facts point to the conclusion that we are dealing here with nuclei, in 
which a strong coupling exists between the rotational and vibrational modes 
of collective motion. i 

The 2+’ level in Os'** has been interpreted as the (2, 2)-+ state associated 
with a gamma vibration (Alder et al. 1956; Nielsen and Nielsen 1958). This 
interpretation is supported by comments (b) and (d) above and by the 
predominantly #2 character of the 2’ — 2 transition. 

An alternative description of the 0, 2, 4, 2’ levels in all four nuclei is given 
by Davydov and Filippov (1958), who have investigated the energy levels 
to be expected from the rotation of an asymmetric top. They show that the 
rotational levels, instead of being a simple 0, 2,4... sequence, now include 
the spin values 0, 2, 2, 3, 4, 4, 4, 5,5... Thespacing of these levels is evaluated 
in terms of the parameter y, which determines the deviation of the nucleus 
from axial symmetry. In particular, the authors express the ratio of the energies 
of the 2’ and 2 states and various reduced transition probabilities in terms 
of the parameter y. Table II presents a comparison between the results of 
their theory and experiment for the four osmium nuclides. In this table y is 


TABLE II 


[b( 2:2’ — 2)]/ 


[b(E2:2’ — 0)]/ 





[b(E2:2’ > 0)] [b(£2:2 — 0)] 
E(2’)/E(2) 7 Theor. Expt. Theor. Expt. 
Os!86 5.6 16.6° 3.2 2.7 0.062 
Os!88 4.1 19.2° 4.7 2.6 0.072 0.16 
Os!99 3.0 22° 8.2 0.067 0.047 
Os! 6 24° >12 0.055 0.059 


| no 


evaluated from the ratio of the excitation energies of the two spin 2 levels 
and the values so obtained used to obtain the other ratios in the table. The 
experimental results in the last column are obtained from Coulomb excitation 
data (Barloutaud et al. 1958). 

In a spherical nucleus with its pure harmonic structure, the one-phonon 
state has spin 2, the two-phonon state is a degenerate triplet with three 
possible spins, 0, 2, or 4, and so on. The addition of nucleons to an unfilled 
shell introduces interactions between these particles and the collective motions 
which remove these degeneracies and gradually alter the level spacings so 
that the description in terms of phonon excitation becomes inappropriate. For 
these nuclei, it is useful to consider the collective motions in terms of rotational 
and vibrational levels, each with their characteristic form of level spacing. 
Since the strength of the coupling between collective and particle motions 
does not change the number of collective levels but merely alters their spacing, 
one can still use the phonon label on a level in a spheroidal nucleus to indicate 
the level with which it is associated in a spherical nucleus. This labelling is 
particularly useful in transition nuclei where the level structure reveals 
qualities of both spheroidal and spherical nuclei. Thus, in even-even nuclei, 
the first member of the ground-state rotational band (K = 0 J = 2) can be 
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labelled a one-phonon excitation, while the K = 0 J = 4 member of the 
same band, the K = 2 J = 2 (gamma vibration), and the K = 0 J=0 
(beta vibration) can all be labelled as two-phonon excitations. 

In the language of the last paragraph, the 2+ state in all four nuclei of 
Fig. 7 would be called a one-phonon state while the 2’+ and 4+ levels would 
both be labelled two-phonon states. It is attractive to consider the 0+ state 
of Os!*§ at 1.086 Mev as the third member of the two-phonon triplet and 
label it as a two-phonon beta vibration. The level of: excitation of this beta 
vibration should decrease with increasing mass number much as the 2’+ 
gamma vibration does and the expected trend has been indicated in Fig. 7. 
Unfortunately, none of the processes studied in the other nuclides of this 
series would be expected to populate such a 0+ level with appreciable intensity. 
It is, however, interesting to note that a 0+ state at about the same excitation 
has been identified in Er!®* (Fraser and Milton 1955), which is located in the 
transition region just above the 82-neutron closed shell. 

If the 0+ state at 1.086 Mev is accepted as the candidate for the two- 
phonon beta vibration, the 1.306-Mev level is probably the first-excited state 
of a rotational band based on this vibrational level. This assignment for the 
1.306-Mev level agrees with the preferred experimental spin value and leads 
to a reasonable value for the moment of inertia of the nucleus in the ground 
state. The ft values for beta transitions to these two levels have a ratio of 
1.6, in agreement with the value 2.0 expected for a pure rotational structure. 
However, one would expect the transition from the 0, 2+ level of this rotational 
band to the 2, 2+ level, associated with the two-phonon gamma vibration 
to be predominantly E2 in character. In fact, the observed transition is 
predominantly M1 (whether the angular correlation data be interpreted in 
terms of spin 2 or 3 for the 1.306-Mev level). This relaxation of the K-selection 
rule is another indication that the collective model represents an oversimpli- 
fication of the situation in this nuclide. An alternative choice for the second 
member of the rotational band based on the beta-vibrational level is the 2+ 
state at 1.461 Mev. This choice has nothing to recommend it since it leads 
to a very small moment of inertia of the nucleus. Moreover, the ratio of the ft 
values for beta transitions to the 1.461- and 1.086-Mev levels is only 0.16 and 
the 0.828-Mev transition de-exciting the 1.461-Mev level is predominantly M1. 

It is impossible to be definite about the interpretation of the four remaining 
levels of the decay scheme of Fig. 7 since the excitation level is now high 
enough that particle excitations are likely to be in competition with collective 
motions. It does seem likely that the levels at 1.765 and 1.941 Mev are the 
0, 0+ and 0, 2+ members of a third rotational band since these levels 
show much the same pattern of de-excitation as that from the pair of 
levels at 1.086 and 1.306 Mev. Support for this view comes from the value 
of ~ 1 for the ratio of the ft values for beta transitions feeding this pair of 
levels. On the other hand, the de-excitation of the 1.941-Mev level is pre- 
dominantly by M1 radiation, which is again K-forbidden. 

One striking feature of the decay scheme of Fig. 1 or Fig. 7 is the absence 
of ground-state transitions and the intensity of transitions feeding the 2’+ 
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state at 0.633 Mev. There are four levels in the decay scheme which feed 
both the 2’+ and the 2+ levels of Fig. 7 and in every case, the reduced 
transition probability for the transition to the 2’+ level is from 3 to 40 
times greater than that for the transition from the same level to the 2+ 
state. The reduced transition probabilities for the unobserved ground-state 
transitions are of course smaller still. This indicates that the internal structure 
of the nucleus in the upper levels is more like that in the 2+’ level than in 
the ground or 2+ states. 

Wilets and Jean (1956) and Mottelson (1957) have pointed out that the 
transition nuclei are “‘softer’’ against shape deformation than either the 
truly spherical nuclei or the strongly deformed spheroidal nuclei. The observa- 
tions of the last paragraph would suggest that Os'** is such a “‘soft’’ nucleus 
and that a gamma vibration of energy 0.633 Mev is sufficient to bring the 
nucleus into a new equilibrium shape which it retains for many of the higher 
excitations. 

This view seems to deny the possibility of regarding the Os!** nucleus as 
an asymmetric rigid rotator since the basic premise of Davydov and Filippov’s 
treatment is that the asymmetry is built into the nuclear shape in the ground 
state and is unaltered by the excitation of the rotational levels. Of necessity, 
their treatment can say nothing about the location of the 0+ beta vibration 
since such a vibration by definition is a shape distortion which does not 
affect the axial symmetry of the nucleus. Despite these conclusions, Davydov 
and Filippov’s theory is very attractive since it leads for the first time to 
detailed numerical predictions concerning level spacings and _ transition 
probabilities for these transition nuclei. In contrast, the approach of the 
Bohr-—Mottelson theory through beta and gamma vibrations leads only to 
order of magnitude estimates—what might be facetiously called a ‘‘soft 
theory” for ‘soft nuclei’. 
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ON THE CALCULATION OF THE VELOCITY OF SOUND IN SEA WATER 


CHARLOTTE FROESE* 






INTRODUCTION 

An accurate knowledge of the velocity of sound in sea water as a function 
of salinity, temperature, and pressure is essential in oceanography. Tables are 
available but the values quoted may differ by as much as 3 m/sec from one 
table to the next, when the desired accuracy is to within .1 m/sec. Of these 
tables, the two most frequently used are by Kuwahara (1939) and by Del 
Grosso (1952). 

Kuwahara computed the sound velocity, Vs, as a function of salinity, 
temperature, and pressure. His calculations were based on theoretical con- 
siderations which express the sound velocity in terms of specific heat, specific 
volume, and compressibility of sea water. Empirical formulae were derived 
for these physical properties, and so the accuracy of his tables will depend 














mainly on the accuracy of these formulae. 

More recently Del Grosso (1952) measured the velocity of sound at atmos- 
pheric pressure experimentally using modern interferometric techniques. From 
these results he derived an equation for V,. which he believes to be accurate 






to within .2 m/sec. 

Sound velocities at pressures up to 1000 atmospheres are currently being 
measured. These direct measurements are capable of greater precision than 
the indirect method of Kuwahara, and at the moment a considerable dis- 
crepancy exists between the two. For example, for V35,0,0 Del Grosso obtained 
a value 3.1 m/sec higher than the one given by Kuwahara. This discrepancy 
is an indication of the error present in the empirical formulae used by Kuwa- 
hara. Some of them were based on measurements made as far back as 1889 
and could be determined now with greater accuracy. 

The purpose of this investigation was, first of all, to check Kuwahara’s 
tables by using the same basic formulae but employing a different method. A 
program was written for ALWAC III-E, the automatic computer at the 
University of British Columbia, to perform the calculations. This program 
would also serve another purpose; as new measurements are made of the 
specific heat, specific volume, or compressibility, the same program may be 
used to recompute the sound velocity and compare the results with direct 
measurements. This paper deals entirely with the results of the first stage of 
the investigation. A description is given of the method by which the sound 
velocity was computed and the results compared with Kuwahara’s tables. 
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SOUND VELOCITY COMPUTATION 


The velocity of sound in sea water is given by the equation 


(1) 
where 
eg (=) 
ap DG xo 
a = specific volume, 
C, = specific heat at constant pressure, 


and 7 = temperature in degrees absolute. 


The specific volume, as,,, Was computed using the same formulae as those 
employed by Kuwahara (1939) but they were rewritten in a form more 
suitable for computation on an automatic computer. 

The specific heat of sea water was divided into two parts so that 


T (’ da 
C, = C,(st0) 5 . OF dp 
where the first is the value at zero pressure, and the second is the pressure 
correction term obtained from thermodynamic considerations. Kuwahara’s 
formula was used to compute the first but none was given for the second. In 
order to obtain the pressure correction term, the second derivative 0?a/0d¢? 
had to be evaluated. Though a is given by an analytic formula, it is too com- 
plex to warrant analytic differentiation; instead the derivative was computed 
numerically according to the formula 


(3) On = (Cin4.1— 2a +Oy-1) /h? +0(h*) 


with h = 1/2° C. A formula was then fitted to the results using a least squares 
approximation; values from the formula differed from the original data by 
less than 2 in 1000. From this formula the pressure correction was obtained. 

According to equations (1) and (2), values of a, (da/dp), and (da/dt) are 
required in order to compute the sound velocity. As before, the partial deriva- 
tives were replaced by differences using the formula 


(4) a, = (an41—On—1) /2h+0(h”), 


where h: was 100 decibars for the derivative with respect to pressure and 1/2°C 
for the derivative with respect to temperature. These partial derivatives occur 
in equation (2). The most important contribution is from the da/dp term and 
if V is to be correct to six figures, da/dp should be correct to about the same 
though some accuracy is gained in taking the square root. Therefore, h = Ap 
had to be chosen large enough so that the difference, aj41—Qp,—1, retained 
six significant figures, but small enough so that the truncation error in (4) 
did not affect these figures. A nearly optimum value was Ap = 100 decibars, 
which gave a difference of almost six figures and a truncation error of one in 
the sixth figure. The significance gained by increasing Ap would be offset by 
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an increase in the truncation error. The value of h = At is not as critical 
since the contribution to (2) from the second term is considerably smaller. 

Equations (1) and (2) were rewritten in terms of differences and the sound 
velocity computed for a series of values of s, ¢, and p. 


COMPARISON OF RESULTS 


The results of these calculations were compared with the values given by 
Kuwahara. In his tables, the sound velocity is divided into four parts so that 


V ssp = V35,0,.pt AV: +AV,+AV sip, 


and each of these quantities tabulated separately. In Table I, the velocities 
V35,0,p are compared; the differences are hardly more than those due to round- 
off, the largest being 0.07. The same applies to Tables II and III where the 
temperature correction, AV,, and salinity correction, AV;, respectively are 


TABLE I 


Comparison of the velocities, V35,0,p in m/sec 


p, decibars 


Kuwahara Difference 


ALWAC 


0 - 1445.44 1445.5 — .06 
500 1454.45 1454.5 — .05 
1000 1463 .43 1463.4 03 
2000 1481.30 1481.3 00 
3000 1499 .02 1499.0 02 
4000 1516.57 1516.6 — .03 
5000 1533 .92 1533.9 02 
6000 1551.05 1551.1 05 
7000 1567 .94 1568.0 06 
8000 1584.55 1584.6 .05 
9000 1600. 86 1600.9 04 
10000 1616.86 1616.9 04 


TABLE II 


the temperature correction, AV;, in 
m/sec 
ALWAC Kuwahara Difference 
—9.32 -—9. — .02 
0.00 0 .00 
8.99 gy .O1 
17.62 17.6 02 
21.79 Zi. .O1 
41.25 41.; 05 
58.34 58 06 
73.22 73. .02 
86.18 86. .02 
97.65 97. 05 


TABLE III 


Comparison of salinity correction, AV,, in m/sec 


Salinity, %@ ALWAC 

30 —6.54 —6.5 — .04 
33 —2.62 —2.6 — .02 
35 0.00 0.0 .00 
37 2.61 2.6 01 
40 6.53 6.5 03 


Kuwahara Difference 
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NOTES 








given. Greater differences occur in Table IV where AV,,,, the correction for 
simultaneous changes in any two or more variables, are compared. These 
usually occur at the limits of the range where ¢ is large. In Fig. 1, plots of 
AVs:) as a function of temperature for several different pressures and salinities 
are given both as computed by ALWAC and by Kuwahara. These indicate 
the relative magnitude of the differences. 


TABLE IV 


Comparison of AVstp = Vetp— Va5,0.»>p—A4Vi—AV, in m/sec 



























Salinity, %o 
30 


35 
p, deci- $$$ 


bars t,°C ALWAC K. Diff. ALWAC K. Diff. ALWAC K. Diff. 


40 




























0 —2 -.07 —-.1 .03 .00 .O .00 .09 1 —.01 
0 .00 0 .00 .00 .0 .00 .00 0 .00 
5 .21 2 01 .00 .0 00 —.27 —.2 —.07 
10 45 5 —.05 .00 0 00 —.57 —-.56 —.07 
15 .70 8 —.10 .00 0 .00 —.87 -.8 —.07 
20 .92 1.0 —.08 .00 0 00 -1.14 -1.0 —.14 
25 1.09 12 —:11 .00 0 .00 -1.37 -—1.2 —.17 
30 1.19 14 —.21 .00 0 00 -1.53 -1.3 —.23 
2000 0 -—.24 —-.2 —.04 .00 .0 .00 24 2 -04 
5 — .02 OO —.02 —.08 .0 —.038 —.10 —i .00 
10 .08 1. —.02 —.22 —-.2 —.02 —.64 —.6 —.04 
15 18 23 —.12 —.4l1 —-.4 —.01 -—1.18 -1.1 —.08 
20 41 6 —-0 —.4 —.5 06 —-1.55 —-1.4 —.15 
25 .88 10 —.12 —.18 —.2 02 —1.57 —1.4 —.17 
4000 0 —.32 —-.3 —.02 .00 .O .00 .32 3 .02 
5 —-.40 —.4 .00 —.36 —.4 04 —.39 —.4 .O1 
10 -.638 -—-6 -—.08 —.90 -1.0 10 -—1.32 -1.3 —.02 
15 —-.70 —.7 .00 -—1.30 -—1.4 10 -—2.11 —2.1 —.01 


| 




























6000 —2 .09 1 —.01 .23 2 .03 .37 3 .07 
0 -.09 -.1 01 .00 0 .00 .09 1 —.01 
2 —.37 —.4 .03 —.33 —-.3 —03 —.29 —.3 01 
4 -.70 —.7 0 —.71 -—-.7 —01 —.7% —-.7 —.03 
8000 —2 .52 5 02 54 5 04 56 5 .06 
0 .02 0 .02 .00 0 .00 —.02 0 —.02 
2 -.56 -—-.5 -—.06 —.62 -.6 -—.02 —.70 —.7 .00 
4 -1.16 —1.1 —.06 —1.28 -1.2 —.08 -1.40 -1.4 .00 














.08 Re .02 91 9 .O1 74 ad .04 
0 21 2 -O1 .00 0 00 —.21 —-.2 —.01 
2 -.70 —7 oO —-.9 -—1.0 .05 -1.20 -1.3 .10 
4 -1.58 —1.6 .02 -1.87 -—1.9 03 —2.17 —2.2 .03 









The author is grateful to Dr. N. P. Fofonoff for many helpful discussions 
and suggestions. 
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A SPHERICAL WAVE EXPANSION FOR DOUBLE $-DECAY* 
S. P. ROosEN 


In a recent article in this journal, the spherical wave solutions of the Dirac 
equation for a central field and the properties of angular momentum coupling 
coefficients were used to derive a general formulation of the theory of single 
8-decay transitions of arbitrary forbiddenness (Lee-Whiting 1958). The theory 
of double 8-decay can be formulated in a similar manner and it may be of 
interest to show briefly how this is done. 

Consider the general nucleon—lepton interaction Hamiltonian: 


(1) H= DO (v3O.n) {C.[V2O.(1 +6275) Wo] +Dz[v20.(1+n275) vs] } +h.c. 


It has been shown (Primakoff 1952; Rosen 1957; Primakoff and Rosen 1959) 
that the matrix element 4(i—/) for the neutrinoless double 8-decay of the 
nucleus (A,Z) to the nucleus (4,Z+2) 


(A,Z+1)+e.4+9 
a 


wv \ 
(2) (4,2)< “ (A,Z+2)+e:+e: 
X 
via D via C 
f 


M 
(1,Z+1)+ei:tv 


is given by 


2\*( me? ‘ 
(3) <a A(t— f) 


ss (way y ct, ag wD TETj}OpO jp 
x[ i (eo 2i( tH) 7 nap EH | 
: "kj Vrs lrj 


x C(Uhy+ Jrv¥s) ower, | W dr 


*This research was supported in part by the United States Air Force through the Air Force 
Office of Scientific Research and Development Command, under contract No. AF 18(603)-108. 
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where mc? H; (t = 1,2) are the electron energies and C is the charge conju- 
gation operator. The forms used for the invariants O, are 
Os = B, Or = (86, Ba), Op = BYs, 


(4) ; 
Oy = (1, 7a), On = (6, t¥s) 


and the indices X, v of (3) refer to the spherical components ot these operators.* 
If Oy is a scalar-type operator (8, ys, 1, etc.), then there is only one term in 
the summation over A, and the factor (—+/3C°}) is replaced by unity; if 
Oy is vector type (6, a, etc.), then no such modifications are needed (simi- 
larly for O,). Also 

Thy = Q\D,(1 —6xn,) +C,D\ (1 —6m), 

Jyv = QD,(6.—1+) + CLD (6, —m) ; 


W;, WV, are the wave functions of the initial and final nuclei respectively; 77, 
tt are the usual isotopic spin operators; r,, f; are the position vectors of the 


(9) 


k®, j™ nucleons; and rm, = f—Ty). 

The spherical wave function y,(r) for the 7™ 
as (Spiers and Blin-Stoyle 1952) 

2) mi \ dili 4 Bia 7mi—o 
(6) ¥i(r) = Vis (r, Ti, Bi) = C mi-—Oi,o% I li ‘(pr) Fe "(6, d) 

(with J; = j;+}a;8;), 
corresponding to an electron with total angular momentum j;, m; and linear 
momentum p. Both parameters a;, 8; take the values +1, and o,, which 
corresponds, in the non-relativistic limit, to the z-component of spin, is +}. 

There are two types of operators acting on the nuclear states in (3): 


(71) OpO jr 4 Yo(ry )OXC[ hot Jrvvs]Or Wilt) 
(712) OpnO jr Po(¥, Oa . PjC[Tyy+ JyvyslOr (85) * 


electron (2 = 1, 2) is written 


(Note that XK means complex conjugate, ~ transpose, and + hermitian 
conjugate.) In both expressions (7), the wave functions (6) are used to form 
eigenfunctions of the total angular momentum J, M/ of the electron pair. 
The operators defined by 

(8) Ey, = (O_,C{ Ih» +JSnvys]O_,), 

Dyys = (O_ya_yC[ Ih» +-Jrvvs5]O_») 

can be written as the outer products of two (2X2) matrices 


(9) (Exy) = (Exv)arat = (Exv)?X (Env) ares 
; (Dyyv) = (Dyuv) 7?" XK (Dyyv) 8081 
where each matrix acts on the variables indicated and both (£y,)8.8,, (Dawv)a.a, 
are independent of Auv (see Spiers and Blin-Stoyle (1952a@) for details). 
1 
"A= —,(detidy), Ao= Ay Aa = (de idy) 


A.B 


(—V3) > C% 4 AyB_y. 
A 
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From (6), (8), and (9), the lepton factor of (72) is equal to 


J i fe fila Jj2l2 
(10) eo Cu mi,m2 Cc O cente Cc Oates 
BiB2 
o102 
mi+m2=M 


x[ VG) VaR) Fi (bars) Fis (bere) (EX) PS 
where 
l,= jit 24;B; 
and 


Y7(j) = Y7T(y $5). 


By means of the relation (Spiers and Blin-Stoyle 1952ca) 


(11) 28 Ce a “ Cc” a seat z Cc” os i er (j) yat-"*(p) 


mi+m2=M 


se i le JL S43 M—¢i—0o2 
os » Biiy(JLS) Cc (M—o1—02), pyiak Cc saat L:il 


and (10), the expression (77) can now be written as 


(12) Wier joax) (J) a » Wher s202)(JLS) = z. p> (Htj10181 20262) (JLS)]* 
‘ 182 


LS 


XK (410181 20282) us) 


where 


(13) Hjresb: saa p.)(JLS) 
oe 1) On EO Gen nr O ERE. 


0192 
A» 


(14) K (410181 j20262) (JLS) ae [Byia(JLS, un (Pars) Fii"*(porr) (EXS)p2811"% 
and (Blatt and Weisskopf 1952, Appendix A) 


(15) Dern = 2) C’mm: Yn) Yn'(e). 


mi+m2=M 


Now expand 
(16) (Ex) = DE ev) (—1)"C4,.21.82 V2). 


For O,, O, scalar-type operators* 


(17) er (Av) = —dbr,oby,0 


and for O,, O, vector-type operators 


(18) ér(dv) = (—1)"" 4/ a Ti Cra. 


*All scalar-type operators have the same dependence on i, a2; the differences between them 
are in their 82, 8; dependence. The same is true for vector-type operators. 
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Therefore, from (14), (16), (17), and (18), 


(19) Hij.a161 420262) (JLS) = iv/2(0,0;) Den wd7, 18,0 (scalar Oy, O,) 
and for vector-type O,, O, 


aes L+S—J+M 2/3 csi M-y , 
= 1(— 1) /2S+1 x Cc S eey). e C, » » OpO yn Dz: hij - 


For the vector interaction (O, = O; = 1) and the axial vector (OQ, = O; = 6) 
(note that source velocity terms have been ignored), 





(Ex») B26) = — (py) 628,Ln»- 


In a similar way, the expression (777) can be written as 


(20) RGjieje)(J) = » Rojas is 202 »(JPS) = be > [TGrer6s sra262) (JP'S) \ 


PS B Bo 


\ 
XN (410181 J2a282) (JPS) f 





where 
(21) TU,a18; 120262) (JP'S) = a (— 1) oer /2F+1 Csi M— Cla (M—u \aites 
—71—02), 
vies 
“Cr . . (Dy) 7*** (OnnO yn)* 
e je I l, h 
VAC A Ff |e De 
ie “Ge, ee 
jn & A oh 
+{~+ 1)°V/2h +1 oa her o} j2 y l, f rs oo L 
ry P | 


(g = fatjethth+J+P+S) 


and 





(22) N (j,0181 j20282) (JP S) = [/(2P+1)( )(2S+1) (2{14+1) (2je +1) 
Fit"\(pars) Fis (pore) (DXr)axai* 
The 12-7 symbol is defined in Jahn and Hope (1954). Expanding 


(Dyyr)"* = DD (in/2)d (Au) (—1)7CZ, 2, 3, 
- 


(23) dh (Awv) = i (—1)"*7 67, 05,4 (Oy, O, scalar) 


i(—1)" x INT, GC" 5 ee Os (O,, O, vector) 
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where 

D(1,G) = (1-G) ,/ ea (4/2)9 
DEG) ek 


From (23) and (21) 


(24) 


25) Thpneen(JPS) = 2 +1) W/ sj2cy re DE aw Br, Df 
X5s,157,p(O.0;)* 
(O,, O, scalar). 
(the coefficients 4, B are as in equation (21)) 
= —/6 = (—1) 8+? (24-1) y/ (2S+1)W(GSJIf:P)D(S, G) 
njG 
X {C7 Se, Canary C78» OnOplAre Decne +Bry Dean )}* 
(Ox, O, vector). 
For vector and axial vector interactions 
(26) (Dyyv) p28; = (pz) 828,Ln0- 
From equations (3), (7), (12), and (20) the operator appearing between 
WV, and W; in the matrix element 4 (i — f) is 
(27) V = 2d tet} Viger jze2)(J) 
where 


Vejen joa2) (J) a 2 ve a ea FS}; 


s 


ES 2i(H+Hh) 


ej 


(28) vi ik a2 ») ( Js) w= 


| a Mais yea (FPS) 


es =H) » W(5ra1 j202)(JLS). 
rj 
The probability per unit time for neutrinoless double 6-decay, accompanied 
by the emission of electrons with energies mc*H; (i = 1, 2) is 
(29) P(Hh, H2)dH dH = oe te > ify V8 drt es 2 6(Ho —H,—H:) 
JM 


J1@1 j2a2 


XdH dH, 


where p} = H?—1 and mc*H, is the energy released in the process. 

The degree of forbiddenness is determined by the power of 7,; appearing 
in the denominator of V: the higher this power, the lower the degree of 
forbiddenness (Primakoff 1952). For low Z, the radial functions F'" (pr) 
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behave roughly as (pr)' (where / = j+43a8). Therefore, from (20), (21), (22) 
and (12), (14) the following rules can be given to select terms with degree 
of forbiddenness f: 


(i) fox =~ Ribas es (IPS) tile: dbl fie 
kj 


(22) for $e Weise: nas)(JLS) take /,+/, = f—2-+m. 
kj 


In both (2) and (iz) if an operator @ or ys appears r times (e.g. Ox. = aa, 
= 1) then f must be replaced by (f—r). Notice that /;, /2 are always greater 
than, or equal to, zero. 
7M hlg 1 ] ‘li 1 M 
Téieise2)(JPS) hasfactors C*¢*oDP-na, Co oDp-v 2: 
therefore, h = /,+1, h’ = 1:41, and the parity is (—1)''*'2+!, Similarly the 
parity of HV.a:jea2)(JLS) is (—1)"*'*. Therefore the parity selection rules 
are: 
for (2) the parity change is (—1)/*", 
for (iz) the parity change is (—1)/*”. 
If a, y; appear 7 times, f is replaced by (f—r). 
The values of J, M are determined by the total angular momenta (J;, Mj), 
(J,, M,) of the initial and final nuclei. 
The matrix element for two neutrino double 6-decay has also been con- 
sidered from this point of view, but, since there are now four leptons involved, 
it seems too complicated to study in detail. 


Most of the above work is taken from chapter three of S. P. Rosen (1957). 
The author would like to thank Dr. R. J. Blin-Stoyle of the Clarendon Labora- 
tory, Oxford, for suggesting this problem and the Department of Scientific 
and Industrial Research of Great Britain for a maintenance grant. 
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